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Abstract 


An  improved  propulsive  nozzlc/cowl  conc^t  has  been  proposed  which  may  have  some 
potential  for  increasing  the  overall  performance  of  hypersonic  flight  vehicles.  The 
concq>t  consists  of  placing  a  gas  generator  into  the  cowl  to  supply  a  relatively  high 
pressure  boundary  condition  to  the  main  exhaust  plume.  This  computational  study 
investigated  two  nozzle/cowl  getxnetries;  an  experimentally  validated  nozzle/cowl 
configuration  evaluated  at  off-design  ctmditions  and  a  generic  hypersonic  jH-opulsive 
nozzle  evaluated  at  m<xe  realistic  on-design  conditions.  The  flowfields  were  analyzed 
lining  a  combination  of  Van  Leer  flux-vectcx'  splitting  (FVS)  and  Roe  flux-difference 
splitting  (FDS)  finite  volume  computational  algoridims  implemented  in  a  code 
developed  by  Wright  Labwatory  (WL/PIMQ.  The  two-dimensional  Navier-Stokes 
equations  were  solved  assuming  laminar  flow  and  a  perfect  gas  equation  of  state.  This 
investigatimi  highlighted  the  effect  of  the  gas  ^nerator  on  the  nozzle  wall  pressure 
distribution  and  its  effect  on  overall  nozzle  performance.  The  analysis  onphasized  the 
effect  of  gas  generator  mass  flow  and  deHectimi  angle  ^fects.  In  general,  the  gas 
generator  tended  to  show  a  larger  improvement  in  nozzle  performance  at  lower  Madi 
number  and  nozzle  pressure  ratio  (NPR)  off-design  conditions  where  the  flow  is 
signiricandy  underexpanded.  For  the  very  high  Mach,  high  NEK  on-design 
cmifiguratiaQ  analyzed,  the  overall  system  performance  increased  with  the  addition  of 
die  gas  generato*  but  the  nozzle  wall  pressure  recovery  actually  tended  to  decrease 
slighdy.  This  unexpected  decrease  occurred  because  the  gas  generator  exit  pressure 
was  relatively  low  compared  to  the  very  high  recompression  pressure  experienced 


XIV 


behind  the  blunt  cowl  at  these  high  speed  trajectory  points.  Simply  placing  a  nominal 
gas  generator  into  the  flowfield  temled  to  have  a  much  larger  effect  on  nozzle  wall 
pressure  dian  varying  its  mass  flow  even  though  nozzle  performance  did  tend  to 
increase  with  increased  mass  flow.  Deflecting  the  gas  generator  flow  towards  the 
nozzle  wall  provided  a  dramatic  improvement  in  pressure  recovery  and  only  a  small 
penalty  was  paid  f<v  deflecting  it  away  from  the  wall.  For  the  e:q)erimentally 
validated  off-design  cases  analyzed,  the  gas  generator  successfully  acted  as  an 
aerodynamic  cowl  extension.  Gas  gentfatw  ^ects  on  nozzle  wall  pressure  recovery 
were  very  similar  to  cowl  geometry  extension  and  deflection  effects  previously 
analyzed  for  this  configuration.  Nozzle  operating  condition  had  a  significant  impact 
(m  the  gas  ^m-atm*  performance.  As  the  nozzle  exit  pressure  and  trajectory  Mach 
number  increased,  the  flow  became  increasingly  dmninated  by  the  large  initial 
e3q>ansions.  The  increased  initial  expansicm  combined  with  highar  internal  nozzle  flow 
Madi  numbers  decreased  the  overall  effects  of  the  gas  generator  and  tended  to  push 
the  propagation  of  the  gas  generator  effects  to  the  nozzle  wall  fiirtlmr  downstream. 
Grid  rdinement  effects  were  limited  to  the  regitm  near  the  peak  of  the  initial 
recmnfM’ession  where  the  finer  grids  captured  slightly  higher  peak  pressures. 
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COMPUTATIONAL  INYESTIGA'ilON  OF  AN  IMPROVED  COWL  CONCEPT 
FOR  HYPERSONIC  PROPULSIVE  NOZZLES 


I.  INTRODUCTION 


1.1  Badtground 

A  fundamental  desire  to  design  aircraft  that  fly  higher  and  faster  has  been 
evident  in  flight  vehicle  design  since  the  earliest  days  of  manned  flight.  The  ability  to 
develt^  and  build  manned  hypersonic  flight  vehides  was  proven  with  experimental 
flight  vehicles  such  as  the  X-IS  whidi  exceeded  Mach  7  and  350,000  ft.  in  1963  (3). 

A  series  of  manned  hypersonic  e]q)erifflental  aircraft  developed  in  the  1960*8  led  to  the 
development  of  the  first  reusable  manned  hypersonic  vehicle~the  Space  Shuttle,  which 
regularly  achieves  Mach  25  flight  during  its  reentry  glide.  The  next  major  stq>  is  to 
design  a  reusable  manned  single  stage  to  orbit  (SSTO)  flight  vehicle  such  as  the 
{xroposed  Natimud  AeroSpace  Plane  (NASP).  Unlike  the  X-'15  which  was  rodcet 
powered  or  die  Space  Shuttle  which  is  ui^iowered,  a  reusable  SSTO  vehide  must  have 
an  integrated  jx-opulsion  system  that  is  cqiable  of  accelerating  the  vehide  to  orbital 
speeds  in  excess  Madi  25.  These  hypersonic  flight  vehide  designs  will  be 
dominated  by  the  air-breathing  supersonic  combustion  ramjet  (SCRAMjet)  engines 
which  will  be  integrated  into  a  single  airframe-engine  prq>ulsion  system  (3). 

Hgure  1-1  shows  a  typical  hypersonic  vehicle  with  its  airframe-integrated 
propulsicm  system.  The  vehicle  forebody  behaves  as  a  compression  ramp  far  the 
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engine  inlet  flow.  Further  c<Mnpresai<m  of  the  flow  and  combusdon  take  place  in  the 
modular  chambers  below  the  vehicle.  The  flow  is  then  e;q>anded  along  the  vehicle 
afterbody  which  acts  like  a  nozzle  wail. 

At  hyptfsonic  speeds,  the  integrated  [vopulsion  system  is  designed  to  operate  at 
very  high  {xessure  ratios.  The  nozzle  pressure  ratio  (NPR)  is  defined  as  the 
SCRAMjet  total  exit  pressure  divided  by  the  ambient  fieestream  static  pressure, 
NPRsPi^m.  For  hypersonic  flight  ctmditions,  the  combustor  exit  flow  is  typically 
undere]q>anded  and  some  of  the  excess  pressure  can  be  recovered  along  the  nozzle 
wall  to  increase  net  thrust  and  lift  At  low  Mach  numbers,  once  the  flow  exits  the 
combustor,  die  flow  typically  overexpands  resulting  in  a  fvessure  just  downstream  of 
the  combustor  that  is  less  than  tiie  ambient  pressure  which  results  in  increased  drag 
and  decreased  net  thrust 

Any  net  performance  improvements  to  hypersonic  propulsive  nozzles  could 
have  a  significant  impact  on  hypersonic  flight  vehicle  designs  like  the  NASP  or 
NASP-dmived  vehicles.  An  improved  cowl  c<»cq>t  has  been  {voposed  whidi  may 
have  some  potential  for  increasing  the  overall  performance  of  these  vehicles.  This 
omcept  c<msists  of  pladng  a  rocket  or  gas  generator  into  the  cowl  to  supply  a 
relatively  higb  pressure  boundary  condition  to  the  main  exhaust  plume. 

1.2  Purpose 

There  are  many  potential  benefits  for  placing  a  throttleable,  vectcxable  gas 
generator  into  the  cowl  of  a  hypersonic  vehicle.  The  gas  generator  could  potentially 
pravido  weight  reductitm  due  to  a  decreased  cowl  length,  increased  direct  thrust  from 
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the  gu  generator,  end  a  potential  increase  in  nozzle  peffonnance  due  to  higher 
pressure  on  the  nozzle  wall.  In  addition,  die  (rflwt  locadon  of  the  gas  generator  in  the 
cowl  may  make  it  well  suited  for  providing  direct  stability  and  control  m<Mnent  c(»itrol 
to  the  vdiicle. 

This  computadonal  analysis  investigates  die  effect  of  inserting  a  gas  generator 
flow  into  a  nozzk/cowl  flowfield. 

U  Scope 

The  scope  of  this  research  effort  incliKled  implementing  a  gas  generator  flow 
model  into  an  existing  crmiputadonal  fluid  dynamics  (CTO)  code  and  evaluating  the 
potential  (m-design  and  off-design  perfomance  of  this  concqit 

The  effect  of  a  gas  generator  on  nozzle  performance  was  investigated  for  two 
different  nozzleAxiwl  gerxnetries.  An  experimentally  validated  hypersonic  nozzle/cowl 
configuratirm  was  analyzed  at  relatively  low  qieed,  off-design  conditioos  of  Mach  1.9 
and  3.0.  A  more  realistic  generic  nozzle/bowl  configuration  was  analyzed  over  a 
typical  trajectory  which  maintained  a  constant  dynamic  jxessure,  q.  of  1000  psf 
dirou^  <xi-design  flight  crmdidrMis  ranging  between  Mach  10  and  25. 

The  first  nozzle  was  based  on  an  expoimental  hypersonic  nozzle/aAerbody 
modd  developed  by  Cochran  (4).  This  is  also  the  same  geometry  used  by  Hyun  (1)  to 
validate  the  computational  code  used  for  tiiis  research.  Figure  1-2  shows  the 
e]q)eiimental  morfel  and  figure  1-3  shows  the  experimental  apparatus  that  this 
configuration  was  based  on. 
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The  second  nozzle/oowl  geometry  aoelyzed  represents  a  more  realistic  gmieric 
hypersonic  vehicle  afterbody  configuration  which  was  analyzed  over  a  typical  flight 
trajectory.  Figure  1-3  shows  die  generic  nozzle/cowl  geometry  (feveloped  by  Doty  (9) 
that  was  used  for  this  portton  of  die  stikly.  T!»  generic  nozzle/cowl  geometry  was 
optimized  by  Boniparte  (2)  over  this  typical  trajectory.  The  nozzle  cmifiguratimi 
chosen  for  this  analysis  corresponds  to  Brxu^iarte's  optimized  generic  nozzle  geometry 
evaluated  at  a  design  Mach  number  of  15. 

Realistic  gas  generator  exhaust  cmxlitions  were  obtained  using  a  code 
developed  to  detormine  the  equilibrium  flow  solution  for  a  liquid  hydrpgen-oxygen 
rocket  motor.  The  gas  generator  con^ureti<»  was  based  on  a  scaled  down  Space 
Shuttle  Main  Engine  (SSMB)  gemnetry  (22). 

The  gas  generator  flowfield  was  apfdied  to  bodi  of  the  nozzle/cowl 
oonfiguFBtions  and  the  flov  solution  was  evaluated  at  a  variety  of  gas  graerator 
operating  conditicMis.  Not  mily  were  the  effects  of  gas  generator  mass  flow  ai^  exit 
pressure  aiudyzed,  but  the  effects  of  gas  generator  thrust  vector  deflection  angle  were 
also  investigated. 

For  the  eaqpeiimental  nozzle/cowl  configuration,  the  gas  generator  effects  were 
evaluated  over  a  range  of  NFR  at  two  low  ^leed  off-design  Mach  numbers,  Mach  1.9 
and  Mach  3.0.  Tables  1-1  through  1-3  describe  the  initial  cmiditimis  analyzed 
computatimially.  For  each  Mach  number,  the  gas  gmierator  mass  flow  was  varied 
frmn  die  baseline  gas  generator  off  ctmditioo  (0%  throtde  setting)  to  a  200%  of  design 
opo'ating  conditimi  throtde  setting  in  25%  increments.  The  gas  generator  deflection 
effects  were  analyzed  at  a  100%  throtde  setting  over  a  range  of  ±10  degrees  in  2.5 
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(k^ree  increments.  A  negative  gas  generator  deflection  angle  is  defined  as  a 
(teflectimi  towards  the  nozzle  wall  and  a  positive  deflection  is  defined  as  a  deflectimi 
away  frmn  the  nozzle  wall.  The  baseline  gas  generator  <0  cases  were  compared  to 
die  flow  soludmi  using  a  non-deflected  100%  throttle  gas  generator  flow  over  a  range 
of  operating  NPR.  A  total  of  4S  final  flow  soluticMis  were  calculated  for  this 
eiqierifflental  nozzle/cowl  getunetry. 

The  gas  generator  effects  are  also  compared  to  the  effects  of  varying  the  cowl 
wall  geometry.  Cowl  wall  trailing  edge  extensions  (+0.5  inch)  and  deflections  (±5 
degrees)  were  analyzed  experimentally  by  Cochran  and  numerically  by  Hyun.  Gas 
generator  effects  are  compared  direcdy  to  these  results. 

For  the  generic  nozzle/cowl  configuradtm,  the  results  were  analyzed  over  a 
typical  trajectory.  Tables  1-4  and  1-S  describe  the  cases  evaluated  for  this  analysis. 
The  gas  generator  effects  were  analyzed  over  4  points  spanning  a  typical  trajectory. 
Mach  10,  IS,  20  and  25  flight  conditions  were  evaluated  over  a  constant  1000  psf 
dynamic  pressure  trajectory.  Gas  generator  dirottle  setting  was  increased  from  0%  to 
200%  of  the  design  operating  condition  in  50%  increments.  The  gas  gmierator 
deflecdtm  angles  covered  ±10  degrees  in  5  (tegree  increments.  Excursions  were  also 
made  to  analyze  the  trends  associated  with  relatively  large  gas  generator  deflection 
angles  of  ±20  degrees.  A  total  of  44  Enal  flow  solutitms  were  calculated  for  this 
generic  nozzWcowl  geometry. 

Tim  effects  of  the  gas  generator  are  presented  in  terms  of  overall  nozzle 
performance.  Since  a  significant  percental  of  total  thrust  is  generated  through 
pressure  recovery  along  the  airframe-integrated  nozzle  wall,  the  effects  of  the  gas 


generator  are  presented  in  terms  of  the  pressure  distributiim  along  the  nozzle  wall. 

The  net  effect  of  these  pressure  recovery  changes  are  presented  in  terms  of  the 
integrated  nozzle  wall  pressure  contribudons  to  lift  and  thrust.  For  the  trajectory 
analysis  of  the  generic  nozzle/cowl  configuradm,  the  fuessure  contribudons  to  lift  and 
thrust  are  also  presented  reladve  to  the  overall  vehicle  afterbody  performance. 

1.4  Approach 

Hypersonic  vehicle  afterbody  flowfields  contain  a  complex  interacdon  of 
flowfield  phenomena  including  e]q)ansion  waves,  shock  waves,  viscous  crmtact 
surfaces,  and  solid  surface  boundaries  including  the  cowl  and  nozzle  walls.  The 
flowfield  is  complicated  with  the  strong  interacdon  of  all  of  these  components. 

The  literature  shows  that  many  computational  fluid  dynamics  (CPD)  codes 
have  been  developed  to  analyze  the  complex  flowfield  of  hypersonic  propulsive 
nozzles.  Any  CPD  code  used  to  arudyze  this  type  of  complicated  flow  must  be 
validated  against  experimental  data  to  ensure  that  it  can  accurately  capture  the  effects 
of  all  of  the  complicated  flowfield  interactimis.  One  such  code,  developed  by  Wright 
Laboratory,  solves  the  unsteady,  two-dimensional  Navier-Stokes  equations  based  on 
perfect  gas  laminar  flow  assumptions  using  several  mq>lidt  flux-splitting  finite  volume 
schemes.  Hyun  (1)  cmnpared  the  Wright  LaboraUxy  code  to  the  experimental  data  of 
Cochran  (4).  Hyxm  demonstrated  that  this  code  captured  the  complicated  flowfield 
phenomena  quite  well  for  high  speed  freestream  flow  conditions.  Hyun  also  showed 
that  as  NPR  and  freestream  Mach  number  increased,  the  geno-al  crurelation  between 
experimental  and  computational  results  improved.  Fm*  these  reasons,  this  code  was 


chosen  to  analyze  the  effects  of  placing  an  atkiidonal  gas  generator  into  this 
complicated  flowffeld. 

The  flux-vector  flitting  (FVS)  approach  of  Van  Leer  and  the  flux-difference 
splittiqg  (FDS)  approach  of  Roe  were  incorporated  into  the  Wright  Laboratory  code  by 
GaiUmde  (7).  These  robust  algorithms  are  second-order  accurate  in  space  and  time. 

The  Van  Leer  method  was  initially  used  to  achieve  a  steady  state  solution  since 
it  requires  no  more  than  two  internal  cells  to  ci4>ture  internal  shock  structures.  Hyun 
(1)  discovered  that  the  Roe  method  should  not  used  for  the  initial  steady  state  solution 
because  it  may  violate  the  entropy  condition  near  sonic  or  stagnation  points.  When 
the  Van  Leer  soluti<»  converged  to  a  steady  state  solution,  the  flowfield  calculations 
were  then  completed  using  the  Roe  ipproach,  which  more  accurately  captures  wave 
interactions  between  cell  interfaces  (1). 

Several  grid  densities  were  evaluated  to  perform  a  grid  resolution  study  on  die 
generic  nozzle/cowl  configuration  investigated.  Once  a  suitable  grid  was  chosen, 
complete  flow  solutirms  were  calculated  for  each  desired  flow  condition. 


7 


Table  1*1  Eqjenmental  Nozzle.  Mass  Flow  Cases 


Mach 

Gas  Generator  Throttle  Setting  [%]  j 

1.9 

a 

25 

50 

75 

100 

125 

150 

175 

200 

3.0 

0 

25 

50 

75 

100 

125 

150 

175 

200 

Table  1*2  Experimental  Ncozle,  Thrust  Vector  Deflection  Cases 


Mach 


Gas  Gentfator  Thrust  Vector  Deflectitm  Angle  [deg] 


1 

1  -‘OO 

-7.5 

-5.0 

-2.5 

0.0 

2.5 

HI 

m 

10.0 

1 

1 

-7.5 

-5.0 

-25 

0.0 

2.5 

d 

la 

10.0 

Table  1-3  Experimental  Nozzle,  Nozzle  Pressure  Ratio  Cases 


1  Mach 

Nozzle  Pressure  Ratio  (NPR)  | 

1 

3.0 

5.0 

7.0 

12.0 

N/A 

1  3.0 

j 

7.0 

9.0 

12.0 

16.0 
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Table  1-4  Generic  Nozzle,  Mass  Flow  Cases 


Mach 

Gas  Generator  Throttle  Setting  [%]  | 

10 

0 

50 

100 

150 

200 

15 

0 

50 

100 

150 

200 

20 

0 

50 

100 

150 

200 

25 

0 

50 

100 

150 

200 

Table  1-5  Generic  Nozzle,  Thrust  VecUx*  Deflection  Cases 


Mach 

Gas  Generate'  Thrust  Vect(M‘  Deflection  Angle  [deg]  | 

10 

-20 

-10 

-5 

0 

5 

10 

20 

15 

-20 

-10 

-5 

miQiiiii 

5 

10 

20 

20 

-20 

-10 

-5 

0 

5 

10 

20 

25 

-20 

-10 

-5 

0 

5 

10 

20 
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Figure  1*2  Experimental  Nozzle/Cowl/Afterbody  Model  (4) 
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OrOBB  BBCiiOIl 
of  the  Btrut 


Nozzle  wall 


Figure  1-4  Expanded  View  of  a  Generic  Hypersonic  Nozzle/Cowl  Configuration  (9) 
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n.  GOVERNING  EQUATIONS  AND  NUMERICAL  METHODS 
2.1  Governing  Equations 

The  governing  equations  of  motion  for  die  planar,  unsteady,  viscous  flow  of  a 
compressible  fluid  which  account  for  viscous  stresses  and  thermal  conduction  are  the 
Navier-Stokes  equations.  By  neglecting  body  fmces  and  heat  sources,  the  Navier- 
Stokes  equations  may  be  written  (10): 


(Continuity) 

ar  axj 


£^_0  (Momentum) 

at  dXj  axt  dXj 


at  dXj  dXj 


(Energy) 


(Equation  of  state) 

P  2 


(2-1) 


(2-2) 


(^3) 


(2-4) 


where  indicial  notation  has  been  used  (rq>eated  indices  are  summed)  and 

component  of  velocity  in  Xf  direction 
viscous  stress 

V 
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dUi  dU.  dU, 


(2-5) 


where  6^-  the  Kr<necker  delta 

ln^X*  and  second  coefficients  oi  viscosity,  respectively 

qj^-kidDaxp 

where  ik  is  the  coefficient  of  thomal  conductivity. 

The  coefficient  of  viscosity  is  related  to  the  thermodynamic  variables  using 
Sutherlands  Law  (10): 


p=Ci 


3 

Ti 

r+c 


(2-6) 


where  Cj  and  are  constants  for  a  given  gas.  For  air  at  standard  temperatures,  Ci  ~ 
2.27B08  lV(ft-8ec  and  198.9  R. 


The  two-dimensional  Navier-Stokes  equations  i»esented  in  Equatimis  (2-1) 
through  (2-4)  can  be  rewritten  in  vector  divergence  form.  The  vectw  divergence  form 
(tf  diese  equations  simplifies  their  implemoitation  into  finite  volume  computational 
fluid  ^namics  algorithms.  In  vector  divogence  form,  these  equations  may  be  written: 


i!?[+.^+i^=0  (2-7) 

dr  ar  dy 


14 


li^iere  die  consorved  variables  are  given  by: 


(2-8) 


and  the  flux  vectors  are  separated  into  tl»ir  inviscid  (F.G)  and  viscous  terms  iF,.GX 


FA 


pll 

pv 

pu*+p 

G= 

puv 

ptfV 

pv^+p 

(£,+p)a 

(£,+p)v 

(2-9) 


0 

0 

“^V*f 

(2-10) 


The  final  equatitm  required  to  define  the  flow  field  is  die  perfect  gas  equation-of-state. 


p=iy-V)pe 


-|(«*+v*) 


(2-11) 


2J,  Numerical  Methods 

2.2.1  Coordinate  Transformation.  The  Navier-Stokes  equati(»is  are 
transformed  fi'om  die  physical  space  cotmlinate  system  to  a  ctunputational  space 
coordinate  system  in  (xder  to  generate  a  suitable  computational  grid  fcx  the  qiplication 
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of  the  numerical  method.  The  coordinate  transformation  simplifies  the  process  of 
Inlying  boundary  conditions,  clustering  grid  points,  ami  providing  orthogmudity  for 
application  of  numerical  solutions  to  tile  Navier-Stokes  equations  (10).  Equation  (2-7) 
can  be  transformed  into  a  general  coordinate  system 

(2-12) 

dr  dn 

The  details  of  this  transfarmatiai  can  be  found  in  (6)  and  (12).  The  general 
transformed  form  of  the  Navier-Stokes  equations  can  be  written  (1): 

where  the  Jacobian  of  the  transformatitti  is  given  by: 

2.12  Discretaatkm.  This  study  uses  a  modified  version  of  a  code  that  was 
developed  by  Wrigbt  Laboratory  (WL/FIMQ.  The  code  has  three  options  available 
for  the  discfetization  of  the  Navier-Stokes  equations;  Steger  and  Waimiog  flux-vector 
splitting.  Van  Leer  flux-vector  flitting,  and  Roe  flux-difference  flitting.  The  viscous 
flux  vector  terms  are  discretized  using  simple  central  differences.  The  inviscid 
flux  vector  terms  (F,G)  may  be  discretized  using  the  either  of  the  flux-vector  splitting 
or  flux-difference  splitting  methods.  A  detailed  discussion  and  examples  of  these 
three  discretizations  can  be  found  in  Hyun  (1). 


(2-13) 


(2-14) 
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Tbe  numerical  tf^nxMch  ia  aecood-order  accurate  in  time  and  ^mce.  The 
scdutioQ  ia  marched  in  time  uaing  a  two-atage  Runge-Kutta  acfaeme.  The  formulation 
ia  finite  volume  and  aeccmd-ocder  accuracy  ia  obtained  diough  the  uae  of  the  MUSCL 
(MonoUmic  Upstream  Sdiemes  for  Convosirm  Laws)  approach  in  conjunction  with  a 
minmod  limiter  to  prevent  oscillationa.  Details  of  tbe  MUSCL  apf^oach  and  the 
minmod  limiter  i^lied  to  a  one-dimensitmal  example  can  also  be  found  in  Hyun  (1). 

The  basic  concq>t  of  tbe  flux-vector  and  flux  difference  splitting  can  be 
illustrated  using  die  fdlowing  <me-dimensional,  inviscid  model: 

(MS) 

at  ax 

A  first-order  eaqilidt  discietization  of  Equation  (2-lS)  uring  the  forward  Euler  mediod 
gives: 


(at/*), 

At  AX 


(M6) 


where 


(2-17) 


and  where  i*  based  on  some  combination  of  F  at  the  adjacent  grid  pomts.  For 
example,  a  simple  average  value  could  be  used  where  F^,y2KF/fF^,)/2. 

The  two-dunensional  Navier-Stokes  Equation  (2-7)  can  be  discretized  in  a 
similar  maimer.  A  successive  balandng  of  the  fluxes  in  each  direction  is  used  to 
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compute  the  rig^t  hand  tide,  or  residual,  of  die  twoKlimensioiial  versiofi  of  Equation 
(2-16): 


^ - 

'  AX 

(2-l«) 

LHS,^LHS,* — 2 - 2 - 

where  £J75,  is  the  left  hand  side  of  the  two-dimensional  version  of  Equation  (2-16). 
The  left  hand  side  of  Equation  (2-16)  is  calculated  using  a  two  stage  Runge-Kutta 
scheme  known  as  Heun’s  method.  The  details  of  this  Runge-Kutta  ai^licatimi  can  be 
found  in  (24). 


The  ri^t  hand  side  of  Equation  (2-16)  can  be  evaluated  using  either  the  flux- 
vector  ^littiqg  approach  of  Stager  and  Warmiqg  or  Van  Leer  or  die  flux-diflerence 
splitting  approach  of  Roe.  The  Van  Leer  apfx'oach  is  essentially  an  improvement  over 
the  earlia'  St^er  and  Warming  flux-vector  splitting  approach.  Inaccuracies  arise  in 
die  St^er  and  Warming  scheme  whmi  the  eigenvalues  of  die  flux  Jacobians  change 
sign.  The  flux  Jacobians,  A  and  B  are  gi\^  by: 

F^AU ;  G=BU  ;  (2-19) 

dU  dU 

eigenvalue  sign  changes  will  occur  near  shock  structures  or  st^gnatimi  points 
in  the  flowfidd.  In  the  Steger  and  Warming  scheme,  die  split-flux  (terivatives  are 
discmitinuous  at  these  dgenvalue  sign  changes  giving  rise  to  inaccuracies  and 
instabilities  in  these  regicHis.  Van  Leer’s  improvement  over  the  Steger  and  Warming 
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scheme  was  to  split  die  fluxes  so  that  the  forward  and  backward  flux  contributicMis 
transitioned  smoothly  at  the  eigenvalue  sign  changes.  A  cmnpariscm  of  these  two 
methods  is  (xesented  in  Hyun  (1). 

The  Van  Leo*  q)iMroach  was  chosen  as  the  initial  flow  solver  for  this  study 
since  it  is  more  robust  than  the  aiul  Warming  approach.  Van  Leer’s  approach 
to  solving  the  right  hand  side  df  Equation  (2-16)  consists  splitting  the  flux  into 
positive  and  negative  components  for  approixiate  upwind  differencing; 

Fi,i=F*(Uh*F~(V‘^  (2-20) 

where  £/*  and  iF  are  vectors  of  die  omserved  flow  variables  at  each  interface.  The 
Van  Leer  flux-vector  flitting  i^oadi  is  designed  to  make  computational  schemes 
more  robust  and  to  imixove  computational  efficiency.  The  goal  of  a  flux-vector 
splittiQg  scheme  is  to  split  the  flux  vectors  in  such  a  way  that  an  iqiwind  finite- 
difference  scheme  may  be  used  throughout  the  entire  flowfield.  This  is  accomplished 
by  flitting  die  flux  vectors  into  upwind  and  downwind  parts.  For  second-order 
accuracy,  V*  and  f/"  are  obtained  using  Van  Leer’s  MUSCL  approach  to  extrapolate 
the  consoved  variables  to  the  cell  surface.  Van  Leer’s  qifxoach  also  incoqiorates  a 
minmod  lifflit.r  based  on  local  coodidtMU  to  avoid  osdllations  in  shock  regions.  The 
addititm  of  die  mnmod  limiter  causes  the  solution  to  locally  revert  back  to  first-oder 
accuracy  in  the  vicinity  of  the  shock  in  order  to  preserve  monotonicity  within  the 
solution. 
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The  Van  Leer  •|>proech  is  also  very  robu&t  in  regirns  of  supersonic  expansion. 
This  mediod  is  useful  for  starting  the  flow  solution  when  there  is  a  significant  amount 
of  e]q)ansi<m  from  the  initial  crniditioos  in  the  steady  state  flow  solution  (1). 

Roe’s  scheme  is  used  to  achieve  a  final  solution  after  the  Van  Leer  solution 
has  achieved  a  steady  state  solution  throu^  the  initial  flowfield  expansions.  Roe’s 
^^oach  to  solving  the  right  hand  side  of  Equadm  (2*16)  is: 


1 

2 


F(U‘^*FiU*)--A(U^-Uh 

2 


(2-21) 


where  2  is  the  Roe-averaged  flux  vector  Jacobian  matrix.  The  refers  to  Roe 
averaging  given  by: 


x‘f.u.y,and  h  (2-22) 

where  die  R  and  the  L  subscripts  refer  to  the  components  of  I/*  and  respectively. 
The  Roe  ^ifxoach  also  uses  a  MUSCL  q^roach  and  a  minmod  limiter  to  achieve 
second-order  accuracy. 


2.3  Boundary  and  Initial  Conditions 

Since  the  finite-volume  computational  ^diemes  used  for  this  study  are  cell 
centered  methods  which  solve  for  fluxes  across  the  cell  surfaces,  ghost  pdnts  must  be 
added  to  the  grid  in  mder  to  apply  boundary  conditions.  Ghost  points  are  required  at 
all  internal  and  external  boimdaries.  Hgure  2-1  demonstrates  the  triplication  of  ghost 
points  to  a  simplified  grid  and  it  helps  highlight  the  necessity  of  ghost  points.  For 
example,  if  the  surface  boundary  depicted  in  Hgure  2-1  requires  a  no-slip  boimdary 
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coiiditi<m  there  is  no  way  to  specify  zero  velocity  <»  this  surface  with  Um  basic  cell 
centered  grid.  Ghost  grid  point  nodes  are  i^^lied  such  that  a  row  of  new  tictitious 
cells  are  created  as  extrq)olati(»u  of  die  original  grid.  The  extrapolated  grid  has  an 
associated  row  of  fictitious  cell  centers  that  are  outside  of  the  original  surface 
boundary  as  shown  in  Figure  2-1.  Flow  variables  can  then  be  specified  at  the  ghost 
cell  centers  such  that  the  desired  boundary  conditions  are  enforced  on  the  surface 


boundary.  For  this  sin^iUfied  example,  where  the  internal  and  external  cell  centers  are 
the  same  distance  from  the  ed^  of  the  grid,  a  no-slip  boundary  conditimi  can  be 


Figure  2-1  Addition  of  Ghost  Points  to  a  Grid 
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applied  on  the  edge  of  dw  grid  by  simply  setting  the  velocity  components  in  the  new 
fictitious  cells  equal  to  and  opposite  firmn  the  velocity  compcmoits  specified  in  their 
nmg^boring  internal  edge  ceUs.  This  ^)proach  creates  a  velocity  gradient  from  the 
last  internal  cell  to  die  fictitious  ghost  ceU  such  that  the  velocity  is  zero  on  the  surface 
boundary. 

Boundary  conditions  must  be  applied  at  all  of  the  edges  and  surface  boundaries 
associated  with  d^  physical  problem  being  modelled.  The  nozzlc/cowl  flowfield 
requires  inflow,  freestream,  outflow,  and  surface  boundaries  conditions  to  be  applied. 

The  inflow  boundary  conditions  must  be  specified  for  several  difiimnt 
flowfields  fw  the  hypersmiic  vehicle  aft  body  problem.  The  internal  flowfield 
correqxmding  to  the  main  combustion  chamber  exhaust,  the  external  flowfield  that  the 
vehicle  aft  body  is  exposed  to.  and  the  new  gas  generator  flowfield  added  to  die  code 
for  this  study  all  require  inflow  boundary  conditions  to  be  specified.  At  these  inflow 
boundaries,  the  flowfield  is  specified  according  to  given  initial  conditions.  The 
required  flow  variables,  V,  are  calculated  from  the  inflow  initial  cmiditions.  On  the 
inflow  faces,  conesponding  to  lines  AO  and  HI  in  Figure  1-4  for  the  internal  and 
external  flowfields  respectively,  the  boundary  conditions  are  q^lied  by  setting  the 
flow  vector  in  the  boundary  cells  equal  to  die  input  flow  vector. 

(2-23) 

For  the  freestream  boimdary  cmresponding  to  line  U  in  figure  1-4,  the  initial 
omditimis  are  specified  as  the  freestream  values: 
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(2-24) 


The  rero  gradient  boundary  condition  is  then  applied  by  simply  setting  the  flow 
vector  in  the  boundary  cell  equal  to  the  flow  vector  in  the  last  cell  of  the  flow 
solution: 

<*•») 

The  outflow  boundary  condition  corresponding  to  line  CJ  in  Hgure  1-4  also 
has  a  similar  no-diange  boundary  condititm  applied: 

(2-20 

Tiie  solid  boundaries  of  the  nozzle  wall  cmresponding  to  curve  ABC  in  Hgure 
1-4  and  all  of  the  cowl  lip  surfaces  conresponding  to  line  OEFH  in  Figure  1-4  must 
have  the  no-slip  boundary  condition  applied.  For  example,  on  the  nozzle  wall: 


(^27) 


The  pressure  gradient  also  goes  to  zero  on  the  wall: 


i&=0  (2-28) 

For  a  zeroth-order  implementation,  this  gives: 

The  Wright  Laborafoxy  code  used  for  this  analysis  also  requires  a  definition  of 
the  wall  temperature.  An  adiabatic  wall  temperature  assumption  was  briefly 
investigated  for  the  generic  nozzle  configuration  analyzed  in  this  research.  An 
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approximati<m  for  adiabatic  wall  temperaOire  for  compressible  boundary  layers  is 
given  by  (13): 

where  r.  is  the  adiabatic  recovery  fachx-  which*  for  laminar  compressible  boundary 
layers,  is  approximately  equal  to  the  square  root  of  the  Prandtl  number,  Pr.  The 
Prandtl  number  for  air  up  to  moderately  high  temperature  is  Pr  :=  0.72. 

Fw  air  at  the  c(»ditions  specified  fw  Mach  IS  flow  in  Table  3-4,  the  adiabatic 
waU  temperature  would  be  approximately  16610  "R  ftM*  the  internal  nozzle  wall  at  the 
inlet  and  q)proximately  18280  ^  on  the  external  cowl  wall.  These  temperatures  are 
much  higher  than  the  melting  temperature  c/[  most  advanced  materials.  Therefore,  the 
vehicle  requires  active  cooling  on  the  upstream  nozzle  wall  and  the  cowl  wall.  Active 
cooling  was  assumed  such  that  the  wall  temperature  could  be  maintained  below  a 
reasonable  target  design  temperature.  A  typical  upper  limit  on  nozzle  wall 
temperature  for  a  NASP-^pe  vehicle  is  on  the  o-der  of  2000  ”R  (21).  Although  this 
temperature  would  vary  due  to  an  inoease  in  convective  heat  transfer  as  the  flow 
e]q)anded  down  the  nozzle  wall,  a  constant  2000  °R  was  applied  as  a  reasonable 
q>proximation  over  the  entire  nozzle.  The  errors  associated  with  the  unifmm 
temperature  non-adiabatic  wall  assumption  manifest  themselves  primarily  in  the 
calculation  of  viscous  wall  forces.  Since  the  primary  objective  of  this  study  was  to 
investigate  the  pressure  contribution  to  the  vehicle  performance,  the  errrxs  associated 
with  the  viscous  drag  and  lift  terms  were  not  considered  to  be  critical. 
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The  noii-adiabatic  wall  temperature  is  specifled  as  7^  The  q)ecific  intenial 
energy  can  then  be  specified  for  the  boundary  cell  such  that  the  wall  is  maintained  at 


the  given  temperature.  Again,  a  zeroth-odm*  implementation  gives: 

(2-31 

Not  only  does  the  flow  vector  have  to  be  specified  on  the  boundary  cell,  but 
the  fluxes  on  the  edge  cell  must  also  be  specified  for  the  finite  volume  applications 
used  in  this  study.  The  no-slip  condition  on  the  wall  forces  the  inviscid  mass  and 
energy  fluxes  on  the  wall  surface  to  go  to  o.  An  examination  of  Equations  (2-9) 
and  (2-10)  show  that  the  only  terms  remaining  in  the  inviscid  flux  vectors  (F,  G)  are 
the  pressure  terms.  The  invisdd  fluxes  on  the  wall  are  then  specified  by: 


u, 


0 

0 


(^32) 


The  viscous  fluxes  on  the  wall  (F^  G,)  are  simply  calculated  based  on  the 
givmi  wall  temperatures  and  the  zero  velocity  no-slip  condition.  The  same  approach 
described  above  for  the  nozzle  wall  is  also  tq>plied  to  tihe  solid  boundaries  of  the  cowl. 

The  initial  conditions  required  to  start  the  solution  proceuure  are  applied  sudt 
diat  the  entire  grid  is  initially  at  the  input  freestream  conditions  excq>t  at  the  inflow 
boundaries  where  the  iq>propriate  inflow  boundary  conditions  are  applied.  The  nozzle 
internal  and  external  inflow  boundary  conditions  and  freestream  conditions  fix  the 
nozzle  configuration  analyzed  in  Section  3.1  were  chosen  by  Hyun  (1)  to  match  the 
experimental  results  of  Cochran  (4).  A  similar  set  of  initial  conditions  for  the  nozzle 


25 


MXifiguration  analyzed  in  Section  3.2  were  developed  by  Bonaparte  (2)  for  a  generic 
flight  vehicle  over  a  typical  trajectory.  Details  of  the  initial  conditions  used  in  these 
studies  can  be  found  in  Section  3.1  and  3.2,  ie^)ectively. 

The  inflow  boundary  conditions  are  specified  by  the  Reynolds  number,  Mach 
number  and  temperature  of  the  flow.  For  a  calorically  perfect  gas,  all  of  the  required 
flow  variables  and  fluxes  can  be  calculate  given  the  gas  constant  R,  and  specific  heat 
ratio  y.  For  air  at  standard  conditions,  R  =  1716  fP/(s^  °R)  and  ys1.4. 

One  disadvantage  of  the  Wright  Labwatory  code  used  fcx  this  research  is  that 
it  currently  has  no  provisions  for  ftowfield  mixing.  A  hypersonic  propulsive  nozzle 
will  typically  have  different  internal  and  external  flow  gas  pr(^)erties  which  will  mix 
al<mg  a  contact  surface  as  illustrated  in  Figure  1-4.  The  current  computatitmal  code  is 
designed  to  run  only  a  single  calrxically  perfect  gas  mixture  throughout  die  entire 
flowfield.  This  was  considered  an  accqitable  source  of  error  since  this  study  is 
parametric  in  nature  and  the  primary  objective  is  to  determine  the  trends  associated 
with  die  incorporation  of  a  gas  generadx’  into  the  cowl.  Slight  modifications  were 
made  to  the  code,  however,  to  improve  the  consistency  of  the  results.  Separate  values 
for  R  and  y  can  now  be  input  for  each  of  the  three  inflow  boundary  conditions 
(internal  combustix  flow,  external  fireestream  flow,  and  die  gas  generator  flow) 
analyzed  in  this  study.  This  modification  uses  the  input  Reynolds  number,  Mach 
number,  tonperature,  R  and  y  to  calculate  the  correct  pressures  and  fluxes  at  the 
inflow  boundaries.  These  new  inflow  pressures  and  fluxes  are  then  used  in  the  code 
flffijuming  a  single  calraically  perfect  gas. 
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The  following  procedure  is  used  to  convert  the  input  vsrisbles  to  the 
sppro|»iste  pressure  snd  flux  terms.  First,  the  speed  of  sound  and  velocity  are 
calculated  using  die  iq>ut  variables. 


(^33) 

U=aM 

(2-34) 

The  components  of  the  input  velocity  are  then  calculated  using  the  input  flow 
deflectimi  angle,  0. 


u=U‘OosB  (2-35) 

v=f/*8in0  (2-3^ 


Next,  coeffident  of  viscosity,  p,  is  calculated  using  Sudierlands  fmmula. 
Equation  (2-6).  Density  is  then  calculated  using  the  deEnition  of  Reynolds  number. 


Next,  specific  internal  energy  can  be  calculated. 


(2-37) 


(1-38) 

m 

The  final  term  diat  is  required  befme  the  inviscid  flux  vectcHS  (F,G)  given  in 
Equadmi  (2-9)  can  be  determined  is  the  static  pressure.  The  perfect  gas  equation  of 
state  is  used  to  calculate  pressure. 
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p*‘pRr 


(2-39) 


These  updated  flow  variables  and  flux  vectors  are  then  used  as  the  inflow 
boundary  crmditions  for  the  rest  of  the  solution  which  assumes  a  single  calrxically 
perfect  gas  throughout  the  nozzle. 

2.4  Computer  Code  Description 

The  baseline  code  used  to  solve  the  nozzle/cowl  flowfields  analyzed  in  this 
study  was  developed  by  Wright  Laboratory  (WIVFIMC).  This  second-orda*  accurate 
code  was  developed  to  solve  the  two-dimensional,  viscous  Navier-Stokes  equations 
based  oa  a  laminar  flow  assunq>tion. 

The  original  code  was  developed  to  study  uniform  flow  through  a  channel. 
Wright  Laboratory  modified  the  code  in  order  to  solve  the  hypersonic  vehicle  aft  body 
problem  analyzed  by  Hyun  (1).  The  primary  code  modifications  required  for  Hyun’s 
study  included  adding  a  second  flow  stream  to  represent  the  external  flow,  adding  a 
freestream  outer  boundary  condition,  and  inserting  a  cowl  solid  boundary.  One  goal 
of  the  present  study  was  to  further  modify  the  code  to  incorprx'ate  a  third  flow  stream 
in  order  to  model  die  effects  of  a  gas  generator  inserted  into  the  cowl  lip.  A 
description  of  the  required  gas  generate  code  modificatitms  is  presented  in  Section 
2.6.1. 

A  validation  of  the  flux-splitting  algorithms  was  completed  by  Hyun  (1). 
Hyun's  study  compared  the  computational  results  of  this  code  to  the  experimental  data 
of  Codhran  (4)  using  the  same  baseline  geometry  and  grid  presented  in  Section  3.1  of 
this  study.  Hyun  showed  that  as  Mach  number  and  nozzle  pressure  ratio  (NPR) 
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increised,  the  conelttioa  between  the  computational  results  and  the  experimental  data 
generally  improved.  The  cmnputations  based  on  the  assumption  of  planar  laminar 
flow  tended  to  be  more  consistent  widi  the  experimental  data  as  the  Mach  and  NPR 
q){N‘oached  their  on-design  conditions.  The  promising  results  for  the  supersonic  flow 
conditimis  and  the  trend  towards  increased  cr^elatitm  at  higher  Mach  and  NI% 
provided  the  required  justification  for  using  this  code  in  the  present  study. 

The  Wright  Laboratory  code  provides  many  options  for  solving  the  given  flow 
{tfoblem.  The  code  can  be  used  to  solve  the  inviscid  Euler  equations  as  well  as  the 
viscous  Navier-Stokes  equations.  The  code  also  allows  the  user  to  choose  which  type 
of  flux-splitting  algmithm  to  use;  Steger  and  Warming  flux-vector  splitting.  Van  Leer 
flux-vector  splitting,  or  Roe  flux-differeiree  splitting  can  be  selected.  The  code  also 
has  a  restart  feature  which  allows  any  combination  of  diese  solutions  to  be  run  in 
series.  For  example,  the  jHesent  study  uses  the  steady-state  Van  Leer  solution  as  the 
initial  conditions  fcv  the  final  Roe  solution.  A  general  discussion  of  the  solution 
procedure  used  is  presented  in  Section  2.5. 

In  addition  to  choosing  the  type  cS  flow  solver,  the  code  allows  the  user  to 
contrcd  other  aspects  c/[  the  problem  including;  cmitrol  ciS  the  CFL  stability  criteria, 
input  of  flowfield  and  boundary  conditions,  and  cmitrol  of  miscellaneous  fmmat  and 
run  parameters.  A  sample  input  file  and  a  description  of  the  inputs  fw  can  be  found 
in  Ai^>endix  A. 

The  code  is  written  entirely  in  FORTRAN  and  it  was  run  on  a  netwwk  of 
Silicrm  CSraphics  Iris  4D/240  and  Indigo  series  wwkstations. 
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13  Solution  Procedure 


The  e]q>licit  formulatim  of  the  numerical  problem  was  used  for  this 
computaticmal  analysis.  This  i^^oadi  is  second-order  accurate  in  space  and  time. 

The  ejqilicit  formulation  was  desired  because  it  is  well  suited  for  optimization  on 
parallel  processing  machines.  Section  2.6.2  and  Appendix  B  highlight  the  code 
optimization  that  was  completed  for  this  research. 

The  apfH’oadt  taken  to  solve  the  nozzle/cowl  problem  uses  a  combination  of  the 
Van  Leer  flux-vector  q)litting  method  and  the  Roe  flux-difference  splitting  method. 
The  initial  condititms  required  to  start  the  iterative  solution  algorithm  are  initially  set 
to  the  i]:q)ut  freestream  ccmditions.  As  tlm  flow  solutimi  is  marched  in  time,  the  initial 
perturbations  are  very  large.  The  initial  cmiditions  for  the  numerical  solution  are 
significantly  underexpanded  and  as  the  solution  marches  towards  steady-state,  tlm  code 
must  overcome  this  initial  undeiexpansion.  The  code  provides  the  user  with  direct 
contrd  over  the  Courant-Fdedridis  and  Lewy  (CPL)  stability  criterion  defined  as  (1): 

CFL=-^  (2-40) 

where  is  the  largest  eigenvalue  of  the  flux  Jacobian  matrix. 

An  iiutial  CPL  of  0.01  is  used  to  force  very  small  initial  steps  in  time  through 
tile  initial  expansion  of  the  flow.  The  CX^  is  then  gradually  increased  to  a  maTimnm 
of  0.9  as  the  solution  is  marched  towards  steady-state. 

The  Roe  ipp-oadi  has  demonstrated  convergence  problems  ^i^ien  trying  to  start 
extremdy  underexpanded  initial  conditions  for  the  numerical  solutirm  (1).  For  this 
reason,  the  Van  Leer  approach  was  chosen  to  generate  an  initial  steady-state  solution 
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became  it  has  demonstrated  robustness  thrmigh  these  initial  supersonic  e]q)ansion 
regioas.  Oi^  a  steady  solution  was  achieved  using  the  Van  Leer  ^>proach,  a  final 
solution  was  generated  using  the  Roe  f!ux>dif]fereace  q)litting  approach  to  accurately 
capture  the  wave  interactions. 

Convergence  to  the  steady-state  solution  was  determined  e3q)licitly  for  each  run 
in  this  study.  The  explicit  criteria  used  to  determine  the  convergence  are  the  root 
mean  square  (RMS)  beat  transfer  and  jxessure  values  along  the  nozzle  wall; 

Convergence  is  detomined  when  and  do  not  fluctuate  more  than 
0.01%  over  a  large  number  of  iteratiais.  The  ^ical  number  of  iterations  required  for 
a  converged  Van  Leer  solutimi  ranged  frmn  750  to  2000.  An  additimial  500  to  4000 
Roe  iterations  were  then  required  on  top  (rf*  this  to  achieve  a  final  steady-state  solution. 
As  fieestream  Mach  number  and  NPR  were  increased,  die  number  of  iteratirms 
required  for  a  steady  state  solution  tended  to  decrease.  This  trend  also  tends  to 
suf^HXt  the  notitm  diat  the  viscous  laminar  flow  assumption  used  in  this  code  is  well 
suited  for  the  hypersonic  nozzle/cowl  problem. 
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Since  one  of  the  primary  goals  oi  this  research  was  to  evaluate  the  effect  ci 
gas  generator  flowfield  <»  the  pressure  distribution  along  die  nozzle  wall,  several 
rq)reaentative  measures  of  performance  are  i»esented.  The  coefficient  of  pressure  (C^ 
and  die  fnessure  contributimi  to  lift  and  thrust  were  calculated  along  the  nozzle  wall. 

The  pressure  coeffidait  almig  die  nozzle  wall  is  calculated  locally  using: 


(i43) 


where  and  are  the  fireestieam  static  and  dynamic  pressure,  respectively. 
The  pressure  contributitm  to  thrust  along  the  nozzle  wall  is  given  by: 


^  <-i 

vdiere  is  calculated  from  die  wall  grid  spacing  and  p,  is  calculated  at  the  cell 
center. 

The  fxessure  c<mtribution  to  lift  along  the  nozzle  wall  is  calculated  in  a  similar 
manner. 


1-1 

whme  Ax^  is  also  calculated  based  on  tire  wall  grid  spacing. 


(2-45) 


2.6  Computer  Code  Modifications 

In  addition  to  the  small  duuge  in  the  inflow  boundary  conditions  desoibed  in 
Section  2.3,  several  more  significant  changes  were  made  to  the  Wright  Labmatiwy 
code  for  the  incorporation  of  a  vectwable  gas  generahx’  flowfield  into  the  cowl  and 
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die  ofitimizatioa  of  the  code  for  use  (»  parallel  prooeasing  machines.  The  following 
secticms  highli^t  these  code  diangea. 

2.6.1  Addition  of  a  Gas  Generator  Flowfield.  The  data  ii^iut  subroutine, 
DATING,  and  the  commcm  block  routine,  gcommon.f,  were  modified  to  read  in  the 
required  gas  generator  flow  <tefinition  ftx>m  the  input  file  cnldat  found  in  Appendix  A. 
The  new  variables  are  Reynolds  number,  Mach  number,  tmnperature,  R.  y,  and  thrust 
vector  deflection  angle  6.  The  variables  are  then  converted  to  the  desired  flow 
variables  as  described  in  Section  2.3. 

The  initial  condition  subroutine,  INITL.f,  was  modified  to  set  the  flow  variables 
and  die  constfved  variables,  U,  to  the  proptf  initial  conditions  on  the  gas  generator 
boundary. 

The  boundary  condition  subroutine,  BC.f,  required  dm  most  extensive 
modifuxUimis.  The  solid  surface  boundary  condition  at  die  cowl  lip  was  rq)laced  with 
an  inflow  boundary  condition  as  described  in  Section  2.3. 

The  main  calculation  subroutine,  L/,  also  required  a  slight  modification.  In  the 
baseline  code  (no  gas  goierator)  the  cowl  Iq)  was  treated  as  a  solid  boundary.  This 
required  that  the  invisdd  flux  terms  cm  the  surface  be  q>ecified  as  shown  in  Equatirm 
(2-31)  where  the  mass  and  energy  terms  are  set  to  zero  and  only  the  pressure 
contribution  to  the  momentum  terms  remain.  In  the  modified  code,  the  cowl  lip  is  no 
lon^  fflockilled  as  a  solid  boimdary  and  the  flux  vectors  are  not  reset.  These  terms 
are  now  set  equal  to  the  flux  tmms  associated  with  the  gas  geno’ato’  flowfield  initial 
conditirms. 
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2.<U  Code  Optimization.  The  steady-state  solution  of  a  t^ical  nozzlc/cowl 
coafigurati<«  required  oa  die  order  of  4  to  20  hours  of  computer  run  time  using  the 
baseline  code.  A  significant  reducti<»  in  run  time  was  required  in  order  to  complete 
the  large  number  of  runs  required  for  this  research  in  a  timely  manner.  The  run  time 
could  be  reduced  because  the  eiqilidt  formulatioo  of  the  numerical  problem  is  well 
suited  foe  optimization  using  parallel  processing  machines  such  as  the  Silicon  Graphics 
ADPAO  used  for  most  of  this  research. 

The  Silictm  Griqihics  Power  Series  comes  equif^ied  with  extensive  support  for 
parallel  qiplications  at  both  the  operating  system  level  and  the  compilers  and  tools 
level.  Moving  the  applications  from  the  operating  system  level  to  the  compiler  or  tool 
level  greatly  enhances  die  ease  of  developing  parallel  applications. 

Code  can  be  parallelized  relatively  easily  using  compilers  such  as  the  Silicon 
Gri^hics  Power  FORTRAN  compiler.  This  compiler  has  a  scheme  for 
homoparallelization  of  DO  loops  through  in-line  comments  that  provide 
multiprocessing  directives  to  the  compiler.  Interpretatimi  of  these  directives  is  enabled 
through  die  use  of  the  -mp  option  with  this  cmnpiler.  The  basic  construct  is  the 
c$<hMicros5  directive  whidi  informs  die  Power  FORTRAN  compiler  that  the  DO  loop 
immediately  following  is  to  be  parallelized. 

Since  the  c$doacross  directive  can  be  cumbersome  to  insert  manually,  Powa 
FORTRAN  is  equipped  with  code  analyzer  tools  which  heh>s  the  user  identify 
opportunities  for  parallelization.  These  tools  are  called  pfa  and  pea  respectively. 

Bodi  are  source-to-source  (^timizing  processors  that  discover  blocks  of  code  that  can 
be  made  parallel  and  actually  insert  the  directives  necessary  to  implement  parallelism. 
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The  pfii  tool  was  used  for  this  task,  llie  pfa  tool  is  generally  very  good  at  identifying 
which  loops  can  be  parallelized  and  which  have  data  dependencies  or  too  few 
iterations  to  be  considered.  It  generates  a  listing  that  shows  which  portion  of  the  code 
were  parallelized  and  why  it  could  not  parallelize  the  others.  This  listing  helps  the 
develops  by  highlighting  any  dq)endencies  which  require  further  attention. 

Another  tool  that  helps  the  devel(^)er  parallelize  ^>plications  is  an  executitm 
profiler.  An  executitm  profiler  such  as  |Mxie  counts  the  number  of  CPU  cycles  a 
section  of  code  takes  to  execute  and  keeps  track  of  die  total  use  by  that  sectitm  of  the 
code.  By  accumulating  these  statistics  during  a  run  of  an  application,  this  tool  can 
generate  detailed  rc^xxts  on  the  execution  of  the  code.  An  example  of  the  pixie 
{x-ofiler  results  accumulated  from  a  run  of  the  baseline  code  q>plied  to  a  typical 
solution  is  presented  in  ^[jpendix  B.l. 

For  the  baseline  profile  example  in  A|^)endix  B.l,  the  solution  of  1000  Roe 
iterations  for  a  typical  nozzl^cowl  configuration  requires  86.82  minutes  of  CPU  time. 
The  routines  sqrt  and  rjsign  are  matii  lilnrary  routines  and  can  not  be  parallelized. 

The  remaining  routines  whidi  were  parallelized  are  FSIROE,  FSJROE,  LMTRI, 

LMTRJ,  VISEUL,  VJSEUL,  and  L.  Between  them,  diese  subroutines  todc 
q>Ix'oximately  80%  of  the  total  q>plication  run  time. 

The  pfis  tool  was  then  used  to  analyze  and  parallelize  the  code  for  each  of 
these  sulx-outines.  All  of  tiie  subroutines  excq>t  L  wne  parallelized  without  any 
manual  intervention.  With  the  help  of  pfa,  the  problem  areas  in  L  were  highlightad 
and  several  minor  changes  were  then  made  to  further  parallelize  this  subroutine. 

Details  of  these  minw  changes  in  L  and  an  example  of  the  c$doacross  structures  can 
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be  found  in  Aj^ndiz  B.3.  Since  the  other  subroutines  were  successfully  parallelized 
by  phi  without  any  manual  code  changes,  their  listings  were  not  included. 

The  executimi  profile  for  the  parallelized  code  running  this  same  nozzle/cowl 
example  is  presented  in  Appendix  B.2.  Hie  limiting  CPU  time  required  for  the  most 
heavily  used  processor  was  reduced  to  33.85  minutes.  The  speedup  ratio  in  terms  of 
CPU  cycles  was  2.56.  This  particular  test  case  took  158  minutes  with  no 
parallelization  and  69  minutes  with  parallelization,  whidi  conesponds  to  a  time 
speedup  factor  of  2.29.  The  numbers  associated  with  time  are  subject  to  some  amount 
of  variation  due  to  the  competition  for  CPU  time  with  die  operating  system. 
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m.  ANALYSIS  OF  NOZZLE  COWl.  FLOWFIELDS 


3.1  An  Experimentally  Validated  Nozzle/Cowl  Geometry 

The  first  nozzle/cowl  geometry  analyzed  in  this  study  is  based  on  the 
e]q>eiimental  hypersonic  vehicle  afterbody  ctmfiguration  analyzed  by  Cochran  (4). 

The  two-dimensitmal  nozzle/cowl  model  was  designed  to  be  representative  of  the 
external  geometry  of  a  hypersonic  vehicle.  Hgure  3-1  presents  a  cross-sectional  view 
of  this  experimental  nozzle/cowl  assembly.  In  order  to  facilitate  fabrication,  the 
internal  nozzle  consisted  of  a  SO  degree  straight  convergent  ramp.  At  the  throat,  the 
flow  is  turned  at  a  sharp  comer  and  allowed  to  expand  along  a  20  degree  external 
expansion  ramp.  Four  different  cowl  configuratimis  were  aiudyzed  in  Cochran’s 
experimental  study:  a  baseline  short  cowl,  a  long  cowl,  a  long  cowl  deflected  towards 
the  body  5  d^t^>  and  a  long  cowl  deflected  away  from  the  body  5  degrees. 

Cochran  analyzed  these  hypersonic  vehicle  afterbody  configurations  in  the 
Wright  Laboratory  two-foot  trisonic  gasdyiuunic  tuimel  at  subsonic,  trairsonic,  and 
superscmic  off-design  operating  conditions.  The  testing  was  limited  to  relatively  slow 
speed  off-design  operating  conditions  due  to  the  q;)erating  limitations  of  the 
gasdynamic  tuimel  facility.  For  the  subsonic  flow  cases  the  maximum  dynamic 
pressure  was  limited  to  350  pounds  per  square  foot  (psf)  and  the  maximum  Reynolds 
number  was  limited  to  2.5  millimi  per  foot.  For  the  supersoiuc  cases  the  maximum 
dynamic  pressure  increased  to  1000  psf  and  the  maximum  Reynolds  number  ino^eased 
to  5.0  million  per  foot.  The  total  stagnation  tempm'ature  was  maintained  at 
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i4>proximately  S59.7  ”R  fw  all  of  ttie  analysis.  The  nozzle/cowl  wall  temperatures 
were  maintained  at  a  nearly  constant  400.6  ”R  for  this  analysis. 

Cochran’s  goal  was  to  e]q>erimentally  analyze  the  nozzle/afterbody  drag  force 
on  these  models  by  integrating  the  pressure  distribution  on  the  external  ramp. 

Pressure  distributions  on  the  external  ramp  were  obtained  for  each  cowl  configuration 
over  a  range  of  NPR,  freestream  Mach  number,  and  Reynolds  number.  The  maximum 
errors  rep<»led  for  the  experimental  tests  were  0.20%  ftxr  the  hneestream  total  pressure, 
0.12%  fcv  the  freestream  static  pressure,  0.57%  for  the  static  pressure  probes,  and 
0.43%  for  the  internal  nozzle  total  pressure. 

Hyun,  (1),  compared  these  data  to  the  results  of  the  Wright  Laboratory  finite 
volume  computational  code  described  in  Chapter  2.  The  goal  of  Hyun’s  researdi  was 
to  attempt  to  validate  the  numerical  method  using  the  experimental  data  gathered  from 
this  voy  practical  experimental  example.  In  general  the  two-dimensional  flux¬ 
splitting  elgtHithms  performed  very  well  fcx*  the  supersonic  external  flow  cases. 
Three-dimensional  effects  in  the  e7q>erimental  data  wo’e  the  primary  reasons  for  most 
of  the  small  discr^ancies.  Fw  the  supersonic  external  flow  cases,  small  differences 
occurred  at  the  far  downstream  locations  on  the  expansion  ramp  where  the  thickening 
boundary  layer  tended  to  interact  with  the  external  flowrield.  Hyun  demonstrated  tbat 
as  NPR  and  freestream  Mach  increased,  the  correlation  between  the  experimental  data 
and  the  computational  results  generally  increased. 

The  two-dimensional  Wright  LaboratCMry  code  could  not  be  validated  for  the 
subsonic  flow  case  because  of  very  strong  three-dimensional  effects  on  the 
nozzle/cowl  experimental  data.  The  experimental  data  contained  significant  three- 
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AWjmional  effects  due  to  the  interactioa  between  the  exhaust  plume,  die  expansitm 
ramp  subsonic  boundary  layer,  and  the  subsonic  external  flow.  For  these  reasons,  the 
scope  of  the  cuirent  analysis  is  limited  to  the  validated  supersonic  off-design 
conditions  analyzed  by  Hyun  and  Cochran. 

The  initial  conditions  used  for  this  analysis  were  chosen  to  match  the 
experimental  conditions  reported  by  Cochran  and  used  by  Hyim.  Table  3-2  describes 
the  Mach  1.9  flowfields  analyzed  at  nozzle  pressure  ratios  of  approximately  3,  S,  7, 
and  12.  Table  3-3  describes  the  Mach  3.0  flowfields  analyzed  at  nozzle  pressure 
ratios  of  5.  7.  9,  12,  and  16.  The  test  point  number  (TPN)  presented  in  these  tables 
corresponds  to  Cochran’s  mi^nal  test  points.  For  this  experimental  nozzle 
configuration,  the  high  pressure  internal  nozzle  flow  is  subsonic  and  it  is  accelerated 
to  supersonic  speeds  tiu-ough  the  converging-diverging  nozzle  section. 

The  computational  grid  develq)ed  by  Hyun  was  used  as  a  baseline  for  the 
current  analysis.  Hyun  used  an  approach  that  is  very  similar  to  the  approach  taken  for 
the  generic  nozzle  grid  generation  described  in  Section  3.2.2.  Figure  3-2  shows  tiie 
final  grid  that  was  used  for  the  baseline  e]q>erimental  nozzle  analysis.  This  grid  has 
dimensions  of  101  (Q  by  71  (11).  An  al^nuc  distribution  of  grid  points  was  used 
along  each  edge  and  a  combinaticm  of  algebraic  and  elUptic  grid  generattxs  were  used 
to  distribute  tiie  internal  nozzle  grid  points. 

The  primary  objective  of  the  first  part  of  this  research  was  to  insert  a  gas 
generator  flowfield  into  the  cowl  lip.  Using  Hyun’s  grid  definition  and  Cochran’s 
initial  conditions,  the  overall  effects  of  the  gas  generaUx  could  be  compared  to  the 
experimentally  validated  baseline  flow  solutions.  This  also  allows  for  a  direct 
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comparism  of  gas  generator  mass  flow  and  tteflection  angle  effects  to  the  effects  of 
nozzle  cowl  length  extensions  and  deflections. 

3.2  A  Generic  Hypersonic  Vdiide  Configuration 

For  the  second  part  of  the  current  research  effort,  a  generic  hypersonic  vehicle 
nozzle/cowl  geometry  was  evaluated  at  more  realistic  on-design  flight  conditions. 
Since  the  Wright  Laboratory  code  showed  good  crxrelation  to  the  experimental  data 
and  a  trend  towards  increasing  accuracy  with  increasing  freestream  Mach  number  and 
NPR,  it  was  determined  that  the  same  code  could  be  used  to  determine  trends  at  m<»re 
realistic  on-design  flight  conditions. 

The  generic  nozzle  design  used  for  this  study  was  develq)ed  by  Doty  (9). 

This  same  nozzle  definition  was  used  in  the  nozzle  geometry  optimization  research  of 
Boni^MUte  (2).  The  purpose  of  Bonaparte’s  study  was  to  implemrat  a  two  parameter 
optimization  for  nozzle  attachment  angle  and  cowl  deflection  angle  using  an  invisdd 
flux-difference  splitting  code  devel(^)ed  by  Doty.  Doty’s  nozzle  definition  and 
Bonaparte’s  trajectory  were  used  fw  the  present  analysis.  Trends  associated  with  a 
gas  generator  flowfield  were  evaluated  for  the  generic  nozzle/afterbody  over  a  diis 
typical  flight  trajectoty  coresponding  to  a  constant  dynamic  pressure,  q,  of  1000  psf 
through  flight  Madi  numbers  ranging  between  10  and  25. 

3.2.1  Nozzie/Cowl  Modd.  There  are  many  rations  available  fo  defining  a 
generic  nozzle/cowl  assembly.  Doty,  (9),  defined  a  family  of  genoic  nozzles 
consisting  of  a  circular  arc  section  extending  from  the  combustor  exit  to  the  nozzle 
wall  and  a  parabolic  arc  section  which  defines  the  nozzle  wall.  Figure  3-3  (not  drawn 
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to  scale)  defines  this  generic  configuration.  An  internal  height,  h^  an  extonal  height, 
h^  a  nozzle  length,  and  a  nozzle  wall  attachment  angle  (Og)  «e  all  that  is  required  to 
parametrically  describe  this  family  of  nozzles.  The  circular  arc  is  defined  relative  to 
the  migin,  O,  by  the  following  relationships: 


x=JCo+ftsin(e)l 

(3-1) 

(3-2) 

where  0  extends  from  0  to  0^. 

The  general  fcxm  of  the  parabolic  wall  is  given  by: 

y^*c^x+c^+c^=0 

(3-3) 

where  the  three  unknowns,  c^,  Cj,  and  C3  can  be  evaluated  using  known  information. 

The  following  relationships  define  the  nozzle  wall: 

^  _  -[(yc^-(y^*2(y^£fiXc-Xg^ 
(yc-yj^’OCc-x^iy^ 

(3-4) 

<^i=-2(ya)(y«)-(y^2 

(3-5) 

(3-6) 

where. 

ys(4y/dbo,=tan(0^ 

(3-7) 

and  the  x  and  y  locations  are  known  at  points  A  and  C.  For  a  given  nozzle  inlet 
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height,  nozzle  exit  height,  and  nozzle  length,  a  complete  family  of  nozzle 
ctmfigurations  is  defined  by  the  nozzle  attachment  angle,  6^. 

Boni^arm’s  optimization  effort  showed  that  the  nozzle  perfcxmance  is  fairly 
sensitive  to  the  nozzle  attachment  angle.  For  this  generic  nozzle  configuration 
analyzed  over  the  ^ical  trajectory  presented  in  Sectira  3.2.3,  the  optimum  nozzle 
attachment  angle,  6b,  varied  from  nearly  36  degrees  at  a  trajectory  Mach  number  of  10 
down  to  almost  16  degrees  at  a  trajectcxy  Mach  number  of  25.  Ihe  nozzle 
configuration  chosen  for  the  current  analysis  was  based  on  a  design  traject<M7  Mach 
number  of  IS  whtfe  the  q)timum  nozzle  attachment  angle  is  q>proximately  24 
degrees.  The  Mach  IS  design  condition  and  the  24  degree  nozzle  attadunent  angle 
represent  a  reasrmable  design  trade-off  ova*  the  Mach  10  to  Mach  2S  flight  trajectory 
since  the  24  degree  nozzle  attachment  angle  falls  just  less  than  half  way  between  the 
optimum  values  for  Madt  10  and  Mach  2S. 

Table  3-1  completely  defines  the  generic  nozzle/cowl  gernnetry  used  in  this 
investigation.  The  total  e3q)ansion  ratio  is  2S:1,  the  overall  nozzle  length  is  100  times 
the  combustitm  chamber  exit  height,  and  the  cowl  wall  spans  1/10^  of  the  overall 
nozzle  length.  The  cowl  wall  is  modelled  as  a  flat  plate  with  a  blunt  trailing  edge. 
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Table  3*1  Generic  Nozzle  Configuration  Definition 


Nozzle  Data 

Value 

Nozzle  Length,  L  (in) 

100 

Inlet  Height,  h^  (in) 

1 

Exit  Height,  h^  (in) 

25 

Circular  Arc  Radius  of  Curvature,  r  (in) 

1 

Circular  Arc/Nozzle  Attachment  Angle,  6^  (deg) 

24 

Cowl  Length  (in) 

10 

Cowl  Thickness  (in) 

0.25 

3.2J2  Grid  Generation.  Successful  grid  generation  is  a  very  meticulous 
iterative  procedure.  The  process  of  develq;>ing  the  grids  used  in  this  research  was 
simplified  using  the  GRIDGEN  software  package  developed  under  contract  by  the 
Pot  Worth  Division  of  the  Lockheed  Coporation  (ftxmerly  General  Dynamics)  for 
the  Flight  Dynamics  Lalxxatmy,  Aeromechanics  Division,  Wright  Labwatory 
(WL/PIM)  (16,17). 

A  quasi  two-dimensicmal  surface  grid  is  required  to  detine  the  generic 
nozzlc/afterbody  ctmtiguration  investigated  in  this  research.  The  tinal  grid  witii 
dimensions  of  (101x71)  used  for  this  analysis  is  presented  in  Hgure  3-4.  The  grids 
used  for  the  grid  resolution  study  are  shown  in  Figures  3-5  and  3-6  corresponding  to 
grid  dimensions  of  (151x101)  and  (201x151)  respectively.  Details  of  the  grid 
resolution  study  can  be  found  in  Section  4.1. 
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Only  the  two<diniensional  grid  generedrm  portion  of  GRIDGEN  called 
GR1DGEN2D  was  required  for  this  effort.  GRIDGEN2D  is  designed  for  use  cm  the 
Silicon  Gnq>hics,  Inc.  (SGI)  3000  and  4D  series  IRIS  workstations  and  the  IBM  RISC 
6000  workstatirms.  The  GRIDGEN  software  was  installed  oa  a  NAIC  SGI  network 
for  this  effort 

For  tiiis  problem,  a  single  block  was  defined  witii  the  desired  dimensimis 

n _ t _ was  set  to  1  since  only  a  two-dimensional  surface  grid  was  required. 

£ _ and  Tt _ were  set  to  the  overall  grid  dimensions  of  101x71  respectivdy.  The 

generic  nozzk/cowl  definitimi  is  well  suited  for  grid  generation  in  Cartesian 
coordinates. 

The  edges  were  defined  by  importing  a  high  resolution  edge  definitimi  file 
which  omtained  the  nozzle  geometry  definitimi  as  specified  in  Equations  3-1  through 
3-7  and  Table  3-1.  The  grid  points  were  then  distributed  along  the  edges  using 
algebraic  distributitms. 

The  main  domain  was  broken  down  into  several  subtmits  in  order  to  he4> 
cluster  grid  points  and  control  the  grid  definition.  Grid  clustering  was  jx-ovided  in  the 
regions  near  the  solid  boundaries  to  provide  suffideat  resolution  of  the  boundary  layer 
flow  and  in  tiie  regirms  whwe  strong  e]q>ansi(ms  and  shocks  were  expected. 

The  external  flow  grid  was  calculated  algebraically  using  the  spacing  defined 
on  the  edges.  The  algebraic  solver  used  transfinite  inteipolation  (TFO  with  arclength 
based  inteipolants  to  generate  the  external  flow  grid.  This  method  tended  to  keep  the 
internal  grid  clustering  proportional  to  die  defined  edge  clustering  as  was  desired. 
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Tbe  intenial  flow  grid  was  calculated  uaing  an  elliptic  partial  differential 
equatkn  (PDE)  grid  solver.  Thomas-Middlecoff  control  functirms,  with  a  relaxation 
factor  of  0.3,  were  used  for  smoothing  since  they  ^^>eared  to  give  the  best  results. 
The  ell4>tic  solver  was  set  up  to  force  the  grid  to  be  perpendicular  to  the  solid 
boundaries  and  the  fieestream  inflow  boundary  and  was  set  up  to  force  a  constant 
gradient  geometry  at  the  downstream  outflow  boundary.  Approximately  400  iterations 
were  run  to  generate  the  final  grid. 

Once  the  surface  grid  was  created,  a  small  FORTRAN  routine  was  used  to  add 
the  requited  boundaiy  ghost  points  as  described  in  Section  2.3.  This  routine  also 
ctmvated  die  grid  to  a  cell  centered  grid  and  it  inverted  the  grid  for  evaluation  widi 
the  Wright  Laboratory  code. 

3JL3  The  Tr^|ectory.  A  hypersonic  vehicle  will  typically  travel  along  a 
design  trajectory.  A  cmistant  dynamic  {xessure,  q,  trajectory  can  be  used  to  rqxesent 
typical  aerodynamic  and  structural  loads  on  a  hyperstmic  vehicle.  A  typical  value  of 
dynamic  pressure  for  a  NASP  type  vehicle  is  g=1000  psf  (2).  The  trejecUvy  is  then 
specified  by  the  definition  of  dynamic  pressure: 

A  pressure  altitude  corresp<»ding  to  p.  is  determined  as  a  function  of  design 
dynamic  pressure  and  vehicle  Mach  number.  Figure  3-7  defines  the  trajectory 
associated  with  a  amstant  1000  psf  dynamic  pressure  flight  path.  The  pressure 
altitude  was  based  cm  the  1976  standard  atmosphere. 
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Table  3-4  defines  the  internal  and  eicternai  Pi  iwfields  used  for  this  analysis. 
The  internal  flow  parameters  were  based  on  the  engine  analysis  conducted  by 
Bonaparte,  (2).  B<ui^>arte  obtained  the  en^e  exhaust  parameters  by  analyzing  a 
generic  scramjet  configuration  over  this  trajectory  using  the  Johns  Hopkins  RamJet 
Porformanoe  Analysis  (RJPA)  {xogram  (24). 

The  external  flowfield  presented  in  Table  3-4  was  based  on  another  set  of 
simplifying  assumptions.  Since  the  forebody  of  a  hypersonic  vehicle  typically  acts  as 
a  large  compression  ramp,  assumptions  had  to  be  made  regarding  the  effect  of  this 
flowfield  compressitm  on  die  afterbody  freestream  flow.  Bonaparte  assumed  a  total  of 
eight  d^rees  of  flowfield  turning  was  accomplished  by  passing  the  external  flow 
through  an  oblique  shock.  The  oblique  shock  solution  was  calculated  using  an 
iterative  solutirm  which  accounted  for  temperature  changes  in  the  caloric  models  used 
for  oithai^y.  h.  and  ^ledfic  heat  c^iacity  at  ctmstant  pressure.  Cj,.  Although  the 
Wright  Laboratory  code  used  for  this  research  is  not  capable  of  mixing  rfissimiiiir 
flowfields,  a  realistic  freestream  boundary  txmdititm  is  desired.  The  trajectory  analysis 
completed  using  the  caloric  model  provided  more  realistic  freestream  inflow  pressure 
and  flux  terms  which  were  then  used  in  the  Wright  Laboratory  code  as  described  in 
Section  2.3. 

33  The  Gas  Generator  Model 

The  primary  goal  of  the  cunent  research  was  to  determine  the  effects  of 
inserting  a  vectorable,  throttleable  gas  generator  flow  into  the  cowl  Up.  The  gas 
generator  flowfield  was  obtained  using  a  computer  program  developed  by  the  authix  to 
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determine  the  equilibrium  flow  sdutka  for  a  liquid  hydrageoHnygea  rocket  motor 
(22).  The  exit  stage  sdutioo  from  diis  model  was  used  as  the  iiqmt  initial  conditions 
for  the  improved  cowl  concept  gas  generator. 

Table  3-S  d^nes  die  gas  generator  flow  used  for  this  analysis.  The  baseline 
nozzle  geometry  and  operating  conditions  were  based  largely  on  a  scaled  down  Space 
Shuttle  Main  &igine  (SSldB).  The  oxygen  to  fuel  mixture  ratio  was  6.026:1,  the 
combustion  chamber  inlet  pressure  at  100%  throttle  was  204.SS  atm,  and  the  exit  to 
throat  area  ratio  was  773  for  the  baseline  gas  generator.  In  order  to  analyze  the 
effects  of  varying  mass  flows  and  exit  jnessures,  dus  engine  configuration  was 
analyzed  operating  at  a  range  of  throttle  settings  varying  from  25%  to  200%.  Engine 
throttliug  is  accomplished  by  increasing  or  decreasing  the  combustor  chamber  pressure. 
In  addition  to  tbrottlihg,  die  gas  gcoenux  flow  was  also  vectored  throug)i  a  range  of 
4>20  to  -20  d^ees  using  the  qifxoach  described  in  Section  2.3. 
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TaUe  3-2  Eaq)erimental  Nozzle,  Mach  1.9  Flow  Definitioii 
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Table  34  Generic  Noezle  Plow  Definitioa 


Table  3-5  Gas  Generator  Plow  Definition 


Figure  3-1  Experimental  Nozzle  Geometry  (1) 
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Figure  3-3  Generic  Nozzle  Wall  Definitioa  (9) 
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30.0  {-  Oeneric  Nozzle/Cowl  Cell  Centered  Grid 
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Figure  3-4  Generic  Nozzle/Cowl  Coarse  Cell  Centered  Grid  (101x71) 
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Generic  Nozzle/Cowl  Medium  Cell  Centered  Grid  (IS 1x101) 
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Figure  3-6  Generic  Nozzle/Cowl  Pine  Cell  Centered  Grid  (201x151) 


Typical  Trajectory 


Figure  3-7  Generic  Nozzle  Trajectory 
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IV.  RESULTS  AND  DISCUSSIONS 


The  two-dimensional  Navier  Stdces  equatimis  were  evaluated  using  the  Wri^t 
Laboratory  flux-splitting  finite  volume  code  for  two  sq>arate  geometries  over  a  wide 
raqge  of  operating  conditions.  The  nozzle  operating  conditions  analyzed  were  all 
initially  undere]q>anded  since  the  total  nozzle  exit  pressure,  Pn.  was  always 
significantly  higher  than  the  ambimit  static  pressure,  P«.  The  flowfield  was  dominated 
by  tile  initial  eiqiansion  of  the  flow  from  tins  umleiexpanded  state. 

A  schematic  of  a  typical  overmqianded  and  undermqianded  exhaust  is  [xesented 
in  Figure  4-1.  The  wall  pressure  distribution  is  typically  dominated  by  two  sq>arati<m 
rej^ons.  The  first  and  largest  separaticm  Inibble  occurs  just  past  the  nozzle  throat 
where  the  flow  undergoes  a  dramatic  initial  mqiansion.  The  secmid  aqiaratimi  region 
may  qipear  after  tiie  shock  wave  associated  witii  an  overmqianded  flow  impinges  mi 
the  nozzle  wall  boundary  layer.  These  separation  re^ons  are  characterized  by  a 
recompressimi  of  the  flow  and  a  subseqimnt  increase  in  nozzle  wall  pressure  followed 
by  a  re-eiqiansimi  of  the  flow  do^^  .  i  s  i  «am. 

For  tile  mqierimeatal  nozzle  evaluated  at  tiie  relatively  low  qpeed  off-design 
cooditimis,  the  flow  was  initially  undermqianded.  The  dramatic  initial  mqiansimi 
caused  by  tiie  sudden  increase  in  nozzle  area  typically  resulted  in  an  overexpanded 
flow  just  downstream  of  the  nozzle  throat.  The  location  of  the  second  sqiaration 
region  associated  with  the  lar^  ovmeiqiansion  and  its  associated  recompression  shock 
is  highly  dependent  on  the  freestream  Mach  numbo*,  the  NPR,  and  gas  generator  setting. 
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For  Uie  generic  nozzle  evaluated  over  the  more  realiatic  trejeciory,  die  initial 
flow  was  also  significantly  underexi>anded.  For  this  nozzle,  the  sudden  initial 
expansion  was  also  very  severe  but  the  flow  did  not  overexpand  to  a  pressure  less  than 
ambient  and  the  flow  typically  remained  undereiqMnded  throughout  die  nozzle.  There 
was  typically  no  characteristic  downstream  sqiaration  bubble  assodated  widi  an 
overeiqiansion  shodr  wave. 

Assessing  the  im(Mct  of  a  ^  geimrator  on  these  characteristic  flow  sdudoos 
was  the  immary  objective  of  diis  research.  The  overall  impact  of  adding  a  gas 
generator  flowfield,  increasing  die  gas  generator  mass  flow,  and  deflecting  die  gas 
generator  thrust  are  presented  in  the  following  sections.  Sectirm  4.2  discusses  die 
effects  of  the  gas  generator  flowfield  on  the  experimentally  validated  nozzle/cowl 
geometry  presented  in  Rgure  3-1  and  Section  4.3  presents  the  effects  of  the  gas 
generator  on  die  generic  nozzle/cowl  configuration  presented  in  Figure  3-2. 

The  flow  solutions  were  evaluated  with  the  Flow  Analysis  Software  Toolkit 
(FAST)  developed  by  NASA  Ames  Research  Center  (IS).  This  software  was  also 
installed  on  the  NAIC  SGI  network  for  this  research. 

4.1  Grid  Resolution  Effects 

The  effects  of  grid  density  must  be  evaluated  before  the  flowfield  can  be 
analyzed  with  c<Mifidence.  Any  finite  difference  or  finite  volume  scheme  will  be 
affected  by  the  truncation  errors  assodated  with  the  chosen  numerical  mediod  (10).  A 
grid  resolution  study  must  be  completed  to  assure  that  the  truncation  errors  do  not 
build  up  and  dominate  the  flow  solution.  Furthermcxe,  the  grid  must  be  dense  enough 
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to  cqjture  all  ol  die  critical  flow  stnicturea.  For  die  coofigiiratioiia  analyzed  in  thia 
study,  a  series  of  grid  refinements  was  made.  The  effects  of  a  coarse  grid  (101x71),  a 
medium  grid  (ISlxlOl),  and  a  fine  grid  (201x151)  were  evaluated  at  sevotd  opaxdng 
conditimis. 

Hyun  (1)  completed  this  grid  resolutioa  study  for  the  eaqierimental  nozzle/cowl 
configuration  jxesented  in  Section  3.1  and  Figure  3-1  at  the  H.=  1.9,  NPRs7.0,  and 
M.=i3.0,  NPR=16.0  test  omditimis.  Hyun  showed  that  the  dominant  effects  of  grid 
resolution  occurred  in  the  region  of  the  first  sqiaratioii  bubble  just  downstream  of  die 
sharp  throat.  As  grid  density  increased,  the  flow  characteristics  in  the  initial 
reoompression  r^on  were  resolved.  Grid  refinement  did  not  significandy  change  die 
magnitude  of  die  peak  reconqxession  and  it  had  almost  no  effect  on  die  re-mqiansion 
of  die  flow  past  the  first  separation  bubble.  Even  as  die  grid  resolutimi  was  increased 
and  die  recomixession  structure  was  resolved,  diere  was  still  a  noticeable  diffuence 
between  the  overall  recompression  calculated  numerically  and  die  eiqioimental 
measurements.  The  primary  reasrm  for  the  small  difference  between  the  experimental 
and  numerical  solution  was  attributed  to  the  laminar  flow  a8sun^>ti(m  used  in  the 
numerical  mediod.  Overall,  die  agreement  between  die  mqierimental  data  and  the 
numerical  sdution  was  very  good.  Since  the  results  were  not  significandy  imfxoved 
and  the  penalty  paid  in  overall  conqiutation  time  was  so  large  for  the  finer  grids,  the 
coarse  grid  was  determined  to  be  acceptable  for  this  analysis. 

A  similar  study  was  done  for  the  generic  nozzl^cowl  configuration  presented 
in  Section  3.2,  and  Figure  3-2  at  the  design  Mach  15  condition  with  the  gas  generator 
operating  at  100%  throtde  with  no  deflection.  The  three  grids  used  for  this  analysis 


60 


were  presented  in  Figures  3-3  diroug^  3-S.  Figure  4-2  shows  die  efifects  of  grid 
resolution.  As  with  the  previous  geometry,  the  Isrgest  differences  occur  in  the  r^on 
ai  the  initial  recompression.  For  this  case,  the  most  significant  changes  occur  near  die 
peak  of  the  reccmipression  and  they  tend  to  disappear  during  die  subsequent 
downstream  eiqiansion.  The  effect  of  dm  comixession  waves  associated  with  the  wall 
curvature  were  resolved  with  the  finer  grids  in  this  region. 

Since  there  is  no  eiqierimental  data  to  compare  to,  the  grid  resolution  effects 
had  to  be  evaluated  in  terms  of  their  relative  impact  (» this  gas  generator  study. 

Figure  4-3  shows  the  net  effect  of  turning  off  die  gas  generator  mi  the  nozzle  wall 
pressure  distribution  at  this  same  Mach  15  trajectory  point  using  the  coarse  (101x71) 
grid.  Note  that  the  gas  generator  is  located  in  the  cowl  at  a  downstream  location  of  10 
indies.  Since  the  flow  is  supersonic,  thoe  is  no  iqistream  effect  of  the  gas  generator. 
Furthermore,  any  pressure  differences  associated  with  the  gas  generator  must  {xopagate 
through  the  internal  nozzle  flowfield  befme  diey  can  have  an  impact  mi  die  nozzle 
wall.  Since  the  internal  nozzle  flow  is  expanded  to  fairly  hi^  Mach  numbers  for  diia 
geometry,  the  propagation  of  the  gas  geimratm  effects  are  limited  to  the  downstream 
nozzle  wall  locations. 

The  cmivergenoe  of  die  final  solution  is  determined  achieving  a  steady  RMS 
heat  transfer  and  pressure  along  the  nozzle  wall  as  defined  in  Sectimi  2.5  and 
Equations  2-41  and  2-42.  Figure  4-4  shows  the  explicit  cmivergence  history  of  die 
Madi  15  Irajectmy  point  using  the  coarse  101x71  grid.  Fm  this  example,  convergence 
fm  the  Van  Leer  solution  was  achieved  in  a  relatively  quick  500  itmations  and  an 
acceptable  convergence  for  the  Roe  solution  was  achieved  by  920  iterations. 
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For  this  same  example,  convergence  was  achieved  in  1120  iterations  for  die 
medium  grid  and  1300  iteradons  for  the  fine  grid.  The  code  required  apf^oximately 
787.9  microsecmids  per  itcraticMi  per  node  for  the  coarse  grid,  954.43  miaosecmids  per 
iteration  per  node  for  the  medium  grid  and  929.6  microseconds  per  iteration  per  node 
for  die  fine  grid.  Using  a  Silicon  Graphics  Iris  Indigo  workstation,  the  total  CPU  time 
required  for  this  example  solution  was  1.444  hours  with  the  coarse  grid,  4.529  hours 
with  the  medium  grid  and  10.188  hours  widi  the  fine  grid. 

As  illustrated  in  Figure  4-2,  the  initial  recompression  was  not  accurately 
captured  with  the  coarse  grid.  The  magnitude  o£  the  peak  recompression  tends  to 
increase  slightly  widi  increased  grid  resolution.  However,  Hyun  (1)  noted  that  the 
numerical  vs.  esqierifflental  discrepancy  in  diis  r^on  was  due  as  much  to  the  laminar 
flow  assumption  as  it  was  to  grid  resolutitm.  Note  from  Figure  4-3  that  the  cowl  and 
gas  generator  effects  are  not  evident  in  the  flowfield  nor  nozzle  {xessure  until 
significandy  downstream  of  the  initial  expansimi  r^on.  Additionally,  since  the 
downstream  influence  of  die  cowl  is  not  affected  significandy  by  the  grid  resolutiem, 
the  computatimial  time  penalty  associated  with  the  increased  grid  dmisity  was  not 
justified.  Ther^ore,  die  coarse  grid  (101x71)  was  used  for  the  remainder  of  this 
analysis. 
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Figure  4-1  Schonatic  of  Over  and  Underexpanded  Exhaust  at  Siq)arsonic  Freestream 
Mach  Number  (1) 


Grid  Resolution  Effects 


F%ure  4>2  Grid  Resolution  ^ects.  Gas  Genorator  On  100% 


Gas  GeneraUMT  Effects 


Figure  4-3  Gas  Generahx'  Effects,  Mach  15,  Coarse  Grid 
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Explicit  Stability  Check 
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F^re  4>4  Explicit  Stability  Check,  Mach  15,  Coarse  Grid  (101x71) 
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4.2  Application  of  a  Cowl  Gas  Generator  Flowfield  to  an  Experimentally 
Validated  Off-Design  Configuration 

The  experimental  hypersonic  nozzl^cowl  geometry  presented  in  Figure  3-1  was 
evaluated  with  and  without  the  gas  generator  for  two  off-design  supersonic  flow 
conditions;  Mach  1.9  and  Mach  3.0.  The  basic  gas  generator  effects  were  analyzed 
over  a  range  of  NPR  at  these  two  Mach  numbers.  The  gas  generatm-  mass  flow  and 
deflection  angle  effects  were  evaluated  fot  two  typical  conffgurations,  Mach  1.9, 
NPR=7,  and  Mach  3  0,  NPR=16  since  these  were  the  same  conditions  analyzed  by 
Cochran  (4)  and  Hyun  (1)  for  their  cowl  geometry  studies  described  in  Section  3.1. 
These  baseline  flow  conditions  correspond  to  Cochran’s  test  point  numbers  148  and 
184,  respectively. 

The  nozzle  wall  pressure  distributions  for  the  Mach  1.9,  NPRs7,  TPN  148  case 
are  presented  in  Figure  4-S  fcx  both  the  baseline  gas  geirerator  off  condition  and  tiie 
gas  generator  100%/no  deflection  operating  condition.  The  nozzle  throat  is  located  at 
a  downstream  location  of  3.25  inches.  A  is  obvious  that  the  internal  flow  upstream  of 
the  throat  is  significantly  undererqumded.  As  the  flow  passes  the  sharp  comer  at  tite 
throat,  it  experiences  a  significant  increase  in  flow  area  resulting  in  a  dramatic 
expansion.  The  flow  actually  overexpands  to  a  pressure  less  than  ambient.  Both  with 
and  without  the  gas  generator,  the  flow  immediately  begins  to  recompress  due  to  tire 
presence  of  the  first  separation  bubble  past  the  sharp  throat.  Fm-  the  baseline  case 
with  the  gas  generatm*  off,  the  recompression  is  not  very  strong  and  the  flow  remains 
overexpanded  along  the  entire  nozzle  wall  downstream  of  the  throat.  With  the  gas 
generator  on,  the  size  of  the  first  separation  bubble  increases  signiffcantly  as  indicated 
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by  the  increased  region  of  recom(Mres8i(m  along  the  nozzle  wall.  The  gas  generator  not 
only  inaeases  the  magnitude  of  the  recom{Hession,  but  it  also  moves  the  location  of 
the  peak  magnitude  and  the  subsequent  re-expansion  downstream.  The  relatively  high 
pressure  gas  generaux'  acts  like  an  aerodynamic  cowl  extension.  The  relatively  high 
pressure  boundary  between  the  internal  and  external  flowfields  created  by  die  gas 
generator  reduces  the  expansion  of  the  hig^  pressure  internal  nozzle  flow  near  the 
cowl  waU  to  the  relatively  low  pressure  external  freestream  flow.  The  effects  of  this 
increased  pressure  then  prt^agate  through  die  relatively  low  speed  supersonic  flow  just 
downstream  of  the  nozzle  throat  where  it  appears  as  an  increase  in  pressure  along  the 
nozzle  wall. 

The  shock  wave  associated  with  the  overeiqianded  baseline  case  impinges  the 
nozzle  wall  downstream  of  the  first  separation  bubble.  This  shock  impingement 
causes  a  second  smaller  sqiaration  region  to  f(xm  and  the  flow  is  again  recompressed. 
This  flow  does  not  reattach  to  the  nozzle  wall  and  the  flow  remains  separated  to  the 
end  of  the  nozzle.  With  the  gas  generator  on,  the  second  sqiaratimi  does  not  fonn. 

For  both  cases,  the  flow  overexpands  and  never  recovers  to  ambient  pressure. 

The  nozzle  wall  pressure  distribution  fcx  the  Mach  3.0,  NPR=16,  TPN  184  case 
is  presented  in  Figure  4-6.  Fw  this  increased  nozzle  pressure  ratio  and  external  flow 
Mach  number  case,  the  initial  expansion  past  the  sharp  throat  is  not  as  strong  as  the 
previous  case  and  the  flow  does  not  overexpand  past  the  throat.  The  magnitude  of  the 
first  sq>aration  bubble  increases  fcx  both  the  gas  generator  on  and  off  as  indicated  by 
the  increased  pressure  recovery  along  the  nozzle  wall.  Again,  the  gas  generator  tends 
to  increase  the  overall  size  of  the  first  separation  bubble  resulting  in  a  delayed  re- 
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oqMuision.  The  seccmd  8q>aratioa  does  not  form  until  the  very  end  of  die  nozzle  wall 
for  the  gas  generator  off  case  and  therefore,  the  downstream  pressure  for  the  gas 
generator  on  case  remains  higher  than  the  gas  generator  off  case  throughout  die  entire 
length  of  the  nozzle.  As  with  the  previous  case  neitfaer  of  these  flows  recover  to 
ambient  pressure. 

^pendix  C  presents  the  pressure  and  Mach  number  contours  for  these  two 
baseline  configurations  with  the  gas  geno-ator  operating  at  100%  throttle  with  no 
deflection  and  with  the  gas  generator  off.  Figures  C-1  and  C>2  ixesent  the  flowfield 
solution  fcM*  the  baseline  Mach  1.9,  NPR=7.0.  gas  generator  off  case.  In  addition  to 
the  nozzle  wall  eiqiaiision  and  recompression  associated  with  the  sqiaration  regions, 
the  other  flow  phenomena  can  be  seen.  Figure  C-1  clearly  dommstrates  how  die 
internal  nozzle  flow  is  dominated  by  the  initial  eiqiansimi  just  downstream  of  die 
throat  As  the  high  fxessure  internal  flow  eiqiands  past  the  cowl,  a  plume  shodc 
appears  as  a  pressure  rise  in  the  extonal  flowfield.  The  boundary  between  the  internal 
and  external  flowfield  can  also  be  seen  as  lai^  Mach  number  gradioits  in  Figure  C-2. 
In  addition,  the  internal  shock  extending  from  the  cowl  wall  and  die  impingement  of 
diis  shock  on  the  nozzle  wall  is  also  evident  in  Figure  C-2.  At  qipraximately  6.S 
inches  downstream,  the  internal  shock  impinges  on  the  wall  and  the  adverse  {xessure 
gradient  acting  on  the  viscous  boundary  layer  causes  a  flow  sqiaraticm  whidi  can  also 
be  seen  in  figure  C-2. 

Figures  C-3  and  C-4  show  this  same  cases  with  the  gas  generator  operating  at 
100%  throttle  with  no  deflection.  The  interaction  of  the  nozzle  flow  structures  with 
the  gas  generate  flow  can  now  be  seen.  Figure  C-3  shows  how  the  gas  generator  acts 
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•8  an  aerodynamic  cowl  extenaioo  by  delaying  the  large  initial  expansioo  just 
downstream  of  the  cowl  wall.  The  increased  pressure  recovery  along  the  nozzle  wall 
is  also  evident  in  this  figure.  Rguie  C-4  clearly  shows  the  location  of  the  gas 
generator  boundaries  as  very  strong  Mach  number  gradients.  Furdrermcx'e,  Figure  C-4 
shows  how  the  internal  nozzle  shock,  the  shock  impingement  <mi  the  nozzle  wall,  and 
the  subsequent  flow  sq>aration  associated  with  an  initial  overexpansion  does  not  form 
with  the  gas  generator  on.  These  figures  also  show  that  there  is  a  slight  increase  in 
the  external  flow  plume  shock  strength  dim  to  the  expansion  of  the  high  faessure  gas 
generauv  boundary  into  the  external  flowfield. 

Figures  C-S  through  C-8  show  the  same  information  for  the  Mach  3.0, 
NPR=16.0  case.  Except  for  the  decreased  plume  shodc  angle  associated  widi  die 
higher  external  flow  Mach  number,  and  the  increased  internal  nozzle  flow  eiqMnsioii 
associated  with  the  higher  nozzle  (xessure  ratio  operating  conditimis,  the  same  basic 
trends  are  obvious.  A  comparisixi  of  Hgures  C-4  and  C-8  show  the  effect  of  nozzle 
(^lerating  conditim  <m  the  gas  generator  plume  sluqie.  The  plume  is  pushed  out 
significantly  due  to  the  higher  internal  nozzle  i»essure  assodated  with  the  increased 
Madi  number  and  increased  nozzle  pressure  ratio  operating  condititm. 

An  i^preciation  for  the  basic  effects  of  the  gas  generator  can  be  gained  by 
investigating  a  cross  section  of  the  flowfield  just  downstream  of  the  gas  generator. 
Figures  4-7  and  4-8  show  the  pressure  distribution  across  the  height  of  the  nozzle  at  a 
constant  axial  location,  indicated  by  index  T',  just  downstream  of  the  cowl  edge.  A 
constant  location  of  1=52,  corresponding  to  an  approximate  downstream  location  of 
3.27  inches,  was  used  to  generate  these  profiles.  For  reference,  note  that  the  edge  of 


the  cowl  and  die  gas  generator  is  located  between  2.375  and  2.S  inches.  Pdr  both  of 
these  cases,  the  flow  along  die  nozzle  wall  (hsl.914  inches)  has  already  undergone  a 
dramatic  eiqiansion  due  to  the  rapid  increase  in  flow  area.  Hie  internal  nozzle  flow 
that  was  traveling  along  the  cowl  wall  (h  <  2.375)  is  still  at  relatively  high  pressure 
relative  to  the  ambient  pressure,  (h  >  2.5).  For  these  operating  conditions,  the  gas 
generator  opo'ates  with  an  exit  pressure  that  is  higher  than  the  expanded  cowl  base 
jxessure  associated  with  the  gas  generat(x  off  cmidition.  This  higher  pressure  behind 
the  cowl  tends  to  keq>  the  internal  nozzle  pressure  higher  and  ddays  the  expansiw  of 
this  flow  to  the  external  freestream  conditions.  This  higher  pressure  dien  propagates 
through  the  nozzle  resulting  in  a  higher  pressure  distribution  on  the  downstream  nozzle 
wall. 


4.2.1  Gas  Generator  Mass  Flow  Effects.  The  effects  of  varying  the  gas 
gmierator  mass  flow  were  investigated  fcx-  the  Madi  1.9  and  Mach  3.0  flow  ccmditimis. 
The  gas  generator  throttle  was  increased  from  0%  to  200%  in  25%  increments.  For  a 
fixed  geometry  gas  generator,  varying  the  dnotde  efiectively  changes  die  gas  generator 
exhaust  (xessure  and  mass  flow  rate.  Figure  4-9  presents  the  dfect  of  varyiqg  the  gas 
generator  mass  flow  rate  on  the  nozzle  wall  pressure  distributicm  for  die  Nfach  1.9, 
NPRb7,  TPN  148  case.  As  the  throtde  setting  is  increased,  die  overall  size  of  the 
initial  sqiaration  bubble  increases  as  indicated  by  an  increase  in  recompressimi  along 
the  nozzle  wall.  In  general,  increasing  the  gas  g^ierator  mass  flow  rate  and  exhaust 
pressure  increases  the  initial  rate  of  recompression  and  moves  the  location  of  the  peak 
recompression  slightly  upstream.  The  subsequent  downstream  re-expansi(m  is  also 
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deUyed  by  the  increue  in  the  size  of  the  separstioa  bubble.  At  the  lower  mass  flow 
rates,  the  effect  of  an  e]q>ansion  and  subsequent  recompressim  can  be  seen  within  the 
first  sqparation  bubble.  In  general,  as  the  gas  generator  mass  flow  is  increased,  the 
effects  of  the  sec<Hid  sqjaration  region  disappear  resulting  in  a  decreased  pressure  near 
the  end  of  the  nozzle  wall. 

The  nozzle  wall  pressures  were  integrated  using  Equations  2-44  to  calculate  the 
nozzle  wall  contributimi  to  thrust.  Figure  4-10  ]xesents  the  effects  of  varying  the  gas 
geuf-intor  mass  flow  in  terms  of  the  nozzle  wall  pressure  contributions  to  thrust  The 
thrust  values  are  negative  for  this  configuration  indicating  there  is  a  net  pressure  drag 
acting  on  the  nozzle  wall  for  this  operating  condition.  This  is  caused  by  the  domiiumt 
high  (vessute  on  the  internal  SO  degree  converging  wall.  This  figure  highlights  diat  at 
this  operating  condition,  simply  placing  a  gas  generator  into  the  cowl  wall  gives  a 
larger  net  effect  than  varying  the  mass  flow.  Increasing  mass  flow  toads  to  iiKrease 
the  performance  due  to  the  rise  in  pressure  associated  with  the  huger  initial 
recmnpression  region.  Decreasing  mass  flow  does  not  appear  to  significantly  decrease 
the  performance  at  diese  operating  comlitions.  This  is  caused  by  fee  slight  increase  in 
downstream  nozzle  wall  fuessure  feat  occurs  when  the  gas  generator  is  operating  at 
low  ferotde  settings. 

Figures  4-11  and  4-12  show  the  same  information  presented  above  for  the 
increased  Mach  and  NPR  case;  Mach  3.0,  NPR=16.0,  ITN  184.  These  figures  show 
fee  same  basic  trends  that  were  observed  for  fee  Mach  1.9,  NPR=7.0,  TPN  148  case 
described  above.  Since  this  case  did  not  have  a  significant  second  8q>aration  region, 
the  downstream  recomiu'ession  associated  with  it  did  n(^  significantly  affect  the  trends 
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•t  tibe  low  throttle  settiqgs.  The  chenge  in  the  int^ated  pressure  dirust  ctmtribution 
varies  nearly  linearly  wiA  changes  in  gas  generator  throttle  setting  at  this  operating 

OMiditirxi. 

4.2J  Gas  Generator  Deflection  Angle  Effects.  The  effect  of  deflecting  the 
gas  generator  flow  was  also  investigated  for  the  Mach  1.9  and  Mach  3.0  flow 
conditiiMis.  The  gas  generator  deflection  angle  was  varied  fnnn  -10  to  -i-lO  degrees  in 
2.5  degree  inaenients.  Figure  4-13  shows  the  dfect  of  deflecting  the  gas  generator 
for  the  Madi  1.9,  NPRs7.0,  TPN  148  case.  The  gas  generator  was  held  at  a  constant 
100%  throttle  setting  for  all  of  these  profiles.  The  overall  trends  associated  with 
deflecting  the  gas  generator  flow  are  quite  different  than  the  trends  associated  with  die 
change  in  mass  flow.  Like  the  increasing  mass  flow  trends,  as  die  flow  is  deflected 
towards  the  nozzle  wall  (negative  deuection  angte),  the  rate  of  recomi»e8si(»  increases 
and  the  peak  recmnpression  increases  and  moves  upstream.  Unlike  die  increasing 
mass  flow  trends,  as  the  flow  is  deflected  towards  the  nozzle  wall,  the  initial 
recompression  occurs  much  faster  and  the  magnitude  of  the  peak  cmnixessioo  is  much 
larger.  The  effects  of  gas  generator  deflectimi  are  primarily  limited  to  the  vicinity  of 
die  first  sqiaraticm  bubble.  The  downstream  jnessure  distributions  lend  to  converge  to 
the  same  solution  calculated  for  the  n<»  deflected  gas  generator. 

Figure  4-14  presents  the  gas  generator  deflection  aqgle  effects  on  the  nozzle 
wall  {uessure  ccmtiibutitm  to  thrust  Sharp  increases  in  the  peak  recmnpression  with 
ddlections  towards  the  wall  tend  to  increasingly  improve  the  {Hessure  contribution  to 
thrust.  Ddflections  rway  from  die  wall  decrease  the  peak  compression  and  decrease 
the  performance  improvement  associated  with  the  gas  generator.  The  performance 
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decreuet  with  deftocdcMu  twty  from  the  wall  are  not  as  great  as  the  poformanoe 
increases  associated  with  deflections  towards  die  wall.  This  occurs  because  the 
deflections  towards  the  wall  increase  both  die  magnitude  of  the  peak  compression  and 
the  rate  at  which  die  recompression  occurs.  Deflections  away  from  the  wall,  tm  the 
other  hand,  tend  to  decrease  die  peak  magnitude  of  the  recompression  but  diey  do  not 
signiflcantly  decrease  the  rate  of  initial  recompression. 

Figures  4-15  and  4-16  show  the  gas  generator  deflection  angle  effects  for  the 
Mach  3.0,  NPR=16.0,  TPN  184  case.  The  overall  trends  are  similar  to  the  trends 
observed  for  the  Mach  1.9  operating  condition.  The  gas  generabx  deflecticm  effects 
are  characterized  by  increases  in  the  magnitude  of  the  peak  recompression  and  the  rale 
of  initial  recomi»res8i<Mi. 

4  J  J  Gas  Generator  Comparison  to  Nozzle/Cowl  Geometry  Effects. 

Cochran  (4)  and  Hyun  (1)  evaluated  the  effects  of  modifying  the  cowl  edge  gemnelry 
for  this  conflguraticm  at  these  same  two  operating  conditicms.  An  extended  cowl  and 
an  extended  cowl  with  +  5  and  -5  d^ree  deflections  were  evaluated  experimentally  by 
Cochran  and  numerically  by  Hyun.  Figure  4-17  shows  the  results  of  Hyun’s 
numoical  analysis  compared  to  tiie  baseline  gas  generator  flow  case  for  the  Madi  1.9, 
NPRs7.0,  TPN  148  case.  The  basic  effect  of  a  cowl  wall  extension  is  very  similar  to 
the  effect  of  placing  a  gas  generator  into  tiie  baseline  cowl.  Increasing  the  cowl  wall 
length  increases  the  overall  size  of  tiie  first  sqiaration  bubble.  Unlike  the  gas 
generator,  however,  the  effect  of  extending  the  cowl  does  not  minimize  the  effects  of 
the  second  sqiaratitm. 
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Overall,  the  basic  effects  of  cowl  wall  deflections  are  very  similar  to  gas 
generator  deflecti<n  effects.  As  the  cowl  wall  extension  is  delected  towards  the  wall, 
the  peak  recompression  increases  and  moves  upstream.  Like  a  gas  generator 
deflection,  the  overall  effect  of  the  cowl  extension  deflectimi  is  limited  to  the  first 
8q>arati<»i  bubble  and  does  not  significantly  alttf  the  downstream  flow. 

For  the  conditions  analyzed,  the  gas  generaUx  behaves  like  an  aerodynamic 
cowl  extension,  hi  general,  gas  generator  mass  flow  and  exhaust  (xessure  inoeases 
can  be  used  to  simulate  the  effects  of  cowl  extensions.  Furthermore,  gas  generator 
deflections  can  be  used  to  simulate  the  effects  of  mechanical  cowl  deflections. 

Figure  4-18  compares  the  available  cowl  gemnetry  data  to  the  baseline  gas 
generator  flow  soluticm  for  the  Mach  3.0,  NPR=16.0,  TPN  184  case.  The 
eiqperifflental  and  numerical  data  was  only  available  for  the  -5  degree  deflection 
towards  the  nozzle  wall  The  overall  trends  are  very  similar  to  die  jxevious  case 
because  the  gas  generator  still  provides  a  relatively  high  pressure  boundary  between 
the  internal  and  extanal  flowfields  which  acts  like  an  aerodynamic  cowl  extension. 

4.2.4  Nozzle  Pressure  Ratio  Effects.  The  effect  of  varying  nozzle  pressure 
ratio  was  also  investigated.  Figure  4-19  presents  the  nozzle  jxessure  ratio  effects  for 
the  Mach  1.9  operating  conditicms.  For  this  case,  the  nozzle  pressure  ratio  was 
investigated  at  NFR=  3.0,  5.0,  7.0,  and  12.0.  As  the  nozzle  pressure  ratio  increases, 
the  magnitude  of  the  filrst  sqiaration  bubble  increases  because  tiie  higher  nozzle 
exhaust  fxessure  propagates  downstream  resulting  in  a  delayed  eiqiansion  and  a  higher 
jxessure  along  the  nozzle  wall.  The  rate  of  recompression  is  increased  and  the  size  of 
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die  s^Mtfatioo  bubble  is  lengdieued,  resulting  in  t  delayed  eiqMnsioo.  The  locatioa  of 
the  peak  rectmipression  in  the  first  sq>aration  bubble  tends  to  remain  at  a  constant 
location  (3.7  inches  downstream).  The  location  of  die  second  separation  is  affected 
dramatically  by  the  operating  nozzle  pressure  ratio.  As  the  nozzle  pressure  ratio  is 
increased,  the  cowl  shock  angle  associated  with  die  initial  overejqiansion  decreases 
because  the  magnitude  of  the  initial  overespansimi  decreases.  This  deoeased  shock 
angle  moves  the  location  of  the  shock  impingement  on  the  nozzle  wall  downstream. 
Since  the  secoiul  sqiaradon  is  caused  by  the  impingement  oS  the  cowl  shock  mi  the 
nozzle  wall  boundary  layer,  the  resultant  second  sqiaration  moves  downstream.  The 
flow  remains  separated  from  the  nozzle  wall  after  the  shock  impingement  and 
therefoie,  it  maintains  a  relatively  cmistant  pressure  downstream  of  die  sqiaratimi. 

Hgure  4-20  shows  the  Mach  1.9  flow  case  analyzed  with  the  gas  gmierator 
operating  at  100%  throtde.  The  gas  generator  shows  the  same  basic  trends  diroitgbout 
the  range  of  nozzle  pressure  ratio.  The  magnitude  of  the  tirst  separation  bubble  is 
consistently  increased.  The  gas  generator  also  tends  to  minimize  the  secmid  sqiaratioa 
over  this  range  of  coating  cmiditions.  At  the  lower  nozzle  i»ressure  ratios,  however, 
a  second  sqiaraticm  does  af^iear.  The  location  of  the  secmid  sqiaration  is  typically 
determined  by  the  magnitude  of  the  recompression  shock  wave  fumed  off  the  edge  of 
the  cowl  wall.  With  the  gas  guieratu-  (gating,  tlK  overeiqiansion  past  the  cowl  wall 
and  its  associated  shodc  wave  is  eliminated.  The  delayed  second  separation  region 
which  i^pears  at  the  low  nozzle  pressure  ratio  operating  conditions  is  not  caused  by  a 
cowl  wall  recompression  shock,  rather  it  is  formed  by  the  impingement  of  a  reflected 
shock.  The  shock  wave  generated  by  the  recompression  of  the  first  separation  bubble 
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it  reflected  off  die  gtt  generttor  stream  boundary  and  it  impinges  the  nozzle  wall  at 
the  far  downstream  locadmis  for  the  low  nozzle  pressure  ratio  cases.  For  die  higho* 
nozzle  i»essurc  ratio  cases,  the  internal  nozzle  flow  eiqiands  to  higher  Mach  numbers 
and  therefore,  the  internal  and  reflected  shocks  are  more  oblique.  Subsequently,  the 
reflected  shock  does  not  impinge  (»  the  nozzle  wall  for  these  operating  conditions. 

Figures  4-21  and  4-22  show  the  same  basic  information  for  the  Mach  3.0 
freestream  flow  conditions.  For  this  analysis,  the  nozzle  pressure  ratio  was  evaluated 
at  NPR  s  S.O,  7.0,  9.0,  12.0,  and  16.0.  The  overall  nozzle  pressure  ratio  trends  are 
very  similar  to  the  Mach  1.9  flow  trends.  Since  the  freestream  Mach  was  increased 
from  1.9  to  3.0,  the  overall  cowl  shock  angle  decreased,  moving  die  second  separation 
downstream.  In  addition,  the  overall  effects  of  nozzle  pressure  ratio  on  the  first 
sqMration  bubble  were  not  as  severe  as  they  were  for  the  Mach  1.9  cases.  Uke  die 
Mach  1.9  case,  increasing  the  nozzle  pressure  ratio  increased  the  magnitude  of  the  first 
sqiaration  bubble. 

The  nozzle  {xessure  ratio  effects  are  presented  in  terms  of  the  integrated 
pressure  contribution  to  thrust  in  Figure  4-23.  The  overall  magnitude  of  the  integrated 
thrust  tends  to  increase  linearly  with  increases  in  nozzle  pressure  ratio. 

For  a  given  nozzle  pressure  ratio,  the  magnitude  of  the  forces  are  much  greater 
for  the  Mach  1.9  case  than  fra-  the  Mach  3.0  case.  This  occurs  because  the  dominant 
cmitribution  to  nozzle  wall  pressure  thrust  is  the  high  pressure  acting  on  the  internal 
SO  degree  compression  ramp.  The  internal  nozzle  pressure  was  significantly  higher  for 
the  Mach  1.9  cases  than  it  was  for  the  Mach  3.0  cases  and  dierefore,  the  magnitude  of 
the  overall  nozzle  wall  pressure  thrust  is  much  larger  for  the  Mach  1.9  case. 
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The  magnitude  of  the  gaa  generator  contribution  to  thruat  and  lift  increases 
with  nosde  pressure  ratio  due  to  the  increase  in  the  size  of  the  first  8q>aration  bubble 
with  increasing  nozzle  {Messure  ratio.  Figures  4-24  presents  the  nozzle  pressure  ratio 
effects  in  terms  of  the  percent  improvement  in  nozzle  wall  pressure  thrust  develqjed 
by  (q)erating  dte  gas  generator  at  100%  throttle.  Although  the  magnitude  of  the 
pressure  contributions  to  thrust  increased  with  increasing  nozzle  pressure  ratio,  the 
percent  improvement  decreased.  The  much  larger  overall  increase  in  magnitude 
associated  with  increasing  the  nozzle  pressure  ratio  overwhelmed  the  much  smaller 
improvements  associated  with  the  gas  generator  flow.  For  these  reasons,  the  net 
effects  of  the  gas  generatw  on  the  nozzle  wall  pressure  recovery  are  much  more 
significant  at  the  lower  nozzle  pressure  ratio,  off-design  operating  conditions. 
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Figure  4-6  Gas  GenerateH’  Effects,  Mach  3.0,  NPR=16.0 
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Figure  4*7  Constant  Axial  Location  Profile,  Mach  1.9,  NPR=7.0, 1^52 
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Figure  4-8  Constant  Axial  Location  Flow  Profile,  Mach  3.0,  NPR=16.0, 1=52 


Figure  4-9  Gas  Generabx'  Mass  Flow  Effects,  Madi  1.9,  NFR=7.0  Nozzle  Wall  Pressure 
Distiibutioii,  No  Gas  Generator  DeflecHon 


Thrust,  No  Gas  GenerattH’  Deflection 
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Figure  4-11  Gas  Generator  Mass  Flow  Efifects,  Mach  3.0,  NI^si6.0  Nozzle  Wall 
Pressure  Distributirm,  No  Gas  Generatw  Deflection 


Gas  Genoator  Mass  Flow  Effects 


Figure  4-12  Gas  Generate  Mass  Flow  Effects,  Mach  3.0,  NPR=16.0  Nozzle  Wall 
Pressure  Thrust,  No  Gas  Generator  Deflection 
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Gas  Generator  Deflection  Effects 


Figure  4-13  Gas  Generator  Deflectiixi  ^ects,  Mach  1.9,  NPR=7.0  Nozzle  Wall  Pressure 
Distribution 


Figure  4-14  Gas  Generator  Deflection  Effects,  Mach  1.9,  NPR=7.0  Nozzle  Wall  Pressure 
Thrust 
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Gas  Geoerator  Deflecdoa  Effects 
Mach  3.0.  NPRsl6.0 
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Figure  4-15  Gas  Generator  Deflection  Effects,  Madi  3.0,  NPRsl6.0  Nozde  Wal 
I^essure  Distribution 


Figure  4-16  Gas  Generator  Deflection  ^ects,  Mach  3.0,  NPRs:16.0  Nozzle  Wall 
Pressure  Thrust 
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Figure  4-17  Experimental  Nozzle/Cowl  Oemnetry  Efifects,  1.9,  NFRb7.0 
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Figure  4-18  Ejqjerimental  Nozzle/Cowl  Geometry  Effects,  Mach  3.0,  NPRs:16.0 
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Figure  4-20  Nozzle  Pressure  Ratio  Elfects,  Mach  1.9,  Gas  Generator  On  100% 


Figure  4-21  Nozzle  Pressure  Ratio  Effects,  Mocb  3.0,  Gas  Generator  Off 
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Figure  4-22  Nozzle  Pressure  Ratio  Effects,  Mach  3.0,  Gas  Generator  On  100% 


Figure  4-23  Nozzle  Pressure  Ratio  Effects,  Nozzle  Wall  Pressure  Thrust 


Figure  4-24  Nozzle  Pressure  Ratio  Effects,  Percent  Change  in  Thrust 
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43  Application  of  a  Cowl  Gas  Generator  Flow  to  a  Generic  Hypersonic  Vehicle 
Configuration  Over  a  Typical  On-Design  Trajectory 

The  generic  hypersonic  nozzle/cowl  geometry  fxesented  in  Figure  1-4  was 
evaluated  with  and  without  the  gas  generator  over  a  typical  trajectmy.  The  gas 
generator  mass  flow  and  deflection  angle  effects  were  analyzed  at  four  design 
trajectory  points;  Mach  10,  15,  20,  and  25. 

The  baseline  nozzle  wall  pressure  distributions  fn*  all  four  trajectory  points 
with  the  gas  generatiff  off  are  i»esented  in  Hgure  4-25.  For  each  of  these  cases,  the 
flow  is  initially  very  underexpanded.  As  the  flow  exits  the  combusUx-,  it  undergoes  a 
dramatic  expansion  to  nearly  ambient  cmiditions.  As  the  trajectory  Mach  number 
increases,  the  initial  expansion  becomes  more  severe.  For  the  Mach  20  and  25 
trajectory  points,  the  flow  initially  owexpands  past  the  ambient  pressure.  All  of  the 
flows  are  characterized  by  a  small  but  steq>  reconq>ression  immediately  following  the 
initial  e}q)ansion  due  to  a  very  small  initial  flow  sq>aration  just  downstream  of  the 
throat.  The  nq;>id  increase  in  area  caused  by  dicular  arc  transition  from  the  combustw 
exit  to  the  nozzle  wall  creates  an  initial  sq)arati(m  re^on  followed  by  an  almost 
immediate  reattachment  of  the  flow  just  downstream  of  the  combustor  exit.  This  flow 
sqMtration  has  a  negligible  effect  on  the  omall  nozzle  performance.  The  majority  of 
the  nozzle  wall  pressure  distribution  is  dmninated  by  the  nozzle  wall  curvature. 

At  the  Mach  10  trajectory  point,  a  very  significant  recompression  takes  place 
alrmg  the  initial  curvature  of  the  nozzle  wall.  As  the  trajeeUxy  Mach  number 
increases,  the  magnitude  of  the  recompression  associated  with  the  nozzle  wall 
curvature  decreases  but  the  length  over  which  the  recompression  occurs  ino-eases. 
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The  locatiao  of  the  peak  recompression  also  moves  downstream  as  the  trajectory  Mach 
is  increased.  This  occurs  because  the  combustor  exit  pressure  is  relatively  high  and 
the  combustor  exit  Mach  number  is  relatively  low  for  the  lower  speed  trajectory 
points.  As  tbe  trajectory  Mach  number  increases,  the  combustor  exit  Mach  number 
increases  and  the  combustor  exit  pressure  decreases  significantly.  Thereftxe,  the  flow 
entering  this  fixed  geometry  nozzle  is  already  significantly  expanded  for  the  high 
Mach  number  trajectory  points.  For  these  reasons,  a  strong  recompressitm  associated 
with  the  high  pressure  exit  conditions  is  not  evident  fw  the  high  Mach  number 
trajecuxy  points. 

For  the  Mach  25  case  the  magnitude  of  the  recompression  is  very  small.  The 
Mach  25  trajectory  flow  exiting  the  combustor  initially  overe:q>ands  slightly  and  then 
immediately  recompresses  to  nearly  ambient  pressure  just  downstream  of  die 
combustor  exit  A  relatively  weak  recompression  does  eventually  occur  and  it  causes 
an  increase  in  |xxssure  along  a  fairly  large  portion  of  the  nozzle  wall. 

For  the  Mach  10  trajectory  point  die  re-expansion  past  the  initial 
recompression  is  fairly  steep  and  the  flow  eventually  ends  up  being  slighdy 
overmqianded  at  die  end  of  the  nozzle  wall.  Although  the  magiitude  of  the 
recompression  decreases  as  the  trajectmy  Mach  number  increases,  the  delayed  re- 
mansion  causes  a  slight  increase  in  pressure  along  the  downstream  portion  of  the 
nozzle  wall.  The  flow  even  becomes  underexpanded  at  the  end  of  die  nozzle  wall  for 
the  higher  Mach  number  trajectory  points. 

Appendix  C  presents  the  pressure  and  Mach  number  contours  for  this  nozzle 
operating  over  the  four  design  trajectory  points  with  the  gas  generator  (grating  at 
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100%  throttle  with  no  deflection  and  operating  with  the  gas  ^nerator  off.  Figures  C-9 
and  C-10  present  the  pressure  and  Mach  number  cmitours  for  the  Mach  10  trajectory 
point  with  the  gas  generator  off.  The  dominance  of  dw  internal  nozzle  flow  by  the 
tremendous  initial  expansion  can  be  seen  as  evkienced  by  the  dramatic  decrease  in  the 
flowfield  pressure  just  downstream  of  the  combustor  exit  In  addition,  the  external 
plume  shocks,  the  plume  flow  boundaries,  and  the  internal  shocks  can  be  seen. 

Figures  C-11  and  C-12  show  the  same  case  with  the  gas  generator  operating  at  100% 
throttle  with  no  deflection.  The  slight  interactimi  of  the  flow  structures  with  the  gas 
generator  flow  can  now  be  seen  as  evidenced  by  the  eliminati(Mi  of  the  internal  flow 
shock  attadwd  to  the  end  of  the  cowl  when  the  gas  generate  is  on.  With  the  gas 
generator  operating  at  100%,  die  internal  nozzle  pressure  distribution  decreases  slightly 
due  to  the  elimination  of  this  shock.  The  nozzle  wall  pressure  does  not  change  in  the 
upstream  region  near  the  initial  recompression  and  the  nozzle  wall  pressure  actually 
decreases  slightly  at  die  downstream  locations. 

With  the  gas  generator  off,  a  very  strong  {xessure  rise  occurs  in  the  region 
bdiind  the  blunt  trailing  edge  of  the  cowl.  The  flow  initially  undergoes  a  dramatic 
eaqiansion  as  it  turns  around  the  sharp  comers  of  the  blunt  cowl.  As  the  internal  and 
external  flows  meet,  diey  are  then  turned  back  downstream  through  a  recompression 
shock.  The  rec<»npression  shock  is  evidenced  by  a  region  of  increased  pressure  just 
downstream  of  the  cowl  in  Figures  C-9,  C-13,  C-17,  and  C-21  for  the  Mach  10,  IS, 

20,  and  25  trajectory  points  respectively.  The  magnitude  of  this  recompiession 
increases  with  increased  trajectory  Mach  number.  This  relatively  high  pressure  affects 
both  sides  of  the  flow  just  downstream  of  the  cowl.  For  the  external  flow,  this 
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inoease  in  pressure  results  in  a  strengthening  of  the  external  plume  shock  attached  to 
the  external  cowl  wall.  For  the  internal  nozzle  flow,  this  pressure  rise  propagates 
downstream  and  eventually  results  in  a  sli^dy  increased  nozzle  wall  pressure. 

With  the  gas  generator  operating  at  100%,  the  slight  changes  in  the  flow  fxofile 
behind  the  cowl  blunt  trailing  edge  can  be  seen.  For  these  cases,  the  gas  generator 
exhaust  pressure  at  the  end  of  the  cowl  is  much  lower  than  the  fxessure  rise  associated 
with  tl»  r^ompression  shock  fcx  die  gas  generator  off  cases.  The  decrease  in 
pressure  results  in  a  weaker  extonal  flow  plume  shock  and  an  eliminati(Mi  of  die 
internal  flow  shock  near  the  end  of  the  cowl.  This  is  evidenced  by  a  relative  decrease 
in  pressure  just  downstream  of  the  cowl  in  Figure  C-11,  C-IS,  C-19,  and  C-23  for  the 
Mach  10,  IS,  20,  and  25  trajectory  points  respectively.  This  re’rdve  decrease  in 
internal  pressure  shows  as  a  sUght  decrease  in  the  downstream  nozzle  wall 
pressure. 

An  af^xedation  of  the  magnitude  of  these  gas  generator  effects  can  be  gained 
by  investigating  die  flow  (xofile  just  downstream  of  the  gas  generator.  Figures  4-26 
through  4-29  (xesent  the  pressure  distiibutimi  at  a  cmistant  axial  location,  indicated  by 
index  "I”,  just  downstream  of  the  cowl  blunt  edge  for  the  Mach  10,  IS,  20,  and  2S 
trajectory  pdnts  respectively.  A  constant  1  location  of  S2,  correqxmding  to  an 
qiproximatB  downstream  location  of  10.19  inches,  was  used  to  generate  these  profiles. 
Note  that  the  cowl  edge  falls  between  hsl.O  inch  and  h=1.2S  inches.  The  internal 
nozzle  flow  is  the  region  h<1.0  inch  and  die  external  plume^peestream  flow  is  the 
region  h>1.2S  inches.  For  reference,  note  that  the  nozzle  wall  is  located  at  h=-3.903 
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inchei.  For  all  of  these  flow  cases,  the  internal  nozzle  flow  has  already  undergone  a 
significant  initial  expansion. 

For  the  higher  Mach  20  aixl  25  trajectory  points,  the  internal  nozzle  flow  is 
highly  expanded  across  the  entire  nozzle  at  this  location.  For  the  lower  Mach  10  and 
IS  trajectory  points,  the  flow  is  still  highly  expanded  but  die  flow  altmg  the  nozzle 
wall  is  dominated  by  the  nozzle  wall  curvature  recompressimi.  The  internal  nozzle 
flow  along  the  upper  cowl  wall  remains  eiqianded  despite  die  large  pressure  rec-^ 
al(»](g  the  nozzle  wall. 

With  the  gas  generator  off,  the  region  just  downstream  of  the  cowl  wall 
experiences  a  sharp  rise  in  jxessure.  The  pressure  behind  the  blunt  wall  is  actually 
higher  than  die  gas  generator  exhaust  {xessure  for  all  four  trajectory  points.  As  would 
be  expected,  the  magnitude  of  this  pressure  rise  increases  significantly  with  die 
increase  in  operatiiig  Mach  number.  The  internal  flow  reconqxesston  shock  shows  iq> 
as  a  large  spike  near  the  inner  Iqi  of  tlw  cowl  at  qiproximately  1.0  inches.  The 
weaker  external  flow  recompressitm  shock  shows  iq>  as  a  smaller  spike  near  die  outer 
edge  of  die  cowl  Figures  4-26  through  4-29  show  that  the  magnitude  of  diese 
recooiqxessioos  increase  widi  increasing  trajectory  Mach  number. 

With  die  gas  generator  on,  the  internal  flow  eiqiansion  almig  the  cowl  wall 
results  in  a  jxessure  diat  is  almost  identical  to  the  gas  generator  exhaust  pressure  for 
the  Mach  10  trajectory  point.  As  the  cperating  Mach  is  inaeased,  the  internal  nozzle 
fvessure  along  the  cowl  wall  increases.  Per  the  higher  Mach  15,  20,  and  25  trajectory 
points,  the  gas  generator  exhaust  pressure  is  actually  lower  than  the  internal  nozzle 
pressure  and  the  internal  nozzle  flow  expands  out  into  this  region.  The  increased 
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expansion  firom  the  internal  nozzle  flow  to  the  external  nozzle  flow  shows  up  as  a 
slight  decrease  in  nozzle  wall  pressure  as  the  effects  propagate  downstream. 

43.1  Gas  Generator  Mass  Flow  Effects.  Figure  4-30  shows  the  effects  of  die 
gas  generator  on  the  nozzle  wail  pressure  recovery  for  the  Mach  10  trajectory  point. 
The  change  in  pressure  distribution  associated  with  the  gas  genoator  is  relatively 
small.  The  net  effects  are  limited  to  small  pressure  changes  at  the  downstream 
locations  with  virtually  no  change  in  pressure  alopg  die  first  25  percent  of  the  nozzle 
wall.  For  a  fixed  gas  genoator  geometry,  die  gas  generator  eidiaust  pressure  increases 
widi  increased  throttle  setting.  The  gas  generator  exit  pressure  is  high  enough  to  cause 
an  increase  in  downstream  nozzle  wall  pressure  relative  to  die  gas  generator  off  case 
only  when  die  gas  generator  is  operating  at  die  150%  and  200%  tibrotde  settings. 

Figure  4-31  shows  the  magnitude  of  integinated  pressure  contribution  to  thrust 
along  the  nozzle  walL  The  nozzle  perfocmance  increases  sli^dy  with  increasing  mass 
flow  due  to  the  downstream  {vopagatirm  d  die  higher  gas  generator  exhaust  pressure. 
The  overall  performance  is  only  better  dum  the  baseline  gas  generator  off  performance 
for  the  highest  throtde  settings. 

figures  4-32, 4-34,  and  4-36  show  the  gas  generator  mass  flow  effects  on 
nozzle  wall  pressure  distribution  for  die  higher  Mach  number  trajectory  points. 
Similarly,  Rgures  4-33, 4-35,  and  4-37  diow  the  effects  of  gas  ^nerator  mass  flow 
changes  <m  the  integrated  nozzle  wall  pressure  contribution  to  thrust  for  die  higher 
Mach  15, 20,  and  25  trajectory  points  respectively.  The  overall  trends  associated  with 
dm  gas  generate  are  the  same  for  each  case;  as  the  gas  generator  mass  flow  is 
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increased,  die  nooszle  wall  pressure  recovery  and  subsequent  perfonnance  increases 
sli^dy.  Since  the  nozzle  internal  pressure  along  the  cowl  wall  increases  with 
trajectory  K&ch  number,  and  die  gas  generator  exit  jxnssure  is  fixed  for  a  given 
throttle  setting,  the  gas  generator  contribution  to  nozzle  wall  performance  decreases 
with  increasing  trajectory  Mach  number. 

432  Gas  Generator  Deflection  Effects.  Figure  4-38  shows  the  effects 
of  deflecting  the  gas  generator  for  the  Mach  10  trajectory  point.  Deflecting  the  gas 
generator  towards  die  nozzle  wall  cauaes  a  fairly  significant  increase  in  fxessure, 
especially  at  die  very  large  deflection  angles.  A  20  degree  deflection  towards  the  wall 
actually  causes  a  small  second  separation  bubble  to  form.  This  is  caused  by  the 
impingement  of  the  shock  generated  by  die  interaction  of  the  internal  nozzle  flow  widi 
the  sharp  inward  deflection  of  the  gas  generator  flow.  The  effect  of  ddlectug  die  gas 
generator  away  from  the  wall  is  very  small  and  die  net  effect  moves  firrdier 
downstream. 

Kgure  4-39  shows  the  magnitude  ci  int^rated  pressure  contribution  to  thrust 
along  the  nozzle  wall.  The  nozade  performance  increases  significandy  with  deflections 
towards  the  wall  due  to  die  development  of  the  small  downstream  8q;>aration  bubble. 
As  the  gas  generator  is  deflected  away  from  die  wall,  only  a  small  additional  penalty 
is  paid  for  the  increased  erqiansion  near  die  end  of  the  cowl. 

Figures  4-40,  4-42,  and  4-44  show  the  gas  generator  deflection  effects  <m 
nozzle  wall  pressure  distribution  for  die  higtor  Madi  number  trajectory  points. 
Similarly,  Figures  4-41,  4-43,  and  4-45  show  the  effect  of  gas  generator  deflectirms  on 
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the  interned  nozzle  wall  pfeesure  cootributioa  to  thrust  for  the  higher  Mach  IS.  20. 
and  2S  trajectory  points  respectively.  The  overall  trends  associated  with  the  gas 
generator  deflecdoos  are  die  similar  for  ewdi  case.  For  die  same  reasons  discussed 
above,  the  overall  gas  generator  contributimi  to  nozzle  wall  performance  decreases 
widi  increasing  trajectory  Mach  number. 

4JJ  Tnyedory  Perfonnance  Analysis.  The  trends  discussed  in  previous 
sections  need  to  be  put  into  perqiective.  This  section  presents  the  nozzle  wall  fvessure 
contributiois  to  dirust  and  lift  relative  to  the  net  forces  acting  on  the  generic 
nozzk/cowl  afterbody.  The  nozzle  wall  fuessure  contribution  to  thrust  and  lift  were 
calculated  usiag  Equations  2-44  and  2-45. 

Kgure  4-46  presents  die  gas  generator  mass  flow  effects  and  Figure  4-47 
presents  the  gas  generator  deflection  angle  effects  in  terms  of  the  pocent  change  in 
the  pressure  contribution  to  the  nozzle  wall  lift.  Except  for  the  high  mass  flow 
conditioos  and  die  large  inward  deflection  angle  cases,  the  gas  generator  results  in  a 
fairly  significant  decrease  in  the  ovmull  nozade  wall  pressure  cmitribution  to  lift.  The 
overall  trends  are  die  same  as  discussed  in  die  previous  two  sections. 

The  gas  generator  also  has  a  direct  force  contribution  to  lift  when  it  is 
deflected.  This  contribution  must  be  accounted  for  to  assess  die  net  effect  of  gas 
generator  d^ecticms  <n  the  complete  nozzle/cowl  afterbody  lift.  Figure  4-48  shows 
the  perceot  change  in  net  lift  due  to  b<^  the  nozzle  wall  js'essure  contributimi  and  the 
direct  gas  ^nerator  deflection  contributirm  to  lift.  Although  gas  geimrator  deflectirms 
towards  the  wall  increase  the  nozzle  wall  i»%ssure  contribution  to  lift,  the  deflectimis 
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towards  die  wall  ioipait  a  direct  thrust  coatribudoa  opposing  lift.  The  decrease  in  lift 
due  to  the  deflected  dinist  cono^wnent  more  than  cancel  out  the  increased  nozzle  wall 
pressure  lift  The  net  lift  forces  inaease  dramatically  with  deflectimis  away  from  the 
wall  because  the  direct  lift  term  associated  with  the  gas  generator  deflection  is  not 
of^KMod  by  a  very  significant  decrease  in  nozzle  wall  pressure  lift 

Figure  4-49  presents  the  gas  generator  mass  flow  effects  with  no  gas  generator 
deflecti<»  in  terms  of  the  percent  changes  that  the  nozzle  wall  jxessure  crmtributioiis 
have  on  the  net  nozzWcowl  afterbody  thrust  The  overall  magnitude  of  the  ^ 
generator  effects  on  the  nozzle  wall  pressure  thrust  are  very  small  compared  to  the 
total  durust  For  die  worst  case  Mach  25  trajectory  point  the  net  decrease  in  thrust  is 
less  than  0.19  percent  of  the  total  thrust 

Rguie  4-50  presents  die  gas  generator  nutts  flow  effects  with  no  gas  generator 
deflection  in  terms  of  percent  increase  in  net  dmist  due  to  bodi  the  nozzle  wall 
laessure  contribution  and  the  direct  gas  generator  contributions  to  thrust  The  effect  of 
die  gas  generator  direct  thrust  contributicm  to  total  thrust  is  nearly  2  orders  of 
magnitude  greater  than  die  nozzle  wall  jxessure  change  contributi(»s  to  total  dirust 
Rgiffe  4-51  presents  die  gas  generator  deflection  effects  widi  the  gas  generator 
operating  at  100%  throtde  in  terms  of  the  percent  changes  duU  the  nozzle  wall 
(xessure  contributions  have  on  the  net  nozzle/cowl  afterbody  thrust  Again,  for  die 
best  case  Mach  10  trajectory  point,  the  net  increase  in  thrust  is  less  dian  0.26  percent 
of  the  tmal  dirust 

Kgure  4-52  presents  die  gas  generator  deflection  effects  with  the  gas  generator 
operating  at  100%  throtde  in  terms  of  percent  increase  in  net  thrust  due  to  both  the 
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nooszle  wall  pressure  contributicm  sad  the  direct  gta  generator  contributions  to  dirust 
Again,  die  effe^  of  the  gas  generator  direct  thrust  c(»itributi<»  significandy  dominates 
the  nozzle  wall  {»essuie  change  contribudoos  to  total  thrust  Obviously,  the  direct  gas 
generator  cmitributicMi  to  thrust  decreases  widi  deflections  in  either  direction.  The 
direct  gas  generator  thrust  more  than  cmnpensates  for  any  relative  losses  or  gains 
associated  with  the  pressure  changes  on  the  nozzle  wall. 
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Constant  Axial  Location  Flow  Profile 
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Constant  Axial  Location  Flow  Profile 
Mach  20A  1^2 
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Gmstant  Axial  Location  Flow  Profile 
Mach  2S.0, 1=62 
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Figure  4-29  Constant  Axial  Location  Flow  Profile,  Mach  25.0, 1~52 


Figure  4-30  Gas  Generator  Mass  Plow  Effects,  Mach  10  Nozzle  Wall  Pressure 
Distributicni.  No  Gas  Generator  Deflection 


Figure  4-31  Gas  Generatcnr  Mass  Flow  Effects,  Mach  10  Nozzle  Wall  Pressure  Thrust, 
No  Gas  Generator  Deflection 
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Oas  Oenorator  Mass  Flow  Effects,  Mach  IS 


Figure  4-33  Gas  Generate  Mass  Flow  Effects,  Madi  IS  Nozzle  Wall  Pressuie  Thrust, 


No  Gas  Generator  Deflection 
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F^ure  4>34  Gas  Generator  Maas  Flow  Effects,  Mach  20  Nozzle  Wall  Pressure 
Distributiofi,  No  Gas  Generator  Deflectitm 


Figure  4-35  Gas  Generator  Mass  Flow  Ejects,  Madi  20  Nozzle  Wall  Pressure  Thrust, 
No  Gas  Gentfatw  Deflection 
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Oas  Oenerator  Mass  Flow  Effects,  Madi  25 


Figure  4*36  Gas  Generator  Mass  Plow  Effects,  Mach  25  Nozzle  Wall  Pkessure 
Distiibutioa,  No  Gas  Generator  Deflection 


Figure  4-37  Gas  Generator  Mass  Flow  Effects,  Madi  25  Nozzle  Wall  Pressure  Thrust, 
No  Gas  G«ierat(»r  Deflection 


Fil^ire  4-38  G«a  Generator  Deflecti(»  Effects,  Mach  10  Nozzle  Wall  Pleasure 
Distributioii 
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Ots  Qenerator  Deflection  Effects,  Mach  15 


Figure  4-40  Gas  Geneistor  Deflectioa  ^fects,  Mach  15  Noazle  Wall  Pressure 
Distributioii 


Figure  4-41  Gas  Generator  Deflection  Effects,  Mach  IS  Nozzle  Wall  Pressure  Thrust 


Gas  Generator  Deflection  Bfiects,  Mach  20 


Figure  4-42  Gas  Generator  Deflection  Effects,  Mach  20  Nozzle  Wall  Pteasure 
Distribiitioo 


Figure  4-43  Gas  Generator  Deflection  Effects,  Mach  20  Nozzle  Wall  Pressure  Thrust 
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Figure  4^  Ges  Oenerattr  Deflection  Effects,  Madi  25  Nozzle  Wall  Plressure 
Distributioa 


Figure  4>4S  Gas  Generator  Deflection  Effects,  Mach  25  Nozzle  Wall  Pressure  Thrust 
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Figure  4^6  Gas  Generator  Mass  Flow  Effects,  Percent  Obange  in  Nozzle  Wall  Lift,  No  Gas  Generator  Deflection 
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Figure  4^7  Gas  Generator  Deflection  Effects,  Nozzle  Wall  Riessure  Contribution  to  Lift,  Gas  Generator  Throttle  100% 


Figure  4-48  Gas  Generator  Deflection  Effects,  Nozzle  Wall  Pressure  and  Direct  Gas  Generator  Contribution  to  Lift.  Gas  Generiuor 
Throttle  100% 
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Figure  4-51  Gas  Generator  Deflection  j^scts.  Nozzle  Wall  Rressure  Contribution  to  Net  Thrust,  Gas  Generator  Throttle  100% 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


5.1  Suminary 

The  proposed  improved  cowl  concept  was  evaluated  for  two  hypersonic  vehicle 
nozzle/cowl  afterbody  configurations.  The  effect  of  inserting  a  gas  generator  into  the 
cowl  trailing  edge  was  investigated  numerically  using  an  e]q)erimentally  validated 
finite  voltime  computati<mal  code  developed  by  Wright  Laboratory.  The  two 
dimensional  Navier-Stokes  equations  wo'e  solved  assuming  laminar  flow  and  a  perfect 
gas  equation  of  state  using  two  flux-split  methods.  The  flowfields  were  solved  using  a 
cmnbinatitm  of  the  e]q)licit  formulations  of  the  Van  Leer  flux-vector  q>litting  and  the 
Roe  flux-difference  splittiqg  algorithms.  Hie  effects  of  gas  generator  man  flow 
increases  and  thrust  vector  deflecticm  angles  were  evaluated  in  terms  of  the  pressure 
distribution  on  tiie  nozzle  wall  and  its  effect  on  overall  nozzle  performance. 

The  experimentally  evaluated  hypersonic  nozzle/cowl  configuration  was 
analyzed  over  a  raoge  of  nozzle  pressure  ratios  for  two  supersonic  <rff-design  Mach 
numbers.  The  flowfields  were  dmninated  by  the  initial  oqiansitm  of  the  nozzle  flew. 
The  extrat  of  die  fxenure  recovoy  almig  tiie  nozzle  wall  was  affected  by  tiie  size  and 
location  of  the  sqiaration  re^ons.  Mach  number,  nozzle  pressure  ratio,  and  gas 
generator  throttle  setting  all  had  a  significant  impact  on  the  locaticm  and  size  of  the 
sq>arati<»i  regimis. 

The  trends  associated  with  inserting  a  gas  generator  into  the  generic 
nozzWcowl  afterbody  configuration  on  the  high  Mach,  on-design  trajecbxy  were  not 
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as  aignifirant  as  tiiey  wcrc  for  the  experimentally  validated  off-deaign  ctxifiguratioa. 
For  the  generic  nozzle  ctmfiguration.  the  flow  was  also  dominated  by  the  very 
dramatic  initial  mq>an8ioa.  Gas  generator  effects  on  tlw  nozzle  wall  pressure  recovery 
were  limited  to  small  changes  in  nozzle  wall  luessuie  at  the  downstream  locations. 

For  a  majority  of  the  operating  envelope,  the  gas  generator  actually  caused  a  slight 
decrease  in  nozzle  wall  {Hessure  recovery. 

5J2  Condusions 

For  die  experimental  off-design  configuration,  inserting  die  gas  gmerator  into 
the  end  of  the  cowl  increased  die  pressure  throu^  the  first  sqiaration  bubble  resulting 
in  a  higher  pressure  recovery  along  the  nozde.  The  hi^  {xessure  boundary  created  by 
die  gas  generator  exhaust  delayed  the  internal  nozzle  flow  expansion  resulting  in  a 
higher  overall  nozzle  wall  fxessure  distributioii.  Furthermore,  die  jxesenoe  of  the  gas 
genmator  tended  to  minimize  die  internal  nozzle  shock  attached  to  the  cowl  which 
eliminated  the  second  sqiaratkm  regitm  associated  widi  the  low  qieed,  low  nozzle 
{xessure  ratio  off-design  flow  conditions. 

For  the  eiqierimental  off-design  configuration,  the  gas  goimtor  throtde  setting 
had  a  significant  inqiact  on  die  overall  nozzle  wall  {xessure  recovery.  The  primary 
variable  affecting  the  nozzle  internal  flow  was  the  high  pressure  boundary  created  by 
the  gas  generator.  As  die  gas  graerator  throtde  setting  increased,  both  the  gas 
generator  mass  flow  and  exhaust  pressure  increased.  The  higher  gas  genoator  exhaust 
{xcssure  resulted  in  a  significantly  higher  nozzle  wall  pressure  recovery  near  the  first 
separation  bubble. 
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Par  Hkt  ejqperimental  coofiguntion,  the  gu  generator  deflectioa  effects  were 
also  significant.  Defiectioas  towards  the  wall  dramatically  increased  die  nozzle  wall 
inessure  recovery  while  deflections  away  from  the  wall  only  caused  relatively  small 
losses.  Most  of  the  effects  of  gas  generator  deflections  were  isolated  to  the  tqistream 
locations  near  the  first  separaticm  bubble.  Gas  generator  deflections  had  very  little 
impact  on  the  downstream  nozzle  wall  {xessure  recovery. 

The  comparismi  of  gas  generator  effects  to  nozzle  cowl  geometry  chaitge 
effects  was  ]»romisiog.  The  gas  generator  effectively  acts  as  an  aerodynamic  cowl 
extensicm.  The  results  associated  with  the  additicm  of  the  gas  generator  were  very 
similar  to  tiie  effects  of  increasing  the  length  of  the  cowl.  Deflections  of  die  cowl 
extension  effects  also  follow  die  same  trends  associated  with  the  gas  generator 
deflections.  An  int^irated  cowl  gas  generator  has  the  potential  of  saving  overall 
wei^t  due  to  a  decrease  in  reqiiired  cowl  lengdi. 

For  the  genetic  nozzle  cowl  coofiguratitm  analyzed  over  a  typical  trajectory, 
the  effects  of  the  gas  generator  on  nozzle  wall  pressure  recovery  were  not  very 
significant  in  terms  of  overall  nozzle  performance.  The  small  changes  in  nozzle  wall 
pressure  recovery  were  primarily  a  result  cS  die  propagation  of  pressure  differences 
behind  the  blunt  cowl  to  die  nozzJe  wall  downstream  locations.  When  the  gas 
gmierator  was  turned  off,  die  recompression  shock  behind  the  blunt  trailing  edge  of  dm 
cowl  caused  a  very  large  pressure  rise  in  this  r^tm.  Although  diis  fxessure  rise 
associated  with  the  blunt  cowl  increased  the  downstream  nozzle  wall  pressure  recovery 
sU^dy,  it  ^nm-ated  a  strong  external  flow  shock  resulting  in  increased  plume  drag. 

For  the  generic  nozzle/cowl  configuration,  the  overall  effects  of  the  gas 
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oo  wn«aii#t  wall  {xcssure  recovery  were  aecondary  to  the  net  effects  of  the  gas 
gpf^rat^  on  overall  nozzle  performance.  As  of^>oeed  to  the  nozzle  wall  pressure 
recovery  changes,  the  direct  thrust  and  lift  terms  associated  with  the  gas  generator 
monriftntiiin  flux  c<mtributed  significantly  to  the  overall  nozzle  performance. 

The  potential  beneff  ts  for  insertiiig  a  gas  generator  into  the  cowl  of  a 
hypersonic  vehicle  were  demonstrated  for  both  omfiguratioos  over  a  wide  range  of 
off-design  and  (m-design  operating  ccmditioos.  For  the  low  speed  off-design  cases,  die 
gas  gfmMratfw  matte  a  aignififjint  cmitribution  to  oversll  mizzle  wall  jnessure  recovery. 
As  proposed,  the  gas  generator  demonstrated  the  potential  to  act  as  an  aerodynamic 
cowl  extension. 

Although  die  gas  ynwotor  effects  on  nozde  waU  pressure  recovery  were 
minimal  for  the  generic  nozzle  operating  at  the  high  h^adi  number  on-design  trajectory 
points,  die  gas  generator  effects  on  the  overall  int^irated  nozzk/cowl/gas  generator 
poformance  were  significant.  In  particular,  die  direct  dirust  tom  associated  widi  the 
gas  gnnarafaa-  contributes  significandy  to  overall  performance.  Ihe  offset  location  of 
the  gas  generator  also  provides  the  potential  to  make  a  significant  contribution  to 
overall  vehicle  stability  and  cmitrol  moment  control.  Direct  moment  contrtd  could  be 
achieved  throu^  any  desired  combination  of  gas  generator  deflections  and  gas 
ymf-mror  miHM  flow  settings.  This  direct  moment  control  has  the  potential  to  reduce 
the  trmnoidous  trim  drag  penalties  associated  with  crmventional  crmtrol  surface 
deflections  typically  required  for  flight  vdiides  operating  at  hypersonic  speeds. 
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Rccommcmtotionn 


Purther  investigtdoo  of  the  improved  cowl  c(»oept  is  warrinted.  The  Iriggest 
limitatioa  of  die  piesent  investigsdoo  was  that  it  wss  limited  to  an  investigatioa  of 
two-dimensional  effects.  To  gain  a  thorou^  understanding  of  the  complex  flowfield, 
a  full  three  dimensional  model  diould  be  evaluated.  Ideally,  a  three  dimensional 
hypersooic  vehicle  afterbody  pmnt  design  wmild  be  evaluated  over  an  optimized 
trajectory  usiqg  a  realistic  cowl-mt^tvted  gas  generator. 

Short  of  a  full  diree-dimensional  analym,  further  refinements  to  the  two- 
dimensional  i^iproach  would  be  beneficial.  The  incorporation  of  turbulence  into  die 
Wright  Laboratory  code  would  he4>  ovmcome  some  of  the  errors  associated  widi  the 
exact  shape  of  die  first  sqparatioo  bubble.  Additional  modifications  that  wouhl  add 
realism  to  the  code  for  die  high  qieed,  hij^  temperature  flows  would  be  die 
incorporation  of  chemically  dissimilar  flowfields  widi  flowfield  muring  or  die  posadile 
incorporation  of  chemical  reactions  into  the  code. 

Additional  variations  on  die  improved  cowl  conoqit  should  also  be  explored. 
Since  the  gas  genoator  performance  is  highly  dependent  on  the  extent  of  the  intenial 
flow  eiqiansion  aloitg  cowl  wall,  die  on-design,  high  Mach  number  performance  of  a 
gas  goierator  concept  should  be  evaluated  using  more  realistic  nozzle  configurations 
and  trajectories. 

Finally,  a  trade  study  riiould  be  done  to  determine  die  benefits  associated  with 
increasing  die  gas  generator  ochaust  pressure  at  the  oqiense  of  a  decrease  in  exit 
momentum.  This  could  easily  be  accomplished  by  parametrically  decreasing  the  gas 
generator  nozzle  expansimi  ratio. 
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Appendix  A:  Code  Inputs 

This  q)peiidix  cootains  a  sample  iqMit  file  and  a  descr4>tion  of  die  ii^ts  for 


the  Wright  Laboratory  computation  fluid  dynamics  code  used  in  this  research.  The 
ii^Nits  can  be  broken  down  into  several  categories:  algorithm  control  parameters, 
geometric  parametm,  flowfield  initial  conditions,  and  basic  solutimi  contnd 
parameters.  The  most  significant  parameters  are  highlighted  in  the  following  tables. 


TaUe  A-1  Algoiidun  Control  Parameters 


ISWVL 

Solution  method;  Van  Leer  (3),  Roe  (4) 

NEND 

Number  of  iterations 

INS 

Flow  modd;  Euler  (0).  Navier-Stdees  (1) 

CPLMAX 

Maximum  allowable  CPL 

CPL 

Initial  CPL  used  to  start  solution 

CFLBXP 

Number  of  iterations  before  CPL  doubles 

ICFL 

Number  of  iterations  between  CPL  increases 

IMPLT 

Numerical  formulation;  Explicit  (0),  Inqilicit  (1) 

Table  A-2  Geometric  Parameters 


m  JL 

Grid  dimensions  in  ^  and  1\  reflectively 

INS1.INS2 

Grid  location  defining  cowl  edge  ^  location 

JNSl.  JNS2 

Grid  location  defining  cowl  edge  T)  location 

130 


Table  A-3  Plowfield  Gxiditioiii 


lADBWL 

Adiabatic  wall  (0)  or  constant  tm^wrature  wall  (1)  BC 

TWALL 

Wall  temperature 

ALPHA 

Pkeestream  angle  of  attack 

RM 

Freeslream  Mach  number 

REL 

Freestream  Reynolds  number  per  foot 

TINF 

Fireestream  static  temperature 

GASCONSTANTEXT 

Freestream  gas  constant  R 

GAMMAEXT 

Freestream  y 

AINALPHA 

Internal  flow  angle  of  attack 

RINM 

Internal  flow  Mach  number 

RINEL 

Internal  flow  Reynolds  number  per  foot 

TININF 

Riternal  flow  static  tonperatnre 

GASCONSTANTIN 

Intenial  flow  gas  constant  R 

GAMMAIN 

Internal  flow  y 

AGGALPHA 

Gas  generator  tiinist  vector  deflection 

RGGM 

Gas  generator  Mach  number 

RGGEL 

Gas  generator  Reynolds  number  per  foot 

TGGINF 

Gas  generator  static  temperature 

GASCONSTGG 

Gas  generator  gas  cmistant  R 

GAMMAGG 

Gas  generator  Y 

Table  A'4  Basic  C<mtrol  Paramelars 


IREAD 

Deadstart  (0)  or  Restart  (1) 

IGRID 

Grid  Pbrmat 

IP3DOP 

PlotSD  ou^ut  file  fwmat 

MODPR 

Iterations  between  [Minting  convergence  data 

The  following  sainple  input  file  rqxesents  the  Mech  15  flow  cooditioa  with  the 


gu  graentor  operating  at  100%  throttle  widi  no  deflectioo. 


cnldat 

IL,  JL,  IMETRC 
102  72  0 

DSHD 

2000 

ICFL,  CFLEXP,  CFLMAX,  CFL 
5  1  0.9  0.01 

IREST,  CFCRHO,  CFCEI,  CFLPEN,  CEXPPEM,  IMOFRZ 

1  3.00  3.00  1  1  1 

IMPLT,  ILCTST,  ISNVL,  lUtTR,  OHEGR,  DELTEP,  DELTIL,  nSMTB,  IMS,  NSMPS 
01  421  l.E-10  0.3  2  12 

lAOBHL,  ICASE,  TMALL,  ALPHA,  RMIHF,  REL,  RL,  TINF,  IGRID 

0  2  2000.00  0.  9.3347  7.3728E5  1.0  1113.21  2 

IREAO,  IP3DOP,  OGBOG,  MOOPR.  NRST,  IFMRTI,  IFMRTO,  IIMTl 
1  20  10  01  1  1 

IMSl,  IMS2,  JNSl,  JNS2 
1  32  37  46 

TIMNALL,  AIHLPHA.  RIIMIMF,  RIHEL,  RIML,  TIMIHF  IGRID 

2000.00  0.0  3.4443  2.8701E6  1.0  5511.96  2 

AggLPHA,  RggMIHF,  RggEL,  RggL,  TgglMF,  GASCONSTgg,  GAtMAgg 

0.  4.680  0.7S723E6  1.0  2158.29  3308.98  1.26059 

GASCONSTia  GAMMAin  GA5CONS7«xt  GAMMA«xt 

2237.34  1.25  1727.62  1.35 
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Appendix  B:  Code  Porallelixatioa 


This  appendix  contains  an  example  of  parallel  processing  techniques  applied  to 
the  solution  of  a  generic  nozzle/cowl  flowfield  at  the  Mach  IS  design  condition  with  a 
gas  generator  operating  at  100%  throttle  (Section  3.2).  For  this  example.  10(X)  Roe 
iteratimis  were  cmnpleted.  An  example  of  the  code  modifications  are  presented  for 
subroutine  L.f. 


B.1  Profile  of  Non  Parallelized  Code 

This  section  contains  the  execudmi  profile  of  the  original  code  without  any 
parallelizatitM.  This  ouq>ut  was  used  to  i(fentify  die  subroutines  that  consumed  the 
moat  CPU  time  and  were  therefore  the  best  candidates  for  parallelization.  Subroutines 
LMTRI,  LMTRJ,  FSIROE,  FSJROE,  VISEUL,  VJSEUL,  and  L  were  chosen  for 
parallelization  based  on  this  breakdown. 

Profil*  listing  gensrstsd  Fri  Jan  14  08:01:28  1994  with: 
prof  -pixis  -quit  1%  -clock  25  noxzlo.exs 

*  -pCrocadurss]  using  basic-block  counts;  • 

*  sorted  in  descending  order  by  the  number  of  cycles  executed  in  each  * 

*  procedure;  unexecuted  procedures  are  excluded  • 


130230201604  cycles  (3309.8083  seoende  at  25.00  megahertz) 


cycles  ' 

Icycles 

cum  % 

seconds 

cycles 

/call 

11962368802 

13.19 

13.19 

718.4948 

89812 

11844680000 

13.70 

27.50 

713.7872 

127462 

11511206120 

13.45 

40.95 

700.6883 

125123 

11186000000 

13.20 

54.14 

687.4400 

85930 

16364936000 

12.51 

66.11 

654.5974 

16364936 

10890460000 

8.36 

15.01 

435.6184 

77789 

10315555314 

1.91 

83.04 

415.0222 

42 

9522240000 

1.31 

90.35  : 

380.8896 

21 

8616600000 

6.62 

96.97  : 

344.6640 

43083 

1116125600 

1.36 

98.33 

71.0450 

1776126 

bytes  procedure  (file) 
/line 

93  Imtr j  (LMTRJ. f ) 

41  feiroe_  (FSIROE. f) 
90  Imtri  (LMTRI.f) 

40  fsjroe_  (FSJROE. f) 
56  1_  (L.f) 

80  vieeul_  (VISEOL.f) 
6  sqrtf  (sqrtf.s) 

23  r_sign  (r_sign.c) 

62  vjseul  (VJSEUL.  f) 
51  08tdt_~(ESTDT.f) 


*  -p ( rocedures ]  using  invocation  counts;  * 

*  sorted  in  descending  order  by  number  of  calls  per  procedure;  * 

*  unexecuted  procedures  are  excluded  * 


181998451  invocations  total 

calls  %calls  cumt  bytes  procedure  (file) 
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453440000  57.54  57.54 

247037108  31.35  88.89 

84420000  10.71  99.61 


112  r_al4n  (r_sign.c) 
248  aqrtf  (aqrtf.a) 

64  cvmqp_  (CVHGP.f) 


*  -h(aavy]  uaing  baaic-block  counta; 

*  aortad  in  dascending  ordar  by  tha  numbar 

*  llna;  unexacuted  lines  are  excluded 

of  cycles 

executed  in 

each 

* 

* 

procedure  (file) 

line 

bytea 

cycles 

% 

cum  % 

r  sign  (r  sign.c) 

11 

52 

4987840000 

3.83 

3.8.3 

r  sign  (r  sign.c) 

10 

44 

3627520000 

2.79 

6.62 

1  (L.f) 

83 

232 

3248000000 

2.49 

9.11 

lintrj  (LMTRJ.f) 

106 

344 

2435194846 

1.87 

10.98 

Imtri  (LMTRI.f) 

105 

344 

2425980342 

1.86 

12.84 

1  <l..f) 

39 

168 

2352000000 

1.81 

14.65 

Imtrj  (LMTRJ.f) 

91 

284 

1995181060 

1.53 

16.18 

Imtrj  (LMTRJ.f) 

98 

272 

1907973495 

1.47 

17.65 

Iffltrj  (LMTRJ.f) 

94 

272 

1902484691 

1.46 

19.11 

Imtri  (LMTRI.f) 

88 

232 

1618924982 

1.24 

20.35 

Imtri  (LMTRI.f) 

96 

208 

1449538271 

1.11 

21.46 

Imtri  (LMTRI.f) 

92 

208 

1442205804 

1.11 

22.57 
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Ba  ProfUe  of  Paralkliied  Code 

This  sectimi  contains  the  executicm  profile  for  die  parallelized  code.  CPU 
usage  is  brdcNi  down  for  the  four  {vocessors  used  on  die  Silicon  Graphics  4D/240 
machine. 


Processor  number  1: 

Profila  listing  genarated  Mon  Fab  14  14:16:04  1994  with: 

prof  -pixia  -quit  1%  -clock  2S  nozzla.axa  nozzla.exa.Addrs  nozzle. axa.Countsl2635 


*  -plrocedurasl  using  basic-block  counts; 

*  sorted  in  descending  order  by  the  number  of  cycles  executed  in  each 

*  procedure;  unexecuted  procedures  are  excluded 


50776235875 

cycles 

(2031.0494  mmooadm 

at  25.00 

megahertz) 

cycles 

kcycles 

cum  % 

seconds 

cycles 

bytes 

procedure  (file) 

/call 

/line 

9191735840 

18.10 

18.10 

367.6694 

45959 

2812 

_lmt  r  3_15 l_aaad_ 

(LMTRJ.f) 

8835036657 

17.40 

35.50 

353.4015 

63108 

406 

_lmt  r i~14  8_aaad_ 

(LMTRI.f) 

4761120000 

9.38 

44.88 

190.4448 

21 

23 

r_eign  (r_eign.c) 

3336355374 

6.57 

51.45 

133.4542 

42 

6 

sqrtf  (sqrtf.s) 

2955960000 

5.82 

57.27 

118.2384 

21114 

3392 

_fsiroe_90_aaac_ 

(FSIROE.f) 

2920260000 

5.75 

63.02 

116.8104 

803 

1620 

~viseul~5^Z****Z 

(VrSEDL.f) 

2867600000 

5.65 

68.67 

114.7040 

14338 

3332 

Zf • j roeZs  9ZaaacZ 

(FSJROE.f) 

2268000000 

4.47 

73.14 

90.7200 

630 

373 

~v3seul~57”’aaaa~ 

(VJSEOL.f) 

1776125600 

3.50 

76.63 

71.0450 

1776126 

57 

estdt  TESTDT.fl 

1433600000 

2.82 

79.46 

57.3440 

28672 

1804 

1  227  aaae  (L.f) 

1412880000 

2.78 

82.24 

56.5152 

5046 

912 

_f si roe_64_aaab_ 

(FSIROE.f) 

1335021836 

2.63 

84.87 

53.4009 

1299 

24 

mp  master  wait  for 

interleaved  loop  completion 

(mp.c) 

1321600000 

2.60 

87.47 

52.8640 

3304 

868 

_f 8  j  roe_63_aaab_ 

(FSJROE.f) 

1162771203 

2.29 

89.76 

46.5108 

855 

16 

~  mp  simple  sched 

(mp  simple. s) 

1109400000 

2.18 

91.95 

44.3760 

5547 

1776 

1  149  aaab  (L.f) 

881160000 

1.74 

93.68 

35.2464 

6294 

1460 

1  87  aaaa  (L.f) 

744285429 

1.47 

95.15 

29.7714 

744286 

45 

rest  ~(REST.f) 

715600000 

1.41 

96.56 

28.6240 

14312 

944 

_l_252_aaaf_  (L.f) 

*  -p[rocedures]  using  invocation  counts; 

*  sorted  in  descending  order  by  number  of  calls  per  procedure; 

*  unexecuted  procedures  are  excluded 


3S3941539  invocations  total 


calls 

%calls 

cum4 

bytes 

procedure  (file) 

226720000 

64.06 

64.06 

112 

r_sign  (r_sign.c) 

79437108 

22.44 

86.50 

248 

sqrtf  (sqrtf.s) 

21720000 

6.14 

92.64 

64 

cvmgp_  ((TVMCSP.f) 

8024000 

2.27 

94.90 

48 

mp  call  (aux.s) 

3640000 

1.03 

95.93 

3240 

_viseul_54_aaaa_  (VISEDL.f) 

3600000 

1.02 

96.95 

2608 

_vj8eul_57_aaaa_  (VJSEOL.f) 

*  -h[eavy]  using  basic-block  counts; 

*  sorted  in  descending  order  by  the  number  of  cycles  executed  in  each 

*  line;  unexecuted  lines  are  excluded 


procedure  (file) 

_lmtr j_151_aaad_  (LMTRJ.f) 


line  bytes  cycles 


%  cum  % 


?  2812  9191735840  18.10  18.10 
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_lmtri  (UfTRI.f) 

54 

2780 

8835036657 

17.40 

35.50 

~fairo«  90~a««c  (PSIROB.f) 

7 

3392 

2955960000 

5.82 

41.32 

fajroa  89  aaac~  <FSJROC.f) 

7 

3332 

2867600000 

5.65 

46.97 

_viaauOO**»I  (VISKDL.f) 

7 

3104 

2796500000 

5.51 

52.48 

r_aiqn  Tr  aiqin.c) 

11 

52 

2493920000 

4.91 

57.39 

_v j8eul_57_aaaa_  (VJSEDL.f) 

7 

2240 

1936800000 

3.81 

61.20 

r  sign  Tr  algn.c) 

10 

44 

1813760000 

3.57 

64.78 

T  227  aaaa  (L.f) 

7 

1804 

1433600000 

2.82 

67.60 

_f8iroa_64_aaab_  (FSIROE.f) 

7 

912 

1412880000 

2.78 

70.38 

~fB jroe_63_aaab  (FSJROE.f) 

7 

868 

1321600000 

2.60 

72.99 

imp  fflaatar  wait  Tor  intarlaavad 
75757 

,1  oop_cc>aipl«t  i  on 

(mp.c) 

691 

40  1311377836 

1  149  aaab  (L.f) 

7 

1776 

1109400000 

2.18 

77.75 

1~87  aaaa  (L.f) 

7 

1460 

881160000 

1.74 

79.49 

1  252  aaaf  (L.f) 

? 

944 

715600000 

1.41 

80.90 

2.58 
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ProooMr  numba*  2; 


Profil*  liatlnq  9«n«rat«d  Mon  Fab  14  14:16: I?  1994  with: 

prof  -pixia  -quit  14  -clock  23  noxzla.axa  nozzla.axa.Addra  iiozzla.axa.Countal2636 


*  -plrocaduraa]  uaing  baaic-block  counta; 

*  aortad  in  daacanding  ordar  by  tha  nunbar  of  cyclaa  exacutad  in  aach 

*  procadura;  unaxacutad  procaduraa  ara  axcludad 


47603084663  cyclaa  (1904.2034  aaaonda  at  25.00  maqahartz) 


cyclaa 

tcyclas 

cum  % 

aaconda 

cyclaa 

/call 

9201032962 

19.33 

19.33 

368.0413 

46006 

8863339687 

18.62 

37.95 

354.5336 

63310 

4761120000 

10.00 

47.95 

190.4448 

21 

3626888840 

7.62 

55.57 

145.0756 

1519 

2867600000 

6.02 

61.59 

114.7040 

14338 

2842420000 

5.97 

67.56 

113.6968 

20303 

2807000000 

5.90 

73.46 

112.2800 

802 

2383600000 

5.01 

78.47 

95.4240 

42 

2268000000 

4.76 

83.23 

90.7200 

630 

1433600000 

3.01 

86.24 

57.3440 

28672 

1338840000 

2.83 

89.10 

34.3336 

4853 

1321600000 

2.78 

91.87 

52.8640 

3304 

1109400000 

2.33 

94.21 

44.3760 

5547 

881160000 

1.83 

96.06 

35.2464 

6294 

715600000 

1.30 

97.56 

28.6240 

14312 

bytaa  procadura  (fila) 

/lina 

2812  _lintr]_151_aaad_  (LMTRJ.f) 

406  _lmtri_148_aaad_  (LMTRI.f) 

23  r_8ign  (r_aign.c) 

688  mp  slava  wait  for  work  (mp_8lava.s) 
3332  _f83roa_89_aaac_  (FSJROE.f) 

3392  _f8iroa_90_aaac_  (FSIROE.f) 

1620  _vi8aul_54_aaaa_  (VISED!.. f) 

6  aqrtf  (sqrtf.s) 

373  _vj88ul_57_aaaa_  (VJSEOL.f) 

1804  _1  227_aaaa_  (L.f) 

912  _f8itoo_64_.aaab_  (FSIROE.f) 

868  _f8 jroa_63~aaab_  (FSJROE.f) 

1776  _l_149_aaab_  (L.f) 

1460  _l_87_aaaa_  (L.f) 

944  _l_252_aaaf_  (L.f) 


*  -pCrocaduraa]  using  invocation  counts; 

*  sortad  in  daacanding  ordar  by  nuaibar  of  calls  par  procadura; 

*  unaxacutad  procaduraa  ara  axcludad 


323364632  invocations  total 


calls 

%calls 

cum% 

bytaa 

procadura  (fila) 

226720000 

69.68 

69.68 

112 

r_8ign  (r_8ign.c) 

56800000 

17.46 

87.14 

248 

sqrtf  (sqrtf.s) 

21300000 

6.55 

93.69 

64 

cvmgp_  (CVMGP.f) 

7884000 

2.42 

96.11 

48 

mp  call  (aux.s) 

3600000 

1.11 

97.22 

2608 

vjsaul  57  aaaa  (VJSEDL.f) 

3500000 

1.08 

98.29 

3240 

visaul  54  aaaa  (VISEDL.f) 

*  -h(aavy]  using  basic-block  counts; 

*  sortad  in  daacanding  ordar  by  tha  numbar 

*  lina;  unaxacutad  linas  ara  axcludad 

of  cyclas 

axacuted  in 

sack 

• 

* 

* 

procadura  (fila) 

lina  bytas 

cyclas 

t 

cum  % 

Imtri  151  aaad  (UfTRJ.  f ) 

? 

2812 

9201032962 

19.33 

19.33 

Imtri  148  aaad  (LMTRI.f) 

54 

2780 

8863339687 

18.62 

37.95 

mp  slava  wait  for  work  (mp  slava. s) 

99 

688 

3626888840 

7.62 

45.57 

fsiroa  89  aaac  (FSJROE.f) 

7 

3332 

2867600000 

6.02 

51.59 

fsiroa  90  aaac  (FSIROE.f) 

7 

3392 

2842420000 

5.97 

57.56 

visaul  54  aaaa  (VISEDL.f) 

7 

3104 

2688000000 

5.65 

63.21 

r  sign  (r  aign.c) 

11 

52 

2493920000 

5.24 

68.44 

vjsaul  57  aaaa  (VJSEDL.f) 

? 

2240 

1936800000 

4.07 

72.51 

r  sign  (r  sign.c) 

10 

44 

1813760000 

3.81 

76.32 

I  227  aaaa  (L.f) 

7 

1804 

1433600000 

3.01 

79.33 

fsiroa  64  aaab  (FSIROE.f) 

7 

912 

1358840000 

2.85 

82.19 

fsjroa  63  aaab  (FSJROE.f) 

7 

868 

1321600000 

2.78 

84.97 

1  149  aaab  (L.f) 

7 

1776 

1109400000 

2.33 

87.30 

1  87  aaaa  (L.f) 

7 

1460 

881160000 

1.85 

89.15 

_l_252_aaaf_  (L.f) 

7 

944 

715600000 

1.50 

90.65 
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Prooeasor  Bomber  3i 


Profil*  liatlnq  9«n«r«t«d  Mon  Fob  14  14:16:26  1994  with: 

prof  -pixio  -quit  19  -clock  25  nozzlo.oxo  nozzlo.oxo.Addra  nozzlo.oxo.Countsl2631 


*  -ptrocoduroo]  uainq  baaic-bloch  counta; 

*  aortod  In  doaeending  ordor  by  tho  nuabor  of  cycloa  oxocutod  in  oach 

*  procoduro;  unoxocutod  procoduroa  are  excluded 


334S4808749  cyclea  (1338.1*23  eeoonda  at  2S.00  neqahertz) 


cyclea  %cyclea 

cun  % 

aeconda 

cyclea 

/call 

bytes 

/line 

procedure  (file) 

1259444S500 
(mp  slave. a) 

37.65 

37.65 

503.7778 

5275 

688  mp  slave  wait  for  work 

284?420000 

8.50 

46.14 

113.6968 

20303 

3392 

_f si roe_90_aaac_  (FS IROE . f ) 

2807000000 

8.39 

54.53 

112.2800 

802 

1620 

_vi8oul_54_aaaa_  (VISEOL.f) 

2708600000 

8.10 

62.63 

108.3440 

13543 

3332 

_f8 jroa_89_aaac_  (FSJROE.f) 

2335200000 

6.98 

69.61 

93.4080 

42 

6 

sqrtf  (sqrtf.s) 

2271800000 

6.79 

76.40 

90.8720 

632 

373 

vjseul  57  aaaa  (VJSEDL.f) 

1433600000 

4.29 

80.69 

57.3440 

28672 

1804 

_l_227_aaae_  (L.f) 

1358840000 

4.06 

84.75 

54.3536 

4853 

912 

_f8iroe_64_aaib_  (FSIROE.f) 

1248800000 

3.73 

88.48 

49.9520 

3122 

868 

fajroo  63~aaab  (FSJROE.f) 

1109400000 

3.32 

91.80 

44.3760 

5547 

1776 

~l_149_aaaS_  (L.f) 

881160000 

2.63 

94.43 

35.2464 

6294 

1460 

1  87  aaaa  (L.f) 

715600000 

2.14 

96.57 

28.6240 

14312 

944 

_l_252_aaaf_  (L.f) 

*  -pirocedurea]  uainq  invocation  counta; 

*  aorted  in  deacendinq  order  by  number  of  cal la  per  procedure; 

*  unexecuted  procedurea  are  excluded 


96164634  invocatione  total 


calls 

%calls 

cum% 

bytea 

55600000 

57.82 

57.82 

248 

20700000 

21.53 

79.34 

64 

7544000 

7.84 

87,19 

48 

3600000 

3.74 

90.93 

2608 

3500000 

3.64 

94.57 

3240 

2388001 

2.48 

97.05 

688 

1028000 

1.07 

98.12 

560 

procedure  (file) 

aqrtf  (aqrtf.a) 
cvmgp_  (CVMGP.f) 
mp  call  (aux.s) 

_vj8oul_51_aaaa_  (VJSEDl.f) 
3yi8oul_54~aaaa~  (VISEOL.f) 

_ nip_8lave~wait~for_wor)t  (mp^slave . s ) 

_ ^mp_do_interleaved_8cfaed  (jnp.c) 


*  -h(eavy]  U8inq  ba8ic-block  counta; 

*  aorted  in  deacendinq  order  by  the  number  of  cyclea  executed  in  each 

*  line;  unexecuted  llnea  are  excluded 


procedure  (file) 

line  bytes 

cycles 

% 

cum  % 

mp_a  1  a ve_wal  t_f or_work  (Bq>_8l  a ve .  a  ) 

99 

688 

12594445500 

37.65 

37.65 

3?airoe_90~aaac_  (FSIROE.f) 

? 

3392 

2842420000 

8.50 

46.14 

_f 8  jroe_89~aaac~  (FS JROE . f ) 

? 

3332 

2708600000 

8.10 

54.24 

_vi8eul_54~aaaa_  (VISEOL , f ) 

? 

3104 

2688000000 

8.03 

62.27 

“vjseul  57~aaaa  (VJSEUL. f) 

? 

2240 

1940600000 

5.80 

68.07 

1  227  aaae  (L.f) 

? 

1804 

1433600000 

4.29 

72.36 

_f8itoe_64_aaab_  (FSIROE . f ) 

7 

912 

1358840000 

4.06 

76.42 

f8jroe~63  aaab  (FSJROE.f) 

? 

868 

1248800000 

3.73 

80.15 

Il_149_aaaE_  (L.f) 

•> 

1776 

1109400000 

3.32 

83.47 

1  87  aaaa  (L.f) 

1460 

881160000 

2.63 

86.10 

~l~252_aaaf_  (L.f) 

7 

944 

715600000 

2.14 

88.24 
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Prooetsor  number  4: 


Profil*  ll«tinq  qaMratad  Hon  Fab  14  14:16:40  1994  Mith: 

prof  -pixia  -quit  19  -clock  25  nozzla.axa  nozzla.axa.Addrs  nozzle. oxa.Countal2638 


*  -p[rocaduraa]  uainq  baaic-block  counts; 

*  sorted  In  dascandinq  order  by  the  nunbar  of  cycles  executed  in  each 

*  procedure;  unexacut^  procedures  are  excluded 


33S21020588  cycles  (1340.8408  seconds  at  25.00  negahertz) 

cycles  %cycles  cum  %  seconds  cycles  bytes  procedure  (file) 

/call  /line 


12807481280 

38.21 

38.21 

512.2993 

5364 

688  mp  slave  wait 

_for_work 

(mp  slave. s) 

2841400000 

8.48 

46.69 

113.7360 

20310 

3392 

_f s 1 roe_90_aaac_ 

(FSIROE. f) 

2807000000 

8.37 

55.06 

112.2800 

802 

1620 

_viseul_54_aaaa_ 

(VISEOL.f) 

2710000000 

8.08 

63.15 

108.4000 

13550 

3332 

_f s  j  roe_8  9_aaac_ 

(FSJROE.f) 

2318400000 

6.92 

70.06 

92.7360 

42 

6 

sqrtf  (sgrtf.s) 

2142000000 

6.39 

76.45 

85.6800 

630 

373 

_v  j8eul_57_aaaa_ 

(VJSEOL.f) 

1433600000 

4.28 

80.73 

57.3440 

28672 

1804 

_l_227_aaae_  (L. 

f) 

1360800000 

4.06 

84.79 

54.4320 

4860 

912 

_f si roe_64_aaab_ 

(FSIROE. f) 

1251600000 

3.73 

88.52 

50.0640 

3129 

868 

_f  s iroe_63_aaab_ 

(FSJROE.f) 

1109400000 

3.31 

91.83 

44.3760 

5547 

1776 

1  149  aaab  (L. 

f) 

881160000 

2.63 

94.46 

35.2464 

6294 

1460 

1  87  aaaa  (L.f) 

715600000 

2.13 

96.60 

28.6240 

14312 

944 

_l_252_aaaf_  (L. 

f) 

*  -p(rocedures]  using  invocation  counts; 

*  sorted  in  descending  order  by  nixnber  of  calls  per  procedure; 

*  unexecuted  procedures  are  excluded 


93364636  invocations  total 


calls 

%call8 

evank 

bytes 

55200000 

57.88 

57.88 

248 

20700000 

21.71 

79.59 

64 

7344000 

7.70 

87.29 

48 

3500000 

3.67 

90.96 

3240 

3400000 

3.57 

94.53 

2608 

2388001 

2.50 

97.03 

688 

1028000 

1.08 

98.11 

560 

procedure  (file) 

sgrtf  (sgrtf.s) 
cv»gp_  (CVMGP.f) 

_ Bp_call  (aux.s) 

viseul  54  aaaa  (VISEOL.f) 
~v3aeol~572aaaa~  (VJSEOL.f) 

mp  slave~ wait  for  work  (mp_slavo . s ) 
mp~ do  interleaved  ached  (mp.c) 


*  -h[eavy]  using  basic-block  counts; 

*  sorted  in  descending  order  by  the  number  of  cycles  executed  in  each 

*  line;  unexecuted  lines  are  excluded 


procedure  (file) 

aip_slave  wait  for  work  (mp_slave.s) 
^si  roe_90~aaac~  (FSIROE .  f ) 
_f8jroe_89  aaac^  (FSJROE.f) 
viseul  S43aaaa_  (VISEaL.f) 
~v3seul~57  aaaa  (VJSEOL.f) 

~1  227  aaae_  (LTf) 

~f8iroe_64_aaab_  (FSIROE. f) 
fsjroe~63  aaab~  (FSJROE.f) 
~l_149_aaa5_  (LTf) 

1~87  aaaa  ~(L.f' 

_055_aaaT_  (L..: 


line  bytes 

cycles 

% 

cum  % 

99 

688 

12807481280 

38.21 

38.21 

7 

3392 

2843400000 

8.48 

46.69 

? 

3332 

2710000000 

8.08 

54.77 

? 

3104 

2688000000 

8.02 

62.79 

•> 

2240 

1829200000 

5.46 

68.25 

7 

1804 

1433600000 

4.28 

72.53 

7 

912 

1360800000 

4.06 

76.59 

7 

868 

1251600000 

3.73 

80.32 

7 

1776 

1109400000 

3.31 

83.63 

7 

1460 

881160000 

2.63 

86.26 

7 

944 

715600000 

2.13 

88.39 
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B  J  FluvMiiation  of  SabroutfaieL/ 


SOBROOTIIK  L 
C  XMCLOOC  'qcoMon.f' 

PAMMETKR  (ILO  •  102,  OLD  -  12,  ILDKX  -  103,  IVHMC  -  103) 
PMtAMETCR  (IPC  •  2) 

PJOtAMETER  (ItDPl  -  103,  JLOPl  *  13) 

PARAMETER(PI-3. 1413926) 

COtMON  /IX:DEL/FLOIIVR(103,13,4,4) 

COIMON/CVGH/XtIOBM,  RCSID,  PRESTO,  TRES  ID,  QRESID,  FRSTRO,  RLIFTJ,  ORAGI 
COWWM/CFLVL/  ICFL,  CFLEICP,  CFUOX,  CFL,  CFLDGl,  CFCRHO,  CFCEI , 

I  CFLPEM, CEXPPEM, IFREEZE, IMOFRZ,  IREST, 

1  ROG2,  TDG2,  OOG2,  VDG2,  COHSTt, R2BR1,CDC2,  C0DG2,  SHKAN,  SLOPE, 

2  lUNST, lONEM, ISDST, ISDEN 

COtftlON  /CONS/RHO(102,12,2),  RHOU<102,  12,  2) ,  RBOV(102, 12,  2) ,  E(102, 
X  12,2) 

COMHOH  /GASO/D(102,12),  V(102,12),  EI(102,72),  P(102,72) 

COMtON  /MSBPT/X(103, 13),  ¥(103,13),  VOL(103,13) 

COIMON/GASP/GAIMA,  GHl,  GGHl, CV, PR,CVISC1,  CVISC2,  IMETRC,  RC 
COMfOM/MESHV/IL,  JL,  ILMl,  JLMl,  IGRTD,  IMSl,  IHS2,  INS3,  INS4, 

$  JtiSl,  J1)S2,  JNS3,  JMS4 

COtMOH  /GTEMP/M,  M,  SX,  SY,  IP3D.  lAOD,  JADD,  FS1V(103, 4, 3) , 

X  MSTART,  HI,  II,  JJ,  CFLTMP,  ICFLAG 
CCIM40H  /GTIME/HAOV,  NEHD,  DT,  TIKE,  OTDG(103, 13) 

(XJtMOH/ALG/IHPLT,  IIXn-ST,  ISNVL,  ILMTR,OHEGA,  DELTEP, 

C  lEHTH, IKS, OELTIL, NSNPS 

COIMON/FLPRM/IAOBHL,  EIMALL,  RL,  OIHF,  VINF,  TINF,  RHOINF,  CINF, 

C  RM,CHARTIH,REL,TliALL,PINF,ICASE,CBARLEH,RMDINF, 

$  ALPHA, CAPaiNF,RCAPOIHF,REIHF, 

$  TINHALL,  AIHLPHA,  RINK,  RINEL,  RIHL,  TIHINF, 

$  EIHIHALL,  CIHIHF,  CIKAPOINF,  OIHIHF,  VINIMF,  RIIMOINF, 

$  R1HROINF,RIHCAPOIIO',PIHIHF,RIHEIHF, 

$  AGGLPHA,  RGGM,  RGGEL,  RGGL,  TGGINF, 

$  CGGIHF,CGGAPOINF,aGGIMF,VGGIHF,RGa(DtHF, 

$  RGGROIHF, RGGCAPOIMF, PGGINF, RGGEIMF, 

$  xmcbr,  ncISQ 

COHMOH  /CHliG/Da(102,12,4),  008(102,12,4) 

COMHOH  /GAOSS/G11(103),  G12(103),  G13(103),  614(103),  G21(103), 

X  622(103),  623(103),  624(103),  631(103),  632(103),  633(103 

X  ),  634(103),  641(103),  642(103),  643(103),  644(103) 

CHARACTER*20  HORBCOH,  GRAHOFL,  TRANF,  HORBLIH,  STATFL,  F6RI0,  FOAT 
X  ,  FMOHT,  FFLHIH,  FFLHOT,  P30GR0,  P30FLM,  CVGFL 
COI«OH/SYS6H/tSYS,  lOHHRB,  I0H6RH,  lOHTIU),  lOHHL,  lOHSTF,  IUH6R0,  IDHDAT, 
C  lUWQIT,  lOHFIH,  lOHFOT,  IOHP3G,  IOHP3F,  lONCVG, 

C  HORBCOH,  CRAHDFL,  TRAHF,  HORBLIH,  STATFL,  FGRIO,  FOAT, 

C  FMOHT,  FFLHIH,  FFLHOT,  P30CR0,P30FLN,  CVGFL 

COtMOH/ lOXX/IREAO, IP300P, lOGBDG, ITRAHP, 

C  MODPR,HRST,IFHRTI, IFMRT0,IIHT1 
IHTEGER  III,  II2,  II3,  II4,  IIS,  II6,  III,  II8,  II9,  IIIO,  Till, 

X  1112,  Ills,  1114 

00  20  M-1,2 
RLIFTJ  -  0. 

ORAGI  -  0. 

M  -  3  -  H 
II2  -  JLMl 
DO  SI  111-2,112 

IF  (ISHVL  .EQ.  1)  THEM 
CALL  STEGI  (III) 

ELSE 

IF  (ISHVL  .EQ.  3)  THEM 
CALL  FIVLVC  (III) 

ELSE 

IF  (ISHVL  .EQ.  4)  THEM 
CALL  FSIROE  (III) 

EHD  IF 
EMO  IF 
EHD  IF 

cmcb_Cig(Mntin6  out  solid  tioundary  fluxes  on  step  (cowl  lip) 

cacb  will  force  the  use  of  the  conditions  specified  in  initl.f  and  bc.f 

C  RESET  FLDXES  OH  STEP 


cmcb  IF  ((J  .GE.  JMSl)  .AND.  (J  .LE.  JNS2-1))  THEN 

cmcb  SX-+(y(INS2,  J+1)-Y(INS2,  J)) 

cmcb  SY— (X(IHS2,  J+1)-X(IHS2,  J)) 

cmcb  PIJJ-P(IMS2,  J) 

cmcb  FSlV(INS2-l,l,l)-0.0 

cmcb  FSlV(IHS2-l,2,l)-Piaj*SX 
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c»cb  rsiv<iiis2-i,3,i)>PXJJ*sy 

amcb  rsiV(lllS2-l,4,l)-0.0 

cMcb  UDIP 

a»eb_«n^K>d 

IF  (IMS  .EQ.  1)  THEN 
CALL  VISEOL  (III) 

ENDtr 

C$MP  SCHEDTYPE-INTERLEAVE 

C$DOSCROSS  SRAR£(ILMl,IIl,DTDG,VOL,Oa.FSlV),LOCAL(DTV,IIll,KOHY,I) 

DO  3  KOMY-1.4 

III!  -  MOD  (IlMl  -  1,  4) 

DO  2  I-2,II11«1 

DTV  *  DTDGd.IIl)  /  VOL(I,  III) 

DO(I,IIl,KOMy)  -  -DTV  •  (FSlVd.KDHy,  1)  -  FSlVd-1, 

X  KOHT,!)) 

2  CONTIMOE 

DO  3  Z-IIll+2, ILM1,4 

DTV  »  DTDGddll)  /  VOLd.lIl) 

DOd.IIl.KOMY)  -  -DTV  •  (FSlVd.KOMY,  1)  -  FSlVd-1, 

X  KDMY.  1)) 

DTV  »  DTDGd+1,111)  /  VOLd+l.IIl) 

OOd+l.IIl.KDMY)  -  -DTV  •  (FS IV (1+1, KDMY,  1)  -  FSlVd, 

X  KDMY,!)) 

DTV  «  DTDG(I+2,II1)  /  VOL(I+2,IIl) 

Da(I+2,  III, KDMY)  -  -DTV  •  (FSIV (1+2, KDMY,  1)  -  FSlVd+1 
X  ,KDHY,1)) 

DTV  -  DTDGd+3,111)  /  VOL(I+3,  Ill) 

DO (1+3,  III, KDMY)  -  -DTV  •  (FS1V(I+3,KDMY,  1)  -  FSlVd+2 
X  ,KDMY,1)) 

3  CONTIMOE 

31  CONTIMOE 

DO  1011  I-2,ILM1 

IF  dSNVL  .EQ.  1)  THEM 
CALL  STEGJd) 

ELSE 

IF  dSHVL  .EQ.  3)  THEN 
CALL  FJVLVCd) 

ELSE 

IF  (ISHVL  .EQ.  4)  THEN 
CALL  FSJRCSd) 

END  IF 
END  IF 
END  IF 

C  FOR  NOZZLE  FLON:  Jal  haa  now  bawn  modiflad  to  b«  a  solid  surfaco 
SX  -  -(Yd+1,2)  -  Y(I,2)) 

SY  -  (Xd+1,2)  -  Xd,2)) 

FS1V(1,1,1)  -  0. 

FS1V(1,2,1)  -  P(I,1)  ♦  SX 
FS1V(1,3,1)  -  P(I,1)  •  SY 
FS1V(1,4,1)  •  0. 

C  Also  for  lip  -  lower  surfac* 

IF  (d  .GE.  IMSl)  .AMD.  (I  +  1  .LE.  INS2) )  THEN 
SX  -  -(Yd+l, JNSl)  -  Yd,JNSl)) 

SY  -  (Xd+1,JNS1)  -  Xd.JNSl)) 

FS1V(JNS1-1,1,1)  -  0. 

FS1V(JNS1-1,2, 1)  -  P(I,JNS1-1)  •  SX 
FS1V(JMS1-1, 3,1)  -  P(I,JNS1-1)  •  SY 
FS1V(JMS1-1,4,1)  -  0. 

C  OPPER  LIP 

SX  -  -(Yd+1,JNS2)  -  Yd,JNS2)) 

SY  •  (X(I+1,JNS2)  -  X(I,JNS2)) 

FS1V(JNS2-1,1,1)  •  0. 

FS1V(JMS2-1,2,1)  -  P(I,JNS2)  •  SX 
FS1V(JNS2-1,3,1)  -  P(I,JNS2>  *  SY 
FS1V(JMS2-1,4,1)  -  0. 

END  IF 

IF  (IMS  .EQ.  1)  THEN 
CALL  VJSEOLd) 

ENDIF 

CSMP  SCHEDTYPE-INTERLEAVE 

C$DOACROSS  SHARE (JLMl, I, DTDG, VOL, DO, FSIV) , LOCAL (DTV, 1112, KDMY,  III ) 

DO  3  KI»1Y-1,4 

1112  -  MOD  (JLMl  -  1,  4) 

DO  4  111-2,1112+1 

DTV  -  DTDG(I,II1)  /  VOLd.lIl) 

DOd,Ill,Kt»lY)  -  DOd, III, KDMY)  -  DTV  •  (FSIV (III, KDMY 
X  ,1)  -  FS1V(II1-1,KDMY,  D) 

4  CONTIMOE 
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00  3  JLM1,4 

DTV  -  DTDGd.IIX)  /  VOL(I,IIl. 

D0(I,II1,K0MY)  •  00(1,  III, XDMY)  -  DTV  *  (rSlV(lIl,KDMY 
X  ,1)  -  FSlV(IIl-l,KDI«,in 

DTV  -  DTDGd,  Xll-flt  /  VOLd,  111*1) 

DO d,  111*1, KDMY)  -  00(1,111+1, KDMY)  -  DTV  •  (FSlVdIl* 
X  l,KDMY,l)  -  FSlVdIl.KDMY,!)) 

DTV  >  DTDGd,  111*2)  /  VOLd,  111*2) 

OOd,  II1*2,K0MY)  -  DOd.  IIl*2,KOt4Y)  -  DTV  •  (FSlVdll* 
X  2,K0My,  1)  -  PS1V(II1*1,K0MY,1)) 

DTV  -  DTDGd,  111*3)  /  VOLd,  111*3) 

D0(I,II1*3,KDMY)  -  DO  (1, 111*3,  KDMY)  -  DTV  *  (FSlVdIl* 
X  3,KWfY,l)  -  FS1V(II1*2,KDMY,1)) 

3  CONTIMOE 

1011  CONTZNOC 

C  FOR  NOZZLE  FLON,  BLANK  OOT  LIP  -  Ellminatad 
in  -  INS2  -  1 

118  -  JNS2  ~  1 
C$MP  SCaEOTYPE>tNTEBLEAVE 

C$DOACROSS  IF(dNS2  --  INSl)  *  (JNS2  -  JNSl)  .GT.  31),SBARE(JNS1, 118,  II7, 
C$<  INSl,00),L0CALdI13,KDMY,IIl,I) 

DO  ^  KDMY-1.4 

DO  7  II1-JNS1,II8 

1113  -  MOD  dI7  -  INSl  *  1,  4) 

DO  6  I>INS1,II13*INS1-1 

DO(I, III, KDMY)  -  0. 

6  COMTINOE 

DO  7  I-II13+INS1,I17,4 

DO(I,  III,  KDMY)  >  0. 

00(1*1, III, KDMY)  -  0. 

D0(I*2, III, KDMY)  >  0. 

D0(I*3, Ill, KDMY)  «  0. 

7  COMTINOE 

IF  (dREST  .EQ.  1)  .AMO.  (H  .EO.  D)  CALL  REST 
IF  (IMPLT  .EQ.  1)  THEM 
CALL  LIMP 
ELSE 

IF  (IMPLT  .ME.  0)  THEN 

PRINT  *,  'NIERD  VALOE  OF  IMPLT  ',  IMPLT 
STOP 
END  IF 
END  IP 

C  FOR  NOZZLE  FLOW,  BLANK  OOT  LIP  -  ElimiA«t«d 

119  -  INS2  -  1 
IIIO  ••  JNS2  -  1 

C$MP  SCHEDTYPE«INTERLEAVE 

CSDOACROSS  IF(dNS2  -  INSl)  *  (JMS2  -  JNSl)  .GT.  31) ,  SHARE  (JNSl,  IIIO,  119 
C$(  ,IMSl,DD),L0CJadI14,KDMY,IIl,I) 

DO  9  KOMr-1,4 

DO  9  II1-JNS1,II10 

1114  -  MOD  (119  -  INSl  *  1,  4) 

DO  8  I>INS1,II14*INS1-1 

DOd,  III, KDMY)  >  0. 

8  CX>MTIMOE 

DO  9  I-II14+INS1, 119,4 

DOd,  III, KDMY)  >  0. 

00(1*1, III, KDMY)  -  0. 

DO (1*2, I II, KDMY)  -  0. 

DO (1*3, III, KDMY)  -  0. 

9  COMTINOE 
RN-M 

RlMl  >  RN  -  1. 

RES ID  -  0. 

IF  (lOGBOG  .EQ.  0)  THEN 
C 

C  AddAd  parallel  construct  by  hand  -  pfa  does  not  )inow  how  to  use 
C  reduction  variables  -  mcb 
C 

C$MP  SCHEDTYPE-INTERLEAVE 

C$DOXCROSS  REDOCTION (RESID) , LOCAL (I, JX, J,RHOOLD,  DRHO,  RUOLD, DRHOO,  RVOLD, 
C$&  DRHOV,BEOLD,DEVL) 

DO  10  I-2,ILM1 
DO  10  JX-2,JLM1 

RHOOLD  -  RHOd,  (JL-JX+1),1) 

RHOd,  (JL-JX+1),M)  =•  (RMd  •  RHOd,  (JL-JX*1)  ,M)  +  RHO( 
X  I,  (JL-JX*1),N)  *  DOd,  (JL-JX*1),1))  /  RN 

DRHO  -  RHOd,  (JL-JX+1),M)  -  RHOOLD 
ROOLD  -  RHOOd,  (JL-JX+1),1) 

RHOOd,  (JL-JX*1),H)  -  (RNMl  •  RHOO (I,  ( JL- JX+1 ) , M)  + 
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X  tlBOOd,  (JL-JX-fl),M)  *  DO<I,  (Jl.-JX4^1),2))  /  RM 
ORflOO  •  RfiOOd.  <JL-JX>1),H)  -  ROOLD 
RVOLD  -  RBOVd,  (Jt-JX't'l),!) 

RBOVd,  •>  <R1M1  *  RBOVd,  (JL-JX+1),M)  * 

X  RHOVd,  (JL-JX-t-D.H)  *  OOd,  (JL-JX-t-DfS))  /  RN 
DRBOV  -  RBOVd,  (JL-JX>1)»M)  -  RVOLD 
REOLD  -  Ed,  (Jt-JX41),l) 

Ed,  (Jt-JX4^1).M)  •  (RMMl  *  Ed,  (JL-JX41),M)  *  Ed,  (JL- 
X  JX+l),H)  +  OOd,  (Jl-JX*l),4))  /  RM 

DEVL  -  Ed,  (JL-JX-fDfM)  -  REOLD 

RESID  -  RESID  *  ((DRBO  /  RHOIMF)  **  2  +  (DRHOU  / 

X  RCAPOIMF)  •*  2  +  (DRBOV  /  RC^paiHT)  •*  2  +  (OEVL  /  REIKF)  *•  2) 

X  /  (DTOGd,  (JL-JX-fl) )  •  DTDGd,  (JL-JX+X))) 

10  (X>I)TIMDE 
C$MP_SCBEDTYPE- INTERLEAVE 

C$DOACROSS  SHARE (ILHl, JLMl , JL, D, M, RHOO, RHO, V, RHOV, El. E, P, Oil ) . LOCAL ( I, JX 
C$t  ) 

DO  11  1-2, ILMl 
DO  11  JX-2,JLm 

Od,  <JL-JX+1) )  -  RHOOd,  <JL-JX+1),M>  /  RHOd,  (JL-JX+1) 
X  ,H) 

Vd,  (JL-JX+1) )  -  RHOVd,  (JL-JX+1  ),M)  /  RHOd,  (JL-JX+1) 
X  M) 

Eld,  (JL-JX+1) )  -  Ed,  (JL-JX+1  ),M)  /  RHOd,  (JL-JX+1), M 
X  )  -  0.3  *  (0(1,  (JL-JX+1 ))  *  Od,  iJL-JX+X))  +  Vd,  (JL-JX+1) )  *  V( 
X  I, (JL-JX+1))) 

Pd,  (JL-JX+D)  -  GMl  *  RHOd,  (JL-JX+1  ),M)  *  Eld,  (JL- 
X  JX+D) 

11  CONTINOE 
ELSE 


C 

C  Conmentad  out  max  accumulation  coda  -  Just  uaed  for  output  and 
C  difficult  to  parallaliza  -  mcb 
C 

Cmcb  RESMX-0.0 

Ctecb  IRSDHX-2 

Cmcb  JRSI)MX-2 

C$MP  SCHEDTYPE-INTERLEAVE 

CSDOACROSS  REOOCTIOM (RESID)  ,  LOCAL  (I,  JX,  J, RB(X)LD,  DRHO, 

C$t  ROOLD,  DRBOO,  RVOLD,  DRBOV,  REOLD,  DEVL,  TERM) 

DO  12  1-2, ILMl 

DO  12  JX-2,JLM1 

RHCXXLD  -  RHOd,  (JL-JX+1),  1) 

RHOd,  (JL-JX+1), M)  -  (RlWl  •  RHOd,  (JL-JX+1), M)  +  RHO( 
X  I,  (JL-JX+1), N)  +  DOd,  (JL-JX+1), 1))  /  RN 

DRHO  -  RHOd,  (JL-JX+1), M)  -  RHOOLD 
ROOLD  -  RHOOd,  (JL-JX+1),  1) 

RHOOd,  (JL-JX+1), M)  -  (RHMl  •  RHOOd,  (JL-JX+1)  ,M)  + 

X  RHOOd,  (JL-JX+1), M)  +  D0(I,  (JL-JX+1), 2))  /  RN 
DRHOO  -  RHOOd,  (JL-JX+1), M)  -  ROOLD 
RVOLD  -  RHOVd,  (JL-JX+1),  1) 

RHOVd,  (JL-JX+1), M)  -  (RlWl  *  RHOVd,  (JL-JX+X) ,M)  + 

X  RHOVd,  (JL-JX+1), N)  +  DOd,  (JL-JX+1),  3))  /  RN 
DRHOV  -  RHOVd,  (JL-JX+1), M)  -  RVOLD 
REOLD  -  Ed,  (JL-JX+1),  1) 

Ed,  (JL-JX+1  ),M)  -  (RIWl  •  Ed,  (JL-JX+1), M)  +  Ed,  (JL- 
X  JX+1),N)  +  DOd,  (JL-JX+1), 4))  /  RN 

DEVL  -  Ed,  (JL-JX+1), H)  -  REOLD 

TERM  -  ((DRHO  /  RHOINF)  •*  2  +  (DRHOO  /  RCAPDINF)  •*  2 
X  +  (DRBOV  /  RCAPDINF)  +*  2  +  (DEVL  /  REIMF)  ••  2)  /  (DTDGd,  (JL- 
X  JX+D)  •  DTDGd,  (JL-JX+1))) 
cmcb  IF  (RESMX  .LE.  TERM)  THEN 

Cmcb  RESMX-TERM 

Qltcb  IRSDMX-I 

Cmcb  JRSDHX-J 

Cmcb  ENDIF 

RESID-RESID+TERH 
12  CONTINOE 

C$MP  SCHEDTYPE-INTERLEAVE 

CSDOACROSS  SHARE ( ILHl ,  JLHl,  JL,  0, M,  RHOO, RHO, V,  RHOV,  El,  E,  P,  GMl ) ,  LOCAL  ( I,  JX 

) 


DO  13  1-2, ILMl 
DO  13  JX-2,JLH1 

0(1,  (JL-JX+D)  -  RHOOd,  (JL-JX+D, M)  /  RHOd,  (JL-JX+1) 
X  ,M) 

V(I,  (JL-JX+D)  -  RHOVd,  (JL-JX+1), M)  /  RHOd,  (JL-JX+1) 
X  ,M) 

Eld,  (JL-JX+D)  =  Ed,  (JL-JX+1  ),M)  /  RHOd,  (JL-JX+1), M 
X  )  -  0.5  *  (0(1,  (JL-JX+D)  *  0(1,  (JL-JX+1))  +  V(I,  (JL-JX+D)  *  V( 
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X 

X 


X 


I. 


(Jl-JX+l))) 

P<I,  (Jl-JX+D) 
JX+D) 
CONTItnC 

PBIMT  *,  'MfcX  RESID 
JRSDHX 


GMl  •  RHO(I, 

SORT  (RESMX).  • OCCURS  AT 


RESID  -  S<»T  (RESID)  /  ((ILMl  -  D  *  tJUll  -  1>) 
CALL  BC 
20  COHnSOE 
return 

END 


El (I. (JL- 
,  IRSDMX, 
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Appendix  C:  Preaure  and  Madi  Number  Contouri 
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Figure  C-1  Pressure  Contours,  Mach  1.9,  Gas  Generattx-  Off 


Experimental  Nozzle/Cowl  Mach  Number  Distributitm 
Mach  L9,  NPRa7.0,  Gas  Oenerator  Off 
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Figure  C-2  Mach  Number  Ckmtours,  Mach  1.9,  Gas  Generator  Off 
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Figure  C-3  Pressure  Contours,  Mach  1.9,  Gas  Generator  100% 
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Figure  C-4  Mach  Number  Contours,  Madi  1.9,  Gas  Generator  100% 


Figure  C-5  Pressure  Contours,  Mach  3.0,  Gas  Generator  Off 
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Figure  C-6  Mach  Number  Contours,  Mach  3.0,  Gas  Generatm*  Off 
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Figure  C-9  Pressure  Contours,  Mach  10,  Gas  (Aerator  Off 


Figure  OlO  Mach  Number  Gmtours,  Mach  10,  Gas  Generator  Off 
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10.0  20.0  30.0  40.0  50.0  60.0  70.0  80.0  90.0  100.0 

Downstream  Location  (in) 


Oeneric  Nozzl«/Cowl  Pressure  Distribution 
Mach  IS.O,  Gas  Generator  Off 
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10.0  20.0  30,0  40.0  50.0  60.0  70,0  80.0  90.0  100.0 

Downstream  Location  (in) 


igure  C*14  Mach  Number  Contours,  Mach  15,  Gas  Generator  Off 
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Figure  C-IS  Pressure  Contours,  Mach  15,  Gas  Generata- 100% 
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Figure  C-17  Pressure  Contours,  Mach  20,  Gas  Generator  Off 


igure  C-18  Mach  Number  Contours,  Mach  20,  Gas  Generator  Off 
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'igure  C-20  Mach  Number  G)titours,  Madi  20.  Gas  Geno'atw  100% 
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